hdiff output

r15481/config.h 2017-01-21 10:41:16.578824024 +0000 r15480/config.h 2017-01-21 10:41:24.331157048 +0000
 30:  30: 
 31: #  Set the C compiler, etc.  31: #  Set the C compiler, etc. 
 32:  32: 
 33: #          For GNU:  CC-->gcc; LEX-->flex; YACC-->bison -y -t; 33: #          For GNU:  CC-->gcc; LEX-->flex; YACC-->bison -y -t;
 34: #          Note: If your lexer is "really" flex, you need to set 34: #          Note: If your lexer is "really" flex, you need to set
 35: #          LEX=flex below.  For example, on many linux distributions, 35: #          LEX=flex below.  For example, on many linux distributions,
 36: #          /usr/bin/lex is really just a pointer to /usr/bin/flex, 36: #          /usr/bin/lex is really just a pointer to /usr/bin/flex,
 37: #          so LEX=flex is necessary.  In general, gcc seems to need 37: #          so LEX=flex is necessary.  In general, gcc seems to need
 38: #          flex. 38: #          flex.
 39:  39: 
 40: CC=gcc   40: CC=gcc 
 41: ##CFLAGS=-march=pentium4  41: ##CFLAGS=-march=pentium4 
 42: #CFLAGS=-O3 -Mnoframe -Mflushz -Munroll 42: #CFLAGS=-O3 -Mnoframe -Mflushz -Munroll
 43: #OCFLAGS=-O3 -Mnoframe -Mflushz -Munroll 43: #OCFLAGS=-O3 -Mnoframe -Mflushz -Munroll
 44: CFLAGS=-O0 -g  44: 
 45: ##-Mprefetch=w -tp=amd64e   45: ##-Mprefetch=w -tp=amd64e  
 46: ##CFLAGS=-O3 -Mprefetch=w -tp=amd64e 46: ##CFLAGS=-O3 -Mprefetch=w -tp=amd64e
 47: #NABFLAGS= 47: #NABFLAGS=
 48: ##CFLAGS= -march=pentium4  48: ##CFLAGS= -march=pentium4 
 49: ##OCFLAGS=-O3 -march=pentium4   49: ##OCFLAGS=-O3 -march=pentium4  
 50: ##NABFLAGS= 50: ##NABFLAGS=
 51: CFLAGS = $(EXTRA_CFLAGS)  51: CFLAGS = $(EXTRA_CFLAGS) 
 52:  52: 
 53: LEX=    flex 53: LEX=    flex
 54: YACC=bison -y -t 54: YACC=bison -y -t


r15481/prm.c 2017-01-21 10:41:23.175107971 +0000 r15480/prm.c 2017-01-21 10:41:30.891431833 +0000
401: int mme2wrapper_(REAL_T *x,REAL_T *ener,REAL_T *vnew,401: int mme2wrapper_(REAL_T *x,REAL_T *ener,REAL_T *vnew,
402:         REAL_T *temphess,REAL_T *atmass,REAL_T *grad1)402:         REAL_T *temphess,REAL_T *atmass,REAL_T *grad1)
403: 403: 
404: {404: {
405:  405:  
406: /*static MOLECULE_T *mol;*/406: /*static MOLECULE_T *mol;*/
407: /* MOLECULE_T *moltemp;*/407: /* MOLECULE_T *moltemp;*/
408: /*static PARMSTRUCT_T *prm;*/408: /*static PARMSTRUCT_T *prm;*/
409:  REAL_T *x0i;409:  REAL_T *x0i;
410:  int temp;410:  int temp;
411:  int i; 
412:  INT_T descF_PxQ; INT_T descF_1x1; INT_T descG_PxQ; INT_T context_PxQ; 411:  INT_T descF_PxQ; INT_T descF_1x1; INT_T descG_PxQ; INT_T context_PxQ; 
413:  INT_T context_1x1; INT_T context_Nx1;412:  INT_T context_1x1; INT_T context_Nx1;
414:  INT_T gridim;413:  INT_T gridim;
415:  int iteration;414:  int iteration;
416:  iteration = 1;415:  iteration = 1;
417:  context_PxQ = 0;416:  context_PxQ = 0;
418:  context_1x1 = 0;417:  context_1x1 = 0;
419:  context_Nx1 = 0;418:  context_Nx1 = 0;
420:  descF_PxQ = 0;419:  descF_PxQ = 0;
421:  descF_1x1 = 0;420:  descF_1x1 = 0;
430: /* iteration = 1;429: /* iteration = 1;
431:  mme_(ener,x,vnew,&iteration);*/430:  mme_(ener,x,vnew,&iteration);*/
432: /*431: /*
433:  iteration = -3;432:  iteration = -3;
434:  mme2(ener,x,vnew,temphess,atmass,grad1,&descF_PxQ,&descF_1x1,&descG_PxQ,433:  mme2(ener,x,vnew,temphess,atmass,grad1,&descF_PxQ,&descF_1x1,&descG_PxQ,
435:               &context_PxQ,&context_1x1,&context_Nx1,&gridim,&mol->m_prm->Natom,&iteration);*/434:               &context_PxQ,&context_1x1,&context_Nx1,&gridim,&mol->m_prm->Natom,&iteration);*/
436:  iteration = 1;435:  iteration = 1;
437:  mme2(ener,x,vnew,temphess,atmass,grad1,&descF_PxQ,&descF_1x1,&descG_PxQ,436:  mme2(ener,x,vnew,temphess,atmass,grad1,&descF_PxQ,&descF_1x1,&descG_PxQ,
438:               &context_PxQ,&context_1x1,&context_Nx1,&gridim,&mol->m_prm->Natom,&iteration);437:               &context_PxQ,&context_1x1,&context_Nx1,&gridim,&mol->m_prm->Natom,&iteration);
439: 438: 
440: /*439: 
441: for(i=0;i<9*prm->Natom*prm->Natom;i++){ 
442:         printf("iii,g[i] %d, %lf\n",i,temphess[i]); 
443: } 
444: */ 
445: /* mol = moltemp;*/440: /* mol = moltemp;*/
446: /* mme_init(mol, NULL, NULL, NULL, NULL);*/441: /* mme_init(mol, NULL, NULL, NULL, NULL);*/
447: 442: 
448: return (0);443: return (0);
449: }444: }
450: 445: 
451: 446: 
452: 447: 
453: 448: 
454: /*449: /*


r15481/sff2.c 2017-01-21 10:41:23.419118341 +0000 r15480/sff2.c 2017-01-21 10:41:31.219445398 +0000
1297: 1297: 
1298: 1298: 
1299:     ct2 = ct*ct;1299:     ct2 = ct*ct;
1300:     st2 = 1.0-ct2;1300:     st2 = 1.0-ct2;
1301:     st = sqrt(st2)*s/fabs(s) ;1301:     st = sqrt(st2)*s/fabs(s) ;
1302: 1302: 
1303:     if (st >= 0 && st < 1.e-3)1303:     if (st >= 0 && st < 1.e-3)
1304:       st = 1.e-3;1304:       st = 1.e-3;
1305:     else if (st < 0 && st > -1.e-3)1305:     else if (st < 0 && st > -1.e-3)
1306:       st = -1.e-3;1306:       st = -1.e-3;
1307: 1307:     if (st2 >= 0 && st2 < 1.e-3)
1308:       st2 = st * st; 
1309: /*    if (st2 >= 0 && st2 < 1.e-3) 
1310:       st2 = 1.e-3;1308:       st2 = 1.e-3;
1311:     else if (st2 < 0 && st2 > -1.e-3)1309:     else if (st2 < 0 && st2 > -1.e-3)
1312:       st2 = -1.e-3;1310:       st2 = -1.e-3;
1313: */1311: 
1314:   multi_term:1312:   multi_term:
1315:     iper = (int) fabs(Pn[atyp]) + 0.0001;1313:     iper = (int) fabs(Pn[atyp]) + 0.0001;
1316:     assert(iper != (fabs(Pn[atyp]) - 0.0001));1314:     assert(iper != (fabs(Pn[atyp]) - 0.0001));
1317:     assert(iper >= 1 && iper <= 6);1315:     assert(iper >= 1 && iper <= 6);
1318:     switch (iper) {1316:     switch (iper) {
1319: 1317: 
1320:     case 1:1318:     case 1:
1321:       e = ct;1319:       e = ct;
1322:       df = 1.0;1320:       df = 1.0;
1323:       ddf = 0.0;1321:       ddf = 0.0;


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0