On Wed, 28 Mar 2012, David Andrs wrote: > 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?
You're using the system properly; that was just an awful typo (or brain failure? had I been using Matlab and indexing from 1 recently?) in WrappedFunction. It should be fixed in the svn head now; thanks very much for tracking it down. --- Roy ------------------------------------------------------------------------------ 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
