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

Attachment: 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

Reply via email to