Sorry, I have the answer: the Hessian is mass-weighted after, line 2066.
Many thanks for your help,
Marie
> On Tue, Oct 13, 2009, Brut Marie wrote:
>>
>> I am using NAB to compute an Hessian matrix with gb=1.
>> I first minimized my structure (with conjgrad and newton functions).
>> Until
>> here, no problem. Then I thought I could use mme2 alone to get the
>> Hessian matrix, using the command :
>>
>> float mme2( float x[], float g[], float h[], float mass[], int iter );
>>
>>
>> mme2( x, g, h, m, 1 );
>
> OK: the manual is out of date, and new arguments have been added to the
> mme2()
> routine. See the source code itself (sff2.c, probably around line 7131,
> depending on which version you have). If you haven't turned on SCALAPACK
> (which I have never done, myself), I think you can set its variables to
> dummies.
>
> But I'd recommend modifying nmode.c, and printing out the hessian after
> it is computed there (say at about line 1884).
>
> ...good luck....dac
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
--
Marie Brut
LAAS - CNRS
N2IS group
7 Avenue du Colonel Roche
31 077 Toulouse Cedex 4
Tél. : (33) 05 61 33 62 47
Fax : (33) 05 61 33 62 08
Mail : mbrut.laas.fr
http://www.laas.fr/laas/1-4267-Annuaire.php
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Oct 13 2009 - 12:30:02 PDT
