[CKH]

Dear foamers,

This might be more of a C++ question...

I've got a function that is supposed to return List that is made up of 2 scalarLists, i.e. scalarListList dsf(2) whereby dsf[0] = scalarList dsf0 and dsf[1] = scalarList dsf1

My code extract is as follows:

scalarListList testMesh::calcDSF
(
    const scalar xi_ ,
    const scalar eta_
) const
{
    Info<<"Calculating shape functions derivatives"<<endl;

    scalarList dSfdXi(8);
    scalarList dSfdEta(8);
    scalarListList dsf(2);

    // calculate shape func derivative wrt to Xi: dSfdXi
    dSfdXi[0] =  ....................
    ........
    dSfdXi[7] = ........................;

    // calculate shape func derivative wrt to eta: dSfdEta
    dSfdEta[0] = ...............;
    ..............
    dSfdEta[7] = ............;

    dsf[0] = dSfdXi;
    dsf[1] = dSfdEta;

    return dsf;
}

The above code compiles but stops halfway when run.....I believe it has to do with the allocation of memory because by tweaking the order of scalarList / scalarListList initialization, the code runs to different points.

My objective is to create a 2 x 8 matrix where the 1st column stores the dSfdXi values and the 2nd column the dSfdEta values (for calculating jacobian, transforming to parametric coordinates).

Any ideas or alternative implementation (instead of using OF's Lists) would be much appreciated!

Thanks!

ckh

[CKH]

Doing a more thorough search through the forums, I believe an autoPtr might be needed, i.e.

List< autoPtr<scalarList> > dsf(2);

dsf[0].setSize(8);
dsf[0][0] = /* some code to calculate each of the dSfdXi...
...
dsf[0][7] = .... */

dsf[1].setSize(8);
dsf[1][0] = /* some more code  to calculate each of the dSfdEta...
...
dsf[1][7] = .... */

Any ideas / comments would be much appreciated!!

Thanks in advance!

ckh