Re: [Libmesh-users] user defined routine for parameter sensitivity

2013-11-06 Thread Manav Bhatia
Hi, I have a first-cut on the modifications to accommodate user provided sensitivity methods. I pushed my local repository to github, so there might be some irrelevant modifications that you can ignore: https://github.com/manavbhatia/libmesh/tree/user_sensitivity_function I will be ad

Re: [Libmesh-users] user defined routine for parameter sensitivity

2013-11-01 Thread Roy Stogner
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

Re: [Libmesh-users] user defined routine for parameter sensitivity

2013-11-01 Thread Manav Bhatia
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. -Manav On Fri, Nov 1, 2013 at 11:49 AM, Manav

Re: [Libmesh-users] user defined routine for parameter sensitivity

2013-11-01 Thread Manav Bhatia
Hi, I am currently working on adding support for use supplied sensitivity. In the process, I looked at the code in ImplicitSystem::weighted_sensitivity_solve(), and now want to make sure I understand what is going on inside it. The header file describes the idea to be [J] x_weighte

Re: [Libmesh-users] user defined routine for parameter sensitivity

2013-10-31 Thread Manav Bhatia
Thanks Paul! I need this soon and will start hacking tomorrow morning. Will keep you in the loop. Manav > On Oct 31, 2013, at 11:09 PM, "Paul T. Bauman" wrote: > >> On Thu, Oct 31, 2013 at 5:00 PM, Roy Stogner >> wrote: >> >> On Thu, 31 Oct 2013, Manav Bhatia wrote: >> >> >My appli

Re: [Libmesh-users] user defined routine for parameter sensitivity

2013-10-31 Thread Paul T. Bauman
On Thu, Oct 31, 2013 at 5:00 PM, Roy Stogner wrote: > > On Thu, 31 Oct 2013, Manav Bhatia wrote: > > >My applications provide analytical sensitivity of RHS for sensitivity > > analysis. The current code structure seems to assumed finite differencing > > for this, and I am hoping to use my rout

Re: [Libmesh-users] user defined routine for parameter sensitivity

2013-10-31 Thread Roy Stogner
On Thu, 31 Oct 2013, Manav Bhatia wrote: >My applications provide analytical sensitivity of RHS for sensitivity > analysis. The current code structure seems to assumed finite differencing > for this, and I am hoping to use my routines instead. > >Perhaps it would be best to provide for a

[Libmesh-users] user defined routine for parameter sensitivity

2013-10-31 Thread Manav Bhatia
Hi, My applications provide analytical sensitivity of RHS for sensitivity analysis. The current code structure seems to assumed finite differencing for this, and I am hoping to use my routines instead. Perhaps it would be best to provide for a user_sensitivity_function, or user_sensitivit