Output_Data::PLegendre2D Class Reference

#include <export.h>

Public Member Functions
	PLegendre2D (size_t Nl, size_t Nm, size_t Np, float pmin, float pmax, size_t Npx, size_t Npy)
	PLegendre2D (const PLegendre2D &other)
	~PLegendre2D ()
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 266 of file export.h.

Constructor & Destructor Documentation

PLegendre2D() [1/2]

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

Definition at line 886 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)),
            py(Algorithms::MakeAxis(float(-1.0)*pmax, pmax, Npy));

//  Generate the structure to save the polynomials
    sz = ((Nm+1)*(2*Nl-Nm+2))/2;
//  
    plegendre = new vector< Array2D<float>>(sz,Array2D<float>(Npx,Npy))  ;
    //    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);
    // float invpmax(1.0/pmax);
//  Generate the polynomial values for each cos(theta) = px/p 
    for (size_t j(0); j < py.size(); ++j) {

        for (size_t i(0); i < px.size(); ++i) {
            float invp(sqrt(py[j]*py[j]+px[i]+px[i]));
            if (invp > 0.0 ) invp=1.0/invp;
            else if (invp > pmax) invp *= 2;
            else invp = 0.0;
//          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);
                    // if (invp == 0.0) (*plegendre)[k](i,j) = 0.0;

                    // 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:

PLegendre2D() [2/2]

Output_Data::PLegendre2D::PLegendre2D

(

const PLegendre2D &

other

)

Definition at line 944 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:

~PLegendre2D()

Output_Data::PLegendre2D::~PLegendre2D

(

)

Definition at line 954 of file export.cpp.

                                    {
    delete plegendre;
}

Member Function Documentation

dim()

size_t Output_Data::PLegendre2D::dim ( ) const

inline

Definition at line 275 of file export.h.

Referenced by PLegendre2D().

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

Here is the caller graph for this function:

operator()() [1/2]

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

inline

Definition at line 276 of file export.h.

{ return (*plegendre)[i]; }

operator()() [2/2]

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

inline

Definition at line 277 of file export.h.

{ return (*plegendre)[i]; }   

Field Documentation

lmax

size_t Output_Data::PLegendre2D::lmax

private

Definition at line 281 of file export.h.

plegendre

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

private

Definition at line 280 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