[Trilinos-Users] access to the Jacobian postprocess

[Veltz Romain]

Dear John,

Thank you for you help.

I am using this Method because the Jacobian cannot be stored in memory (in fact I have also an implementation of the Jacobian using an Operator).

My question is rather "from where can I call the jacobian method?"…sorry for the lack of precision.

Everything I need is in the LOCA::Epetra::Group that I pass to the stepper. Should I implement a new version of the Stepper? Is there a simpler way?

Romain

On Feb 28, 2013, at 6:20 PM, John T. Foster wrote:

> Can you use the accessor method for NOX::Epetra::FiniteDifference
> 
> getUnderlyingMatrix()  ??
> 
> http://trilinos.sandia.gov/packages/docs/r10.0/packages/nox/doc/html/classNOX_1_1Epetra_1_1FiniteDifference.html#3b86ded520525685db9a2eb2ed1cc189
> 
> 
> 
> John T. Foster
> Assistant Professor of Mechanical Engineering
> The University of Texas at San Antonio
> AET 2.330
> One UTSA Circle
> San Antonio, TX 78249
> 210.458.5521
> john.foster at utsa.edu
> 
> 
> On Thu, Feb 28, 2013 at 6:22 PM, Veltz Romain <romain.veltz at inria.fr> wrote:
>> Hello,
>> 
>> I am performing numerical continuation using LOCA::Epetra and a Matrix Free
>> method for the Jacobian (using NOX::Epetra::FiniteDifference).
>> I would like to compute the following scalar product after each successful
>> continuation step: <e1,Jac^(-1).e1> where Jac is the jacobian and e1 is a
>> constant vector (for all the continuation processes).
>> 
>> Does anyone have some advice to do so please?
>> 
>> Thank you a lot,
>> 
>> Romain VELTZ.
>> 
>> _______________________________________________
>> Trilinos-Users mailing list
>> Trilinos-Users at software.sandia.gov
>> http://software.sandia.gov/mailman/listinfo/trilinos-users
>> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://software.sandia.gov/pipermail/trilinos-users/attachments/20130228/3a7421cc/attachment.html