hdiff output

r24824/sff2.c 2017-01-21 10:41:15.870793159 +0000 r24823/sff2.c 2017-01-21 10:41:16.114803809 +0000
 22:  * gradient and Hessian are distributed, the ptr1d, ptr1f, ptr2d and ptr2f 22:  * gradient and Hessian are distributed, the ptr1d, ptr1f, ptr2d and ptr2f
 23:  * functions are used to write the subvector and submatrix elements by only 23:  * functions are used to write the subvector and submatrix elements by only
 24:  * the process that owns the subvector or submatrix.  The SUPERD_T type is 24:  * the process that owns the subvector or submatrix.  The SUPERD_T type is
 25:  * used to facilitate argument transmission to the ptr2f function. 25:  * used to facilitate argument transmission to the ptr2f function.
 26:  */ 26:  */
 27:  27: 
 28: /* sf344 addition */ 28: /* sf344 addition */
 29: /*#include "nabinterface.c"*/ 29: /*#include "nabinterface.c"*/
 30: /* end sf344 addition */ 30: /* end sf344 addition */
 31:  31: 
 32: /* csw34> Adding math.h to allow float comparisons and fabs */ 
 33: #include math.h 
 34:  
 35: #ifdef SCALAPACK 32: #ifdef SCALAPACK
 36:  33: 
 37: typedef struct superd { 34: typedef struct superd {
 38:   INT_T dtype, ctxt, m, n, mb, nb, rsrc, csrc, lld, row, col; 35:   INT_T dtype, ctxt, m, n, mb, nb, rsrc, csrc, lld, row, col;
 39:   size_t size; 36:   size_t size;
 40:   INT_T *divmb, *divnb, *modmb, *modnb, *divrw, *divcl, *modrw, *modcl; 37:   INT_T *divmb, *divnb, *modmb, *modnb, *divrw, *divcl, *modrw, *modcl;
 41: } SUPERD_T; 38: } SUPERD_T;
 42:  39: 
 43: #endif 40: #endif
 44:  41: 
2825:           /* "standard" (HCT) effective radii:  */2822:           /* "standard" (HCT) effective radii:  */
2826: 2823: 
2827:           reff[i] = 1.0/(ri1i + sumi);2824:           reff[i] = 1.0/(ri1i + sumi);
2828: 2825: 
2829:         } else {2826:         } else {
2830: 2827: 
2831:           /* "gbao" formulas:  */2828:           /* "gbao" formulas:  */
2832: 2829: 
2833:           psi[i] = -ri*sumi;2830:           psi[i] = -ri*sumi;
2834: 2831: 
2835: /* csw34> Catch for people trying to use igb=5 when asking for second derivatives. 
2836:  *        There is an inconsistancy between AMBER and NAB - but as we're planning 
2837:  *        on moving to interface AMBER 12 where this is fixed, for now we just 
2838:  *        want to stop people falling into the trap! 
2839:  */ 
2840:           if ( fabs(gbgamma - 4.850) < 1e-3 ) { 
2841:              fprintf(stderr,"WARNING: You are requesting second derivatives for igb=5!\n"); 
2842:              fprintf(stderr,"         The igb=5 energy and gradient is inconsistant between AMBER and NAB\n"); 
2843:              fprintf(stderr,"         Exiting for your protection!\n"); 
2844:              exit(1); 
2845:           } 
2846:  
2847:           reff[i] =1.0/( ri1i - tanh( (GBALPHA + GBGAMMA*psi[i]*psi[i])2832:           reff[i] =1.0/( ri1i - tanh( (GBALPHA + GBGAMMA*psi[i]*psi[i])
2848:                                  *psi[i])/rborn[i]);2833:                                  *psi[i])/rborn[i]);
2849:         }2834:         }
2850: #endif2835: #endif
2851: 2836: 
2852:       }2837:       }
2853:     }2838:     }
2854: 2839: 
2855:     if (gb2_debug) {2840:     if (gb2_debug) {
2856:       printf("reff calculated, mytaskid = %d\n", mytaskid);2841:       printf("reff calculated, mytaskid = %d\n", mytaskid);


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0