On Fri, 1 Nov 2013, Manav Bhatia wrote:
I think I understand the purpose of that routine. I seems to be an
approximation to the following form:
[J] x_weighted = - sum_i ( w_i {partial R / partial p_i} )
and is only called from the routines dealing with the Hessian in ImplicitSystem.
It's been way too long since I wrote the weighted_sensitivity_solve()
code - the purpose was to enable Hessian-vector products without
requiring full Hessian calculations. It's grossly untested, so if
there's something that looks fishy I'd suspect the code is at fault
more than your understanding of the code.
In fact, I don't think I ever got around to writing an app using this
myself, so I've *never* properly hit it with a benchmark. Vikram, is
there any chance you did?
---
Roy
------------------------------------------------------------------------------
Android is increasing in popularity, but the open development platform that
developers love is also attractive to malware creators. Download this white
paper to learn more about secure code signing practices that can help keep
Android apps secure.
http://pubads.g.doubleclick.net/gampad/clk?id=65839951&iu=/4140/ostg.clktrk
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
