void fun()
{
try {
fxx();
} catch(const std::exception& e) {
std::string s = e.what();
assert(0);
}
}
void funcL(TInt x)
{
fxx();
}
void func(TInt x)
{
fxx();
}
int operator+(TInt x,TInt y)
{
fxx();
}
void TypeOfFE_P2ttdc::FB(const bool *whatd,const Mesh & ,const Triangle & K,const R2 & P1,RNMK_ & val) const
{
R2 P=Shrink1(P1);
// const Triangle & K(FE.T);
R2 A(K[0]), B(K[1]),C(K[2]);
R l0=1-P.x-P.y,l1=P.x,l2=P.y;
R l4_0=(4*l0-1),l4_1=(4*l1-1),l4_2=(4*l2-1);
// throwassert(FE.N == 1);
throwassert( val.N()>=6);
throwassert(val.M()==1);
// throwassert(val.K()==3 );
val=0;
// --
if (whatd[op_id])
{
RN_ f0(val('.',0,op_id));
f0[0] = l0*(2*l0-1);
f0[1] = l1*(2*l1-1);
f0[2] = l2*(2*l2-1);
f0[3] = 4*l1*l2; // oppose au sommet 0
f0[4] = 4*l0*l2; // oppose au sommet 1
f0[5] = 4*l1*l0; // oppose au sommet 3
}
if( whatd[op_dx] || whatd[op_dy] || whatd[op_dxx] || whatd[op_dyy] || whatd[op_dxy])
{
R2 Dl0(K.H(0)*cshrink1), Dl1(K.H(1)*cshrink1), Dl2(K.H(2)*cshrink1);
if (whatd[op_dx])
{
RN_ f0x(val('.',0,op_dx));
f0x[0] = Dl0.x*l4_0;
f0x[1] = Dl1.x*l4_1;
f0x[2] = Dl2.x*l4_2;
f0x[3] = 4*(Dl1.x*l2 + Dl2.x*l1) ;
f0x[4] = 4*(Dl2.x*l0 + Dl0.x*l2) ;
f0x[5] = 4*(Dl0.x*l1 + Dl1.x*l0) ;
}
if (whatd[op_dy])
{
RN_ f0y(val('.',0,op_dy));
f0y[0] = Dl0.y*l4_0;
f0y[1] = Dl1.y*l4_1;
f0y[2] = Dl2.y*l4_2;
f0y[3] = 4*(Dl1.y*l2 + Dl2.y*l1) ;
f0y[4] = 4*(Dl2.y*l0 + Dl0.y*l2) ;
f0y[5] = 4*(Dl0.y*l1 + Dl1.y*l0) ;
}
if (whatd[op_dxx])
{
RN_ fxx(val('.',0,op_dxx));
fxx[0] = 4*Dl0.x*Dl0.x;
fxx[1] = 4*Dl1.x*Dl1.x;
fxx[2] = 4*Dl2.x*Dl2.x;
fxx[3] = 8*Dl1.x*Dl2.x;
fxx[4] = 8*Dl0.x*Dl2.x;
fxx[5] = 8*Dl0.x*Dl1.x;
}
if (whatd[op_dyy])
{
RN_ fyy(val('.',0,op_dyy));
fyy[0] = 4*Dl0.y*Dl0.y;
fyy[1] = 4*Dl1.y*Dl1.y;
fyy[2] = 4*Dl2.y*Dl2.y;
fyy[3] = 8*Dl1.y*Dl2.y;
fyy[4] = 8*Dl0.y*Dl2.y;
fyy[5] = 8*Dl0.y*Dl1.y;
}
if (whatd[op_dxy])
{
assert(val.K()>op_dxy);
RN_ fxy(val('.',0,op_dxy));
fxy[0] = 4*Dl0.x*Dl0.y;
fxy[1] = 4*Dl1.x*Dl1.y;
fxy[2] = 4*Dl2.x*Dl2.y;
fxy[3] = 4*(Dl1.x*Dl2.y + Dl1.y*Dl2.x);
fxy[4] = 4*(Dl0.x*Dl2.y + Dl0.y*Dl2.x);
fxy[5] = 4*(Dl0.x*Dl1.y + Dl0.y*Dl1.x);
}
}
}
int operator+(TInt x)
{
fxx();
}
