Welcome! Log In Create A New Profile

Advanced

How to write a fieldScalar with Palabos?

Posted by mfduquedmfduqued  
How to write a fieldScalar with Palabos?
September 19, 2017 12:49AM
Dear community,

I am in this moment write a class for calculation of strain strength as

template<typename T, template<typename U> class Descriptor>
auto_ptr<MultiScalarField2D<T> > computeShearstrength(MultiBlockLattice2D<T,DESCRIPTOR>& lattice,
Box2D domain, plint component)
{
std::auto_ptr<MultiScalarField2D<T> > density = *computeDensity(lattice);
std::auto_ptr<MultiTensorField2D<T,3> > shear = *computeStrainRate(*computeVelocity(lattice));
std::auto_ptr<MultiScalarField2D<T> > result = extractComponent(*shear, domain, component);
std::auto_ptr<MultiScalarField2D<T> > strain = multiply((T)0.5, *multiply(*density,*result) );
}


but when I like to write this fieldScalar as a matrix It doesn't work,

successiveStrainUp << setprecision(16) <<
computeShearstrength(lattice, profileTapUp, 1)
<< endl;


I dont understand how to write a FieldScalar with palabos, can someone help me with this aspect.

Thanks for your time.



Edited 1 time(s). Last edit at 09/19/2017 01:33AM by mfduqued.
Sorry, you do not have permission to post/reply in this forum.