Output_Data::p2p1x1_1D Class Reference

#include <export.h>

Public Member Functions
	p2p1x1_1D (const Grid_Info &_G)
	p2p1x1_1D (const p2p1x1_1D &other)
	~p2p1x1_1D ()
Array2D< float >	operator() (DistFunc1D &df, size_t x0, size_t s)
size_t	Species () const
size_t	Nl (size_t s) const
size_t	Nm (size_t s) const
size_t	Npx (size_t s) const
size_t	Npy (size_t s) const
float	Pmin (size_t s) const
float	Pmax (size_t s) const
PLegendre2D	PL (size_t s) const
Array2D< float >	p2p1_x1 (size_t s) const
Array2D< float >	prad (size_t s) const
Array2D< size_t >	npcell (size_t s) const
Array2D< float >	dcell (size_t s) const

Private Attributes
vector< float >	pmin
vector< float >	pmax
vector< size_t >	numl
vector< size_t >	numm
vector< PLegendre2D >	Pl
vector< Array2D< float > >	pr
vector< Array2D< size_t > >	nextpcell
vector< Array2D< float > >	distancetothatcell
vector< Array2D< float > >	p2p1x1

Detailed Description

Definition at line 327 of file export.h.

Constructor & Destructor Documentation

p2p1x1_1D() [1/2]

Output_Data::p2p1x1_1D::p2p1x1_1D

(

const Grid_Info &

_G

)

Definition at line 1156 of file export.cpp.

References Grid_Info::axis, Grid_Info::l0, Grid_Info::m0, Algorithms::MakeAxis(), Algorithms::AxisBundle< T >::Np(), Algorithms::AxisBundle< T >::Npx(), Algorithms::AxisBundle< T >::Npy(), Algorithms::AxisBundle< T >::pdim(), Algorithms::AxisBundle< T >::pmax(), and Algorithms::AxisBundle< T >::pmin().

                                                    {

//  Generate the required structures
    for (size_t s(0); s < _G.axis.pdim(); ++s) {
        pmax.push_back( static_cast<float>(_G.axis.pmax(s)) );
        pmin.push_back( static_cast<float>(_G.axis.pmin(s)) );
        numl.push_back( _G.l0[s]);
        numm.push_back( _G.m0[s]);

        Pl.push_back( PLegendre2D( _G.l0[s], _G.m0[s], _G.axis.Np(s), static_cast<float>(_G.axis.pmin(s)),
                                   static_cast<float>(_G.axis.pmax(s)), _G.axis.Npx(s), _G.axis.Npy(s) ) );
        pr.push_back( Array2D<float>( _G.axis.Npx(s),_G.axis.Npy(s)));


        nextpcell.push_back( Array2D<size_t>( _G.axis.Npx(s),_G.axis.Npy(s)));
        distancetothatcell.push_back( Array2D<float>( _G.axis.Npx(s),_G.axis.Npy(s)));


        p2p1x1.push_back( Array2D<float>( _G.axis.Npx(s),_G.axis.Npy(s)) );
    }

    size_t k(0);
    //  Calculate a 2D array for the pr for each px and py   
    for (size_t s(0); s < _G.axis.pdim(); ++s) {

        valarray<float> p( Algorithms::MakeAxis( pmin[s], pmax[s], _G.axis.Np(s) )),
                px(Algorithms::MakeAxis( float(-1.0)* pmax[s], pmax[s], _G.axis.Npx(s) )),
                py(Algorithms::MakeAxis( float(-1.0)* pmax[s], pmax[s], _G.axis.Npy(s) ));
//      --->
        for (size_t i(0); i < _G.axis.Npx(s); ++i){
            for (size_t j(0); j < _G.axis.Npy(s); ++j){
                pr[s](i,j) = sqrt(px[i]*px[i]+py[j]*py[j]);

                while ((pr[s](i,j) > p[k]) && (k < _G.axis.Np(s))) ++k;
                if (k == 0) ++k;
                nextpcell[s](i,j) = k-1;
                distancetothatcell[s](i,j) = (p[k] - pr[s](i,j))/(p[1]-p[0]);
                // std::cout << "\n\n (px,py,pr,p[k],nextpcell(i,j)) = " << px[i]
                // << ", " << py[j] 
                // << ", " << pr[s](i,j) 
                // << ", " << p[k] 
                // << ", " << nextpcell[s](i,j);

                k = 0;

            }
        }
//      <---
    }


//    f0x1 = new valarray<float>(Npx);
}

Here is the call graph for this function:

p2p1x1_1D() [2/2]

Output_Data::p2p1x1_1D::p2p1x1_1D

(

const p2p1x1_1D &

other

)

Definition at line 1211 of file export.cpp.

References dcell(), Nl(), Nm(), npcell(), p2p1_x1(), PL(), Pmax(), Pmin(), prad(), and Species().

                                                       {

    for (size_t s(0); s < other.Species(); ++s) {
        pmin.push_back( other.Pmin(s) );
        pmax.push_back( other.Pmax(s) );
        numl.push_back( other.Nl(s));
        numm.push_back( other.Nm(s));
        Pl.push_back(other.PL(s));
        pr.push_back(other.prad(s));
        nextpcell.push_back(other.npcell(s));
        distancetothatcell.push_back(other.dcell(s));
        p2p1x1.push_back( other.p2p1_x1(s) );
    }
}

Here is the call graph for this function:

~p2p1x1_1D()

Output_Data::p2p1x1_1D::~p2p1x1_1D

(

)

Definition at line 1227 of file export.cpp.

                                {

}

Member Function Documentation

dcell()

Array2D<float> Output_Data::p2p1x1_1D::dcell ( size_t  s ) const

inline

Definition at line 354 of file export.h.

Referenced by p2p1x1_1D().

{return distancetothatcell[s];}

Here is the caller graph for this function:

Nl()

size_t Output_Data::p2p1x1_1D::Nl ( size_t  s ) const

inline

Definition at line 341 of file export.h.

Referenced by p2p1x1_1D().

{ return numl[s]; }

Here is the caller graph for this function:

Nm()

size_t Output_Data::p2p1x1_1D::Nm ( size_t  s ) const

inline

Definition at line 342 of file export.h.

Referenced by p2p1x1_1D().

{ return numm[s]; }

Here is the caller graph for this function:

npcell()

Array2D<size_t> Output_Data::p2p1x1_1D::npcell ( size_t  s ) const

inline

Definition at line 353 of file export.h.

Referenced by p2p1x1_1D().

{return nextpcell[s];}

Here is the caller graph for this function:

Npx()

size_t Output_Data::p2p1x1_1D::Npx ( size_t  s ) const

inline

Definition at line 343 of file export.h.

{ return Pl[s](0).dim1(); }

Npy()

size_t Output_Data::p2p1x1_1D::Npy ( size_t  s ) const

inline

Definition at line 344 of file export.h.

{ return Pl[s](0).dim2(); }

operator()()

Array2D< float > Output_Data::p2p1x1_1D::operator()	(	DistFunc1D &	df,
		size_t	x0,
		size_t	s
	)

Definition at line 1234 of file export.cpp.

                                                                                    {

    p2p1x1[s] = 0.0; // this is a Array2D<float>
    // std::cout << "\n\n x0 = " << x0 << " \n";
    size_t i(0), temploc(0);
    float tempdist(0.0);

    // valarray<float> p( Algorithms::MakeAxis(pmin, pmax, Np));

    for (size_t l(0); l < Nl(s)+1; ++l){
        for (size_t m(0); m<((Nm(s)<l)?Nm(s):l)+1; ++m){

            i = ((l < Nm(s)+1)?((l*(l+1))/2+m):(l*(Nm(s)+1)-(Nm(s)*(Nm(s)+1))/2 + m));

            for (size_t ipy(0); ipy < Npy(s); ++ipy) {
                for (size_t ipx(0); ipx < Npx(s); ++ipx) {
                    //  At each location in px,py
                    if (nextpcell[s](ipx,ipy) == 0)
                    {
                        temploc = 1;
                        tempdist = 1.0;
                    }
                    else
                    {
                        temploc = nextpcell[s](ipx,ipy);
                        tempdist = distancetothatcell[s](ipx,ipy);
                    }

                    p2p1x1[s](ipx,ipy) += Pl[s](i)(ipx, ipy) *
                                          ((float((df(i))(temploc  ,x0).real())*(1.0-tempdist))+
                                           (float((df(i))(temploc-1,x0).real())*(    tempdist)));
                }
            }
        }
    }
    return p2p1x1[s];
}

p2p1_x1()

Array2D<float> Output_Data::p2p1x1_1D::p2p1_x1 ( size_t  s ) const

inline

Definition at line 351 of file export.h.

Referenced by p2p1x1_1D().

{ return p2p1x1[s]; }

Here is the caller graph for this function:

PL()

PLegendre2D Output_Data::p2p1x1_1D::PL ( size_t  s ) const

inline

Definition at line 350 of file export.h.

Referenced by p2p1x1_1D().

{ return Pl[s]; }

Here is the caller graph for this function:

Pmax()

float Output_Data::p2p1x1_1D::Pmax ( size_t  s ) const

inline

Definition at line 347 of file export.h.

Referenced by p2p1x1_1D().

{ return pmax[s]; }

Here is the caller graph for this function:

Pmin()

float Output_Data::p2p1x1_1D::Pmin ( size_t  s ) const

inline

Definition at line 346 of file export.h.

Referenced by p2p1x1_1D().

{ return pmin[s]; }

Here is the caller graph for this function:

prad()

Array2D<float> Output_Data::p2p1x1_1D::prad ( size_t  s ) const

inline

Definition at line 352 of file export.h.

Referenced by p2p1x1_1D().

{ return pr[s]; }

Here is the caller graph for this function:

Species()

size_t Output_Data::p2p1x1_1D::Species ( ) const

inline

Definition at line 340 of file export.h.

Referenced by p2p1x1_1D().

{ return p2p1x1.size(); }

Here is the caller graph for this function:

Field Documentation

distancetothatcell

vector< Array2D<float> > Output_Data::p2p1x1_1D::distancetothatcell

private

Definition at line 362 of file export.h.

nextpcell

vector< Array2D<size_t> > Output_Data::p2p1x1_1D::nextpcell

private

Definition at line 361 of file export.h.

numl

vector<size_t> Output_Data::p2p1x1_1D::numl

private

Definition at line 358 of file export.h.

numm

vector<size_t> Output_Data::p2p1x1_1D::numm

private

Definition at line 358 of file export.h.

p2p1x1

vector< Array2D<float> > Output_Data::p2p1x1_1D::p2p1x1

private

Definition at line 364 of file export.h.

Pl

vector< PLegendre2D > Output_Data::p2p1x1_1D::Pl

private

Definition at line 359 of file export.h.

pmax

vector<float> Output_Data::p2p1x1_1D::pmax

private

Definition at line 357 of file export.h.

pmin

vector<float> Output_Data::p2p1x1_1D::pmin

private

Definition at line 357 of file export.h.

pr

vector< Array2D<float> > Output_Data::p2p1x1_1D::pr

private

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