On Tue, 23 Feb 2010, Yujie wrote:
> Without mesh refinement, how to solve the linear equation with multiple
> rhses. In current libmesh framework,
>
> one can use "equation_systems.get_system("system").solve()" to get the
> solution with one rhs.
>
> In .solve(), there are assemble() and solver() functions. rhs can be
> replaced in assemble(). However, which parameters need to be reinitialized
> after finishing one rhs? Whether it is ok to call equation_systems.reinit()
> to reinitialize necessary parameters for multiple rhses? Thanks a lot.
Depends on how complicated your system is, but in general you
shouldn't need to reset anything except (perhaps) the solution (which
gets used as a starting guess in iterative solvers) before solving
with a new right hand side. Take a look at the parameter sensitivity
methods in implicit_system.C for some examples; since we use multiple
solution and rhs vectors there we don't have to reinitialize anything
at all.
---
Roy
------------------------------------------------------------------------------
Download Intel® Parallel Studio Eval
Try the new software tools for yourself. Speed compiling, find bugs
proactively, and fine-tune applications for parallel performance.
See why Intel Parallel Studio got high marks during beta.
http://p.sf.net/sfu/intel-sw-dev
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
