#ifndef LINSTO_H
#define LINSTO_H

#include "pSTO.h"

using namespace std ;

// linear combination of STO's***********************************************
class linSTO : public vector<pSTO> {
public:

#if 0
	explicit linSTO(int nolin=0) : vector<pSTO>(nolin) ;
#else
	linSTO() ;
#endif

	linSTO(const pSTO & origSTO) ;

	linSTO(const vector<pSTO> & origSTO) ;

	linSTO (const char  *infile) ;

#if 0
	~linSTO() ;
#endif
	linSTO & operator+= ( const vector<pSTO> & toadd) ;

	linSTO operator-() const ;

	linSTO & operator-= ( const linSTO & toadd) ;

	linSTO & operator*= (const complex<double> & fact) ;

#if 0
	linSTO & operator/= (const complex<double> & fact) ;
#endif

	linSTO & operator*= (const pSTO & right) ;

#if 0
	linSTO dup(bool altern) ;
#endif

	linSTO & operator*= (const linSTO & right)  ;

	void shiftOrb(const int p);
private:
} ;

linSTO conj(const linSTO & arg) ;

ostream & operator<<(ostream &os, const linSTO & someSTO) ;

linSTO operator+(const linSTO & left, const linSTO & right) ;

linSTO operator-(const linSTO & left, const linSTO & right) ;

linSTO operator*(const complex<double> & fact, const linSTO & in) ;

linSTO operator*(const double & fact, const linSTO & in) ;

linSTO operator*(const pSTO & left, const linSTO & right) ;

#if 0
linSTO operator/(const linSTO &in, const complex<double> & fact) ;

linSTO operator/(const linSTO &in, const double & fact) ;
#endif

linSTO operator*(const linSTO & left, const linSTO & right) ;

complex<double> TlinSTO(const linSTO & left, const linSTO & right) ;

complex<double> NlinSTO(const linSTO &left, const linSTO &right) ;

complex<double> ClinSTO(const linSTO & left, const linSTO & right) ;

#endif /* LINSTO_H */