Commit 476eb362 authored by COULAUD Olivier's avatar COULAUD Olivier

Fix buf for K_{a,b}(x,y) in FInterpMatrixKernelRH

parent df89a39a
...@@ -121,20 +121,21 @@ struct FInterpMatrixKernelRH :FInterpMatrixKernelR{ ...@@ -121,20 +121,21 @@ struct FInterpMatrixKernelRH :FInterpMatrixKernelR{
static const unsigned int NPOT = 1; //< dim of potentials static const unsigned int NPOT = 1; //< dim of potentials
static const unsigned int NRHS = 1; //< dim of mult exp static const unsigned int NRHS = 1; //< dim of mult exp
static const unsigned int NLHS = 1; //< dim of loc exp static const unsigned int NLHS = 1; //< dim of loc exp
FReal LX,LY,LZ ;
FReal LX,LY,LZ ;
FInterpMatrixKernelRH(const FReal = 0.0, const unsigned int = 0) : LX(1.0),LY(1.0),LZ(1.0) FInterpMatrixKernelRH(const FReal = 0.0, const unsigned int = 0) : LX(1.0),LY(1.0),LZ(1.0)
{ } { }
FReal evaluate(const FPoint& x, const FPoint& y) const FReal evaluate(const FPoint& x, const FPoint& y) const
{ {
const FPoint xy(x-y); const FPoint xy(x-y);
std::cout << "Lx " <<LX <<std::endl ;
return FReal(1.) / FMath::Sqrt(LX*xy.getX()*xy.getX() + return FReal(1.) / FMath::Sqrt(LX*xy.getX()*xy.getX() +
LY*xy.getY()*xy.getY() + LY*xy.getY()*xy.getY() +
LZ*xy.getZ()*xy.getZ()); LZ*xy.getZ()*xy.getZ());
} }
void setCoeff(const FReal& a, const FReal& b, const FReal& c) void setCoeff(const FReal& a, const FReal& b, const FReal& c)
{LX= a ; LY = b ; LZ = c ;} {LX= a*a ; LY = b*b ; LZ = c *c;}
// returns position in reduced storage // returns position in reduced storage
int getPosition(const unsigned int) const int getPosition(const unsigned int) const
{return 0;} {return 0;}
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment