[gmx-users] Normal Mode Analysis

Franzi Hoffgaard f.hoffgaard at gmail.com
Wed May 27 10:44:16 CEST 2009


Hi,

I don't know exactly what went wrong when I tried to patch the file. There
was a > on line 18 in the fix file, actually. I checked it. With the file
you sent me, gmxdump works fine.

Thank you very much for your help.

Best regards, Franzi

---
Franziska Hoffgaard
PhD Student
Bioinformatics & Theo. Biology Group
TU Darmstadt


2009/5/26 Berk Hess <gmx3 at hotmail.com>

>  Hi,
>
> But there is a '>' on line 18, at least in my mail.
> I tested the patch and it worked for me.
> Maybe something went wrong when cut and pasting the patch from the email?
> I can mail you the whole file personally tomorrow, if required.
>
> Berk
>
> ------------------------------
> Date: Tue, 26 May 2009 14:47:52 +0200
> Subject: Re: [gmx-users] Normal Mode Analysis
> From: f.hoffgaard at gmail.com
> To: gmx-users at gromacs.org
>
> Hi
>
> thanks for your reply. I tried patching as you described it, but I got the
> error message:
>
> patch: **** `>' expected at line 18 of patch
>
> Franzi
>
> ---
> Franziska Hoffgaard
> PhD Student
> Bioinformatics & Theo. Biology Group
> TU Darmstadt
>
> 2009/5/26 Berk Hess <gmx3 at hotmail.com>
>
>  Hi,
>
> I just saw that 4.0 uses a new mtx format which can also efficiently store
> sparse matrices.
> Please try if my modified gmxdump for 4.0.5 works?
> Store the data below in a file called fix and then do:
> patch gmxdump.c fix
>
> Berk
>
> 61a62,63
> > #include "sparsematrix.h"
> > #include "mtxio.h"
> 364a367,402
> > static void list_mtx(char *fn)
> > {
> >   int  nrow,ncol,i,j,k;
> >   real *full=NULL,value;
> >   gmx_sparsematrix_t * sparse=NULL;
> >
> >   gmx_mtxio_read(fn,&nrow,&ncol,&full,&sparse);
> >
> >   if (full == NULL) {
> >     snew(full,nrow*ncol);
> >     for(i=0;i<nrow*ncol;i++) {
> >       full[i] = 0;
> >     }
> >
> >     for(i=0;i<sparse->nrow;i++) {
> >       for(j=0;j<sparse->ndata[i];j++) {
> >         k     = sparse->data[i][j].col;
> >         value = sparse->data[i][j].value;
> >         full[i*ncol+k] = value;
> >         full[k*ncol+i] = value;
> >       }
> >     }
> >     gmx_sparsematrix_destroy(sparse);
> >   }
> >
> >   printf("%d %d\n",nrow,ncol);
> >   for(i=0; i<nrow; i++) {
> >     for(j=0; j<ncol; j++) {
> >       printf(" %g",full[i*ncol+j]);
> >     }
> >     printf("\n");
> >   }
> >
> >   sfree(full);
> > }
> >
> 381a420
> >     { efMTX, "-mtx", "hessian", ffOPTRD },
> 409a449,450
> >   else if (ftp2bSet(efMTX,NFILE,fnm))
> >     list_mtx(ftp2fn(efMTX,NFILE,fnm));
>
>
> ------------------------------
> Express yourself instantly with MSN Messenger! MSN Messenger<http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/>
>
> _______________________________________________
> gmx-users mailing list    gmx-users at gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please search the archive at http://www.gromacs.org/search before posting!
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
>
>
>
> ------------------------------
> Express yourself instantly with MSN Messenger! MSN Messenger<http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/>
>
> _______________________________________________
> gmx-users mailing list    gmx-users at gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please search the archive at http://www.gromacs.org/search before posting!
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20090527/7fb177c1/attachment.html>


More information about the gromacs.org_gmx-users mailing list