Hi,
Using spatial coordinates (q_point) in matrix Ke almost as it's
presented in the example 14 for matrix Fe, I may misunderstand
something important. Here is the code.
// ------------------------------------------------------------------ //
const std::vector<Point>& qrule_points = fe->get_xyz();
for (unsigned int qp=0; qp<qrule.n_points(); qp++)
{
Real S = qrule_points[qp](0);
// Real S = 10.0;
for (unsigned int i=0; i<phi.size(); i++)
{
for (unsigned int j=0; j<phi.size(); j++)
{
Ke(i,j) += JxW[qp]*(dphi[i][qp]*dphi[j][qp]) * S;
}
}
}
// ------------------------------------------------------------------ //
When the Real S (line 4) is a constant, the equation solution is the
same as Matlab's results whereas it's different when S depends on the
first spatial factor (line 3). What's the good way ?
Best regards,
Maxime
-------------------------------------------------------------------------
This SF.net email is sponsored by: Microsoft
Defy all challenges. Microsoft(R) Visual Studio 2008.
http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
