[Trilinos-Users] evalModel, dfdp with several parameters

Bartlett, Roscoe A rabartl at sandia.gov
Thu May 27 11:13:57 MDT 2010


Nico,

From: Nico Schlömer [mailto:nico.schloemer at ua.ac.be]
Sent: Thursday, May 27, 2010 11:06 AM
To: Bartlett, Roscoe A
Cc: trilinos-users at software.sandia.gov
Subject: Re: [Trilinos-Users] evalModel, dfdp with several parameters



This just means that if p(0) has five components, then if you say that you support dfdp(0), that means that you have to support the computation of the derivatives for all five of the scalar parameters in p(0).
-- and I still have to explicitly set them by



outArgs.setSupports( OUT_ARG_DfDp, 0, DerivativeSupport(DERIV_MV_BY_COL) );

outArgs.setSupports( OUT_ARG_DfDp, 1, DerivativeSupport(DERIV_MV_BY_COL) );

outArgs.setSupports( OUT_ARG_DfDp, 2, DerivativeSupport(DERIV_MV_BY_COL) );

outArgs.setSupports( OUT_ARG_DfDp, 3, DerivativeSupport(DERIV_MV_BY_COL) );

outArgs.setSupports( OUT_ARG_DfDp, 4, DerivativeSupport(DERIV_MV_BY_COL) );



and fetch them with



RCP<Epetra_MultiVector> dfdp0_out = outArgs.get_DfDp(0).getMultiVector();

RCP<Epetra_MultiVector> dfdp1_out = outArgs.get_DfDp(1).getMultiVector();

RCP<Epetra_MultiVector> dfdp2_out = outArgs.get_DfDp(2).getMultiVector();

RCP<Epetra_MultiVector> dfdp3_out = outArgs.get_DfDp(3).getMultiVector();

RCP<Epetra_MultiVector> dfdp4_out = outArgs.get_DfDp(4).getMultiVector();



unless I want Thyra to complain. Right?



[Bartlett, Roscoe A]



Typically, you would aggregate a bunch of scalar parameters into a single parameter vector like p(0) and then just support a single MultiVector dfdp(0).  Otherwise, if you support Np parameter *vectors* and you want to support derivatives w.r.t. to all of these, you will need to support Np derivative muti-vectors.



- Ross


-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://software.sandia.gov/pipermail/trilinos-users/attachments/20100527/4d655e1b/attachment.html 


More information about the Trilinos-Users mailing list