Output_Data::PLegendre1D Class Reference

#include <export.h>

Public Member Functions
	PLegendre1D (size_t Nl, size_t Nm, size_t Np, float pmin, float pmax, size_t Npx)
	PLegendre1D (const PLegendre1D &other)
	~PLegendre1D ()
size_t	dim () const
Array2D< float > &	operator() (size_t i)
Array2D< float >	operator() (size_t i) const

Private Attributes
vector< Array2D< float > > *	plegendre
size_t	lmax

Detailed Description

Definition at line 241 of file export.h.

Constructor & Destructor Documentation

PLegendre1D() [1/2]

Output_Data::PLegendre1D::PLegendre1D	(	size_t	Nl,
		size_t	Nm,
		size_t	Np,
		float	pmin,
		float	pmax,
		size_t	Npx
	)

Definition at line 816 of file export.cpp.

References Algorithms::Legendre(), and Algorithms::MakeAxis().

                                                                            {

    size_t sz(0);

    valarray<float> p( Algorithms::MakeAxis(pmin, pmax, Np)),
            px(Algorithms::MakeAxis(float(-1.0)*pmax, pmax, Npx));

//  Generate the structure to save the polynomials
    sz = ((Nm+1)*(2*Nl-Nm+2))/2;
//  
    plegendre = new vector< Array2D<float>>(sz,Array2D<float>(Npx,Np))  ;
    //    for (size_t m(0); m < Nm+1; ++m) {
    //     for (size_t l(m); l < Nl+1; ++l) { 
    //         (*plegendre).push_back( Array2D<float>(Npx,Np) );
    //         ++sz;
    //  }
    // }
    size_t k(0);
//  Generate the polynomial values for each cos(theta) = px/p 
    for (size_t j(0); j < p.size(); ++j) {
        float invp(1.0/p[j]);
        for (size_t i(0); i < px.size(); ++i) {

//          For given px/p generate all the polynomial values ...
//             valarray<float> vL( Algorithms::Legendre( px[i]*invp, Nl, Nm ) );
            Array2D<float> vL( Algorithms::Legendre( px[i]*invp, Nl, Nm ) );
            // std::cout << "\n x = " << px[i]*invp << ", sqrt(1-x^2) = " << sqrt(1.0-px[i]*invp*px[i]*invp) ;
//          ... and save them
            // std::cout << "dim = " << vL.dim() << "\n";
            // std::cout << "dim1:" << vL.dim1() << "\n";
            // std::cout << "dim2:" << vL.dim1() << "\n\n";
            for (size_t l(0); l < Nl+1; ++l){
                for (size_t m=0; m<((Nm<l)?Nm:l)+1; ++m){
                    // for (size_t k(0); k < sz; ++k) {
                    k = ((l < Nm+1)?((l*(l+1))/2+m):(l*(Nm+1)-(Nm*(Nm+1))/2 + m));

                    (*plegendre)[k](i,j) = vL(l,m);

                    std::cout << "\n LP(" << l << "," << m << "," << k << ") = " <<  vL(l,m);

                }
                // (*plegendre)[k](i,j) = vL(k);
            }

            // exit(1);
        }
    }

}

Here is the call graph for this function:

PLegendre1D() [2/2]

Output_Data::PLegendre1D::PLegendre1D

(

const PLegendre1D &

other

)

Definition at line 868 of file export.cpp.

References dim().

                                                              {

//  Generate the structure to save the polynomials
    plegendre = new vector< Array2D<float> > ;
    for (size_t i(0); i < other.dim(); ++i) {
        (*plegendre).push_back( other(i) );
    }
}

Here is the call graph for this function:

~PLegendre1D()

Output_Data::PLegendre1D::~PLegendre1D

(

)

Definition at line 878 of file export.cpp.

                                    {
    delete plegendre;
}

Member Function Documentation

dim()

size_t Output_Data::PLegendre1D::dim ( ) const

inline

Definition at line 250 of file export.h.

Referenced by PLegendre1D().

{ return (*plegendre).size(); }

Here is the caller graph for this function:

operator()() [1/2]

Array2D<float>& Output_Data::PLegendre1D::operator() ( size_t  i )

inline

Definition at line 251 of file export.h.

{ return (*plegendre)[i]; }

operator()() [2/2]

Array2D<float> Output_Data::PLegendre1D::operator() ( size_t  i ) const

inline

Definition at line 252 of file export.h.

{ return (*plegendre)[i]; }   

Field Documentation

lmax

size_t Output_Data::PLegendre1D::lmax

private

Definition at line 256 of file export.h.

plegendre

vector< Array2D<float> >* Output_Data::PLegendre1D::plegendre

private

Definition at line 255 of file export.h.

---

The documentation for this class was generated from the following files:

• /Users/archis/Dropbox/work/dev/oshun/oshun-OS/OSHUN/1d_cpp/source/export.h
• /Users/archis/Dropbox/work/dev/oshun/oshun-OS/OSHUN/1d_cpp/source/export.cpp