#ifndef PSTO_H
#define PSTO_H
#include <vector>
#include <complex>

#include "Ylm.h"

using namespace std ;

#define NMAX 10

//primitive STO at origin, Nnm*r^(n-1)exp(-alpha*r)Y_l^m ******************************************
class pSTO : public Ylm {
public:
	int n ;		// n=0,1,2,3,.... <= NMAX
	double alpha;		// exponent
	double Nnm;		// normalization (extra factor besides the factor c of the base class)

#if 0
	pSTO() : Ylm() ;
#endif

	explicit pSTO(const int nn, const int ll, const int mm, const double expo, const complex<double> coef=1.)  ;

	pSTO & operator*= (const double & fact) ;

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

	complex<double> NpSTO() const ;

	static double normf(const double alpha, const int n) ;

	void shiftOrb(const int p) ;
private:
	const static double Nntab[NMAX+1] ;

} ;  /* pSTO */

pSTO conj(const pSTO & arg) ;

bool operator== (const pSTO & first, const pSTO & secnd) ;

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

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

vector<pSTO> operator*(const pSTO & left, const pSTO & right) ;

double radiInt(const int n, const double alpha) ;

complex<double> OpSTO(const pSTO & left, const pSTO &right) ;

complex<double> TpSTO(const pSTO & left, const pSTO &right) ;

complex<double> NlinSTO(const vector<pSTO> & in) ;

complex<double> NpSTO(const pSTO & left, const pSTO & right) ;

complex<double> CpSTO(const pSTO &left, const pSTO &right) ;

#endif /* PSTO_H */