I am using DifferentiableSystem class for some time now, together with
EulerSolver class. Recently EulerSolver::element_residual() function
was changed and now it calls _system.mass_residual() only once, instead
of two
times, as in 0.6.2 version.
This has a side-effect that EulerSolver can now solve only equations
linear in time
derivative, that is, of the form
dx/dt = F(x),
while before one could also solve
d/dt M(x) = F(x),
with nonlinear M. The reason is that now the function forms residual
R = dt F(x_theta) + M(x_old - x)
while before it would form:
R = dt F(x_theta) + M(x_old) - M(x).
Is there any reason for this change except optimization? I would like to
be able to treat nonlinearity in time derivative term.
Regards,
Mladen Jurak
<http://web.math.hr/%7Ejurak>
------------------------------------------------------------------------------
SF.Net email is Sponsored by MIX09, March 18-20, 2009 in Las Vegas, Nevada.
The future of the web can't happen without you. Join us at MIX09 to help
pave the way to the Next Web now. Learn more and register at
http://ad.doubleclick.net/clk;208669438;13503038;i?http://2009.visitmix.com/
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
