/[ascend]/trunk/models/johnpye/fprops/helmholtz.c
ViewVC logotype

Diff of /trunk/models/johnpye/fprops/helmholtz.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 2298 by jpye, Sun Aug 15 10:11:54 2010 UTC revision 2299 by jpye, Fri Aug 20 06:53:38 2010 UTC
# Line 1182  double helm_resid_tautau(double tau, dou Line 1182  double helm_resid_tautau(double tau, dou
1182  */  */
1183  double helm_resid_deldeldel(double tau,double delta,const HelmholtzData *data){  double helm_resid_deldeldel(double tau,double delta,const HelmholtzData *data){
1184      double sum = 0, res = 0;      double sum = 0, res = 0;
1185      double D,L;      double D;
1186      unsigned n, i;      unsigned n, i;
1187      const HelmholtzPowTerm *pt;      const HelmholtzPowTerm *pt;
1188      const HelmholtzGausTerm *gt;      const HelmholtzGausTerm *gt;
# Line 1200  double helm_resid_deldeldel(double tau,d Line 1200  double helm_resid_deldeldel(double tau,d
1200    
1201      /* seem to be errors in the following, still haven't tracked them all down. */      /* seem to be errors in the following, still haven't tracked them all down. */
1202    
1203  #if 0  #if 1
1204      /* wxmaxima code:      /* wxmaxima code:
1205          a*delta^d*%e^(-delta^l)*tau^t          a*delta^d*%e^(-delta^l)*tau^t
1206          diff(%,delta,3);          diff(%,delta,3);
# Line 1235  double helm_resid_deldeldel(double tau,d Line 1235  double helm_resid_deldeldel(double tau,d
1235      //fprintf(stderr,"DELDELDEL fiff = %f, sum = %f  ",fdiff, res);      //fprintf(stderr,"DELDELDEL fiff = %f, sum = %f  ",fdiff, res);
1236  #endif  #endif
1237    
1238  #if 0  #if 1
1239      /* gaussian terms */      /* gaussian terms */
1240      n = data->ng;      n = data->ng;
1241      //fprintf(stderr,"THERE ARE %d GAUSSIAN TERMS\n",n);      //fprintf(stderr,"THERE ARE %d GAUSSIAN TERMS\n",n);
# Line 1264  double helm_resid_deldeldel(double tau,d Line 1264  double helm_resid_deldeldel(double tau,d
1264      }      }
1265  #endif  #endif
1266    
1267  #if 0  #if 1
1268      /* critical terms */      /* critical terms */
1269      n = data->nc;      n = data->nc;
1270      ct = &(data->ct[0]);      ct = &(data->ct[0]);

Legend:
Removed from v.2298  
changed lines
  Added in v.2299

john.pye@anu.edu.au
ViewVC Help
Powered by ViewVC 1.1.22