Hi! We are using scalar variables for solving a system of ODEs along with a system of PDEs. In that ODE system we have a scalar variable of higher order to represent unknown variable \lambda_i. We set an initial condition for this variable to let's say 5000. Now, all the component of the variable are set properly except the first one. If the variable is of a first order, its value is ok.
I tracked the issue down to the one line, see attached patch. It is in WrappedFunction, if the variable has more components the else branch goes over all but the first component. If that is correct, are we not using the system properly? Thanks, -- David Andrs
0001-Fixing-initial-condition-setup-for-scalar-variables-.patch
Description: Binary data
------------------------------------------------------------------------------ This SF email is sponsosred by: Try Windows Azure free for 90 days Click Here http://p.sf.net/sfu/sfd2d-msazure
_______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel