[sharonyue]

Hi guys,

This is my class:

class definition:

namespace Foam
{

/*---------------------------------------------------------------------------*\
			  Class gaussianDistribution Declaration
\*---------------------------------------------------------------------------*/

class gaussianDistribution
{
private:

    // Private data

        //- Parameter sigma
        double sigma_;
	
	//- Moments from previous value
	const scalarDiagonalMatrix moments_;
	
	//- Quadrature weights
	scalarDiagonalMatrix weights_;

	//- Quadrature abscissas
	scalarDiagonalMatrix abscissas_;

	//- Number of weights and abscissas
	const label nWeights_;

	//- Number of moments
	const label nMoments_;

    //- Private member functions

	//- Disallow default bitwise copy construct
	gaussianDistribution(const gaussianDistribution&);

	//- Disallow default bitwise assignment
	void operator=(const gaussianDistribution&);

public:

    // Constructors

	//- Construct from a scalarSquareMatrix
	gaussianDistribution
	(
	    const double& sigma, 
	    const scalarDiagonalMatrix& moments,
	    const label nWeights
	);
};


// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

}

Foam::gaussianDistribution::gaussianDistribution
(
    const double& sigma, 
    const scalarDiagonalMatrix& moments,
    const label nWeights
) 
:
    sigma_(sigma),
    moments_(moments),
    weights_(nWeights),
    abscissas_(nWeights),
    nWeights_(nWeights),
    nMoments_(nWeights*2)  
{
    scalarDiagonalMatrix we(nWeights_, scalar(0));
    scalarDiagonalMatrix abs(nWeights_, scalar(0));
    scalarSquareMatrix A(nMoments_, nMoments_, scalar(0));

    // Construct A matrix
    
    A[2][0] = Foam::pow(sigma_,2);
    A[3][1] = Foam::pow(sigma_,2);
    
    Info << A <<endl;
}

This is my main function:

int main(int argc, char *argv[])
    {
	scalarDiagonalMatrix m(4, scalar(1));
	Info << m << endl;
	gaussianDistribution(4.0, m, 1);
	return 0;
    }

You can see that I just wanna see if this class works. But it said:

a@a:~/OpenFOAM/myCase/momentsinversion$ C
4{1}
Segmentation fault (core dumped)
a@a:~/OpenFOAM/myCase/momentsinversion$

Any ideas? Thanks

Solved..gaussianDistribution(4.0, m, 4) works.

[huangxianbei]

Hi:
Could you please explain a little bit about how to use the new namespace when it's compiled with wmake libso?

Xianbei