/*
 * Programmieren fuer Physiker, SS 2010
 *
 * Vektorklasse: dynamisch anlegen, Index-Bereich ueberpruefen.
 * Destruktor. Copy Constructor. Assignment.
 */

#include <iostream>
using namespace std ;

// ---- CLASS VEC
class Vec
{
private:
    int size ;   // Groese
    double* p ;  // Zeiger auf den Vektor 
public:
    Vec(int) ;  // Konstruktor
    ~Vec() ;    // Destruktor
    Vec(const Vec& v) ;              // "copy constructor"
    Vec& operator= (const Vec& v) ;  // Zuweisung
 
    double & operator[] (int) ;      // Elementzugriff
    // alternativ auch "operator()" moeglich

    // hier sollten nun noch viele Rechenoperationen 
    // fuer Vektoren definiert sein...
} ;

Vec::Vec(int s)  // Konstruktor
{
    size = s ; // nimm die Groesse auf
    p = new double[s] ; // besorge Speicher fuer die Vektorkomponenten
}

Vec::~Vec()  // Destruktor
{
    delete[] p ;  // Speicher freigeben
}

Vec::Vec(const Vec& v)  // copy constructor
// wird bei "Call-by-Value" und Funktionsrueckgaben aufgerufen
{
    size = v.size ;
    p = new double[v.size] ;
    for (int i=0; i<size; i++) p[i] = v.p[i] ;
}

Vec& Vec::operator= (const Vec& v)  // Zuweisung (assignment constructor)
{
    if (this == &v) return *this ;  // Vorsicht bei Selbstzuweisung!
    if (size!=v.size)
	cout << "Zuweisung verschieden grosser Vektoren!!!" << endl ;
    delete[] p ;
    size = v.size ;
    p = new double[v.size] ;
    for (int i=0; i<size; i++) p[i] = v.p[i] ;
    return *this ;
}

double & Vec::operator[](int i)  // Elementzugriff
{
    if (i<0 || i>=size)
    {
	cout << "Illegaler Vektorzugriff. Stop." << endl ;
	exit(1) ;
    }
    return p[i] ;
}

// ---- MAIN
int main()
{
    int k=17 ;
    Vec v(k) ;
    // v[18]=4.4 ;  // gibt Fehler, weil ausserhalb Bereichs
    for(int i=0; i<k; i++)
	v[i] = i ;
    for(int i=0; i<k; i++)
	cout << v[i] << " " ;
    cout << endl ;

    Vec v1(8) ;
    Vec v2(8) ;
    v2[3] = 5 ;
    v1 = v2 ;  // was passiert hier? // ist das gewuenscht?
    cout << v1[3] << " " << v2[3] << endl ;
    v2[3] = 8 ;
    cout << v1[3] << " " << v2[3] << endl ;
    // Frage: Wie kriegt man anderes Verhalten von "=" ?
    // (Antwort: siehe oben: eigener "=" Operator)
}