Dear all,
I want to use the ContinuationSystem class. I noted the following comments:
enum libMesh::ContinuationSystem::RHS_Mode [protected]
There are (so far) two different vectors which may be assembled using the
assembly routine: 1.) The Residual = the normal PDE weighted residual 2.)
G_Lambda = the derivative wrt the parameter lambda of the PDE weighted
residual
It is up to the derived class to handle writing separate assembly code for
the different cases. Usually something like:
switch (rhs_mode) { case Residual: { Fu(i) += ... // normal PDE residual
break; }
case G_Lambda: { Fu(i) += ... // derivative wrt control parameter break; }
My question is that how should I pass the enum rhs_mode to my own assembly
function if I write the following code:
int main (int argc, char** argv)
{
//main program
}
void assemble_jacobian_residual (EquationSystems& es, const std::string&
system_name)
{
assemble the jacobian and residual;
//need some code like this
switch (rhs_mode) {
case Residual: { Fu(i) += ... // normal PDE residual break; }
case G_Lambda: { Fu(i) += ... // derivative wrt control parameter break;}
}
Thank you.
Best wishes,
yunfei
------------------------------------------------------------------------------
Special Offer-- Download ArcSight Logger for FREE (a $49 USD value)!
Finally, a world-class log management solution at an even better price-free!
Download using promo code Free_Logger_4_Dev2Dev. Offer expires
February 28th, so secure your free ArcSight Logger TODAY!
http://p.sf.net/sfu/arcsight-sfd2d
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
