On Thu, 15 Oct 2009, Roy Stogner wrote: > I'm going to need to change out the adjoints API a bit, mostly to > accomodate the possibility of multiple QoIs.
> Give assemble_qoi and assemble_qoi_derivative optional arguments of > type vector<int> telling them which QoIs to assemble "this time" - > typically the argument would just be the default empty vector, which > would be shorthand for "assemble them all". > > The qoi_parameter_sensitivity call (which would now fork off to > forward and adjoint versions depending on which is more efficient, by > comparing qoi->size() and parameters->size()) would now fill a > vector<vector<Number> > of sensitivities. A couple other changes: I'm also going to give that optional indices argument to qoi_parameter_sensitivity and its forward/adjoint versions. Users can already solve with a subset of their systems' parameters by leaving those out of the parameters argument, but the most natural way to solve for a subset of QoIs seems to be with that indices arg. The indices argument is also going in element_qoi, etc. for FEMSystem users. --- Roy ------------------------------------------------------------------------------ Come build with us! The BlackBerry(R) Developer Conference in SF, CA is the only developer event you need to attend this year. Jumpstart your developing skills, take BlackBerry mobile applications to market and stay ahead of the curve. Join us from November 9 - 12, 2009. Register now! http://p.sf.net/sfu/devconference _______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
