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

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

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

revision 2110 by jpye, Mon Dec 7 04:13:05 2009 UTC revision 2111 by jpye, Tue Dec 8 03:05:36 2009 UTC
# Line 18  Line 18 
18  */  */
19    
20  #include "carbondioxide.h"  #include "carbondioxide.h"
21    
22  #define CARBONDIOXIDE_R 188.9241  #define GAS_C 8314.510
23  #define CARBONDIOXIDE_TSTAR 304.1282  #define CARBONDIOXIDE_TSTAR 304.1282
24  #define CARBONDIOXIDE_M 44.0098  #define CARBONDIOXIDE_M 44.0098
25    #define CARBONDIOXIDE_R (GAS_C/CARBONDIOXIDE_M)
26    
27  #define CARBONDIOXIDE_PC 73773e6  #define CARBONDIOXIDE_PC 73773e6
28  #define CARBONDIOXIDE_RHOC 467.6  #define CARBONDIOXIDE_RHOC 467.6
29    
# Line 32  Line 34 
34      Ideal gas data for CO2      Ideal gas data for CO2
35  */  */
36  const IdealData ideal_data_carbondioxide = {  const IdealData ideal_data_carbondioxide = {
37      1.5759375242e+03/CARBONDIOXIDE_R /* constant, adjust to solve s */      -1.1571354956e+003/CARBONDIOXIDE_R /* constant, adjust to solve s */
38      , (-2.1189540883e+05 + 2.6222371192e+05)/CARBONDIOXIDE_TSTAR/CARBONDIOXIDE_R /* linear, adjust to solver h */      , 2.9392750129e+005/CARBONDIOXIDE_TSTAR/CARBONDIOXIDE_R /* linear, adjust to solver h */
39      , CARBONDIOXIDE_TSTAR /* Tstar / [K] */      , CARBONDIOXIDE_TSTAR /* Tstar / [K] */
40      , CARBONDIOXIDE_R /* cpstar / [J/kgK] */      , CARBONDIOXIDE_R /* cpstar / [J/kgK] */
41      , 1 /* power terms */      , 1 /* power terms */
# Line 117  const HelmholtzData helmholtz_data_carbo Line 119  const HelmholtzData helmholtz_data_carbo
119  };  };
120    
121  /*  /*
122      Test suite. These tests attempt to validate the current code using      Test suite. These tests attempt to validate the current code using
123      a few sample figures output by REFPROP 7.0.      a few sample figures output by REFPROP 7.0.
124    
125      To run the test, compile and run as follows:      To run the test, compile and run as follows:
126    
127      gcc ideal.c helmholtz.c carbondioxide.c -DTEST -o carbondioxide -lm && ./carbondioxide      gcc ideal.c helmholtz.c carbondioxide.c -DTEST -o carbondioxide -lm && ./carbondioxide
128    
129      These tests all currently pass with a maximum error of 0.09%.      These tests all currently pass with a maximum error of 0.09%.
130  */  */
131  #ifdef TEST  #ifdef TEST
132    #include "test.h"
133    
134  #include <assert.h>  #include <assert.h>
135  #include <stdlib.h>  #include <stdlib.h>
136  #include <stdio.h>  #include <stdio.h>
137  #include <math.h>  #include <math.h>
138    
 typedef struct{double T,p,rho,h,s,cp,cp0,w;} TestData;  
139  const TestData td[]; const unsigned ntd;  const TestData td[]; const unsigned ntd;
140    
141  int main(void){  int main(void){
# Line 145  int main(void){ Line 147  int main(void){
147      double maxerr = 0;      double maxerr = 0;
148    
149      n = ntd;      n = ntd;
     fprintf(stderr,"Running through %d test points...\n",n);  
   
 /* a simple macro to actually do the testing */  
 #define ASSERT_TOL(FN,PARAM1,PARAM2,PARAM3,VAL,TOL) {\  
         double cval; cval = FN(PARAM1,PARAM2,PARAM3);\  
         double err; err = cval - (double)(VAL);\  
         double relerrpc = (cval-(VAL))/(VAL)*100;\  
         if(fabs(relerrpc)>maxerr)maxerr=fabs(relerrpc);\  
         if(fabs(err)>fabs(TOL)){\  
             fprintf(stderr,"ERROR in line %d: value of '%s(%f,%f,%s)' = %0.8f,"\  
                 " should be %f, error is %.10e (%.2f%%)!\n"\  
                 , __LINE__, #FN,PARAM1,PARAM2,#PARAM3, cval, VAL,cval-(VAL)\  
                 ,relerrpc\  
             );\  
             exit(1);\  
         }else{\  
             fprintf(stderr,"    OK, %s(%f,%f,%s) = %8.2e with %.6f%% err.\n"\  
                 ,#FN,PARAM1,PARAM2,#PARAM3,VAL,relerrpc\  
             );\  
         }\  
     }  
   
 #define CP0(T,RHO,DATA) helmholtz_cp0(T,DATA)  
     fprintf(stderr,"CP0 TESTS\n");  
     for(i=0; i<n;++i){  
         cp0 = td[i].cp0*1e3;  
         //fprintf(stderr,"T = %f --> cp0 = %f\n", td[i].T, cp0);  
         ASSERT_TOL(CP0, td[i].T, 0., d, cp0, cp0*1e-1);  
     }  
 #undef CP0  
150    
151      fprintf(stderr,"PRESSURE TESTS\n");  #if USING_CITED_REFERENCE_STATES
     for(i=0; i<n;++i){  
         p = td[i].p*1e6;  
         //fprintf(stderr,"T = %f, rho = %f --> h = %f\n", td[i].T, td[i].rho, p);  
         ASSERT_TOL(helmholtz_p, td[i].T, td[i].rho, d, p, p*2e-4);  
     }  
   
 #if 0  
152      fprintf(stderr,"REFERENCE POINT CHECK\n");      fprintf(stderr,"REFERENCE POINT CHECK\n");
153      /* solve rho to give p = PREF */      /* solve rho to give p = PREF */
154      T = CARBONDIOXIDE_TREF;      T = CARBONDIOXIDE_TREF;
# Line 204  int main(void){ Line 169  int main(void){
169      }      }
170      exit(1);      exit(1);
171  #endif  #endif
172    
173    
174        return helm_run_test_cases(d, ntd, td, 'K');
175    
176      fprintf(stderr,"ENTHALPY TESTS\n");      //fprintf(stderr,"Tests completed OK (maximum error = %0.2f%% (%5.2e))\n",maxerr,maxerr);
177      for(i=0; i<n;++i){      //exit(0);
         T = td[i].T+273.15;  
         rho = td[i].rho;  
         h = td[i].h*1e3;  
         //fprintf(stderr,"%.20e\n",(h - helmholtz_h(T,rho,d)) );  
         ASSERT_TOL(helmholtz_h, td[i].T+273.15, td[i].rho, d, h, 1E3);  
     }  
   
     /* entropy offset required to attain agreement with REFPROP */  
     fprintf(stderr,"ENTROPY TESTS\n");  
     for(i=0; i<n;++i){  
         T = td[i].T+273.15;  
         rho = td[i].rho;  
         s = td[i].s*1e3;  
         //double err = s - helmholtz_s(T,rho,d);  
         //se += err; ss += err*err;  
         ASSERT_TOL(helmholtz_s, td[i].T+273.15, td[i].rho, d, s, 1e-1*s);  
         //fprintf(stderr,"%.20e\t%.20e\t%.20e\n",T,rho,err);  
     }  
   
     fprintf(stderr,"Tests completed OK (maximum error = %0.2f%% (%5.2e))\n",maxerr,maxerr);  
     exit(0);  
178  }  }
179    
180  /*  /*
181      Some test data generated by REFPROP 7.0 for p=0.1, 1, 10, 20, 100 MPa.      Some test data generated by REFPROP 7.0 for p=0.1, 1, 10, 20, 100 MPa.
182  */  */
183  const TestData td[] = {  const TestData td[] = {
184  /*, {Temperature, Pressure, Density, Enthalpy, Entropy, Cp, Cp0, Sound Speed}*/  /*, {Temperature, Pressure, Density, Int. Energy, Enthalpy, Entropy, Cv, Cp, Cp0, Helmholtz}
185  /*, {(K), (MPa), (kg/m³), (kJ/kg), (kJ/kg-K), (kJ/kg-K), (kJ/kg-K), (m/s)}*/  , {(K), (MPa), (kg/m), (kJ/kg), (kJ/kg), (kJ/kg-K), (kJ/kg-K), (kJ/kg-K), (kJ/kg-K), (kJ/kg)}*/
186   {2.16E+2, 1.E-2, 2.454062144E-1, 4.409287843E+2, 2.91905617E+0, 7.555917727E-1, 7.53186339E-1, 2.331957752E+2}   {2.16E+2, 1.E-2, 2.454062144E-1, 4.0018002E+2, 4.409287843E+2, 2.91905617E+0, 5.653701779E-1, 7.555917727E-1, 7.53186339E-1, -2.303361127E+2}
187  , {2.66E+2, 1.E-2, 1.991337854E-1, 4.80078336E+2, 3.081890365E+0, 8.104414236E-1, 8.094088439E-1, 2.559309622E+2}  , {2.66E+2, 1.E-2, 1.991337854E-1, 4.298608404E+2, 4.80078336E+2, 3.081890365E+0, 6.208933259E-1, 8.104414236E-1, 8.094088439E-1, -3.899219966E+2}
188  , {3.16E+2, 1.E-2, 1.675727826E-1, 5.219266469E+2, 3.225924974E+0, 8.62745442E-1, 8.621948769E-1, 2.764359366E+2}  , {3.16E+2, 1.E-2, 1.675727826E-1, 4.622510848E+2, 5.219266469E+2, 3.225924974E+0, 6.734607926E-1, 8.62745442E-1, 8.621948769E-1, -5.571412072E+2}
189  , {3.66E+2, 1.E-2, 1.446570587E-1, 5.662693649E+2, 3.356110837E+0, 9.10105514E-1, 9.097699996E-1, 2.953724734E+2}  , {3.66E+2, 1.E-2, 1.446570587E-1, 4.971403497E+2, 5.662693649E+2, 3.356110837E+0, 7.209482393E-1, 9.10105514E-1, 9.097699996E-1, -7.311962166E+2}
190  , {4.16E+2, 1.E-2, 1.272587039E-1, 6.128544267E+2, 3.475359028E+0, 9.525151636E-1, 9.522915251E-1, 3.130934511E+2}  , {4.16E+2, 1.E-2, 1.272587039E-1, 5.342743398E+2, 6.128544267E+2, 3.475359028E+0, 7.634284574E-1, 9.525151636E-1, 9.522915251E-1, -9.114750158E+2}
191  , {4.66E+2, 1.E-2, 1.135979244E-1, 6.614481847E+2, 3.58562638E+0, 9.905630878E-1, 9.904041983E-1, 3.298200692E+2}  , {4.66E+2, 1.E-2, 1.135979244E-1, 5.734184073E+2, 6.614481847E+2, 3.58562638E+0, 8.015193523E-1, 9.905630878E-1, 9.904041983E-1, -1.097483486E+3}
192  , {5.16E+2, 1.E-2, 1.02586621E-1, 7.118486226E+2, 3.688334069E+0, 1.024880263E+0, 1.024761888E+0, 3.457079762E+2}  , {5.16E+2, 1.E-2, 1.02586621E-1, 6.143700244E+2, 7.118486226E+2, 3.688334069E+0, 8.358645706E-1, 1.024880263E+0, 1.024761888E+0, -1.288810355E+3}
193  , {5.66E+2, 1.E-2, 9.352189158E-2, 7.638828536E+2, 3.784562017E+0, 1.055991905E+0, 1.055900437E+0, 3.608749549E+2}  , {5.66E+2, 1.E-2, 9.352189158E-2, 6.569560172E+2, 7.638828536E+2, 3.784562017E+0, 8.669955064E-1, 1.055991905E+0, 1.055900437E+0, -1.485106085E+3}
194  , {6.16E+2, 1.E-2, 8.592935538E-2, 8.17401209E+2, 3.875154775E+0, 1.084304186E+0, 1.084231447E+0, 3.754140217E+2}  , {6.16E+2, 1.E-2, 8.592935538E-2, 7.010265433E+2, 8.17401209E+2, 3.875154775E+0, 8.953216241E-1, 1.084304186E+0, 1.084231447E+0, -1.686068798E+3}
195  , {6.66E+2, 1.E-2, 7.947723266E-2, 8.722718632E+2, 3.960786755E+0, 1.110129832E+0, 1.110070631E+0, 3.894004519E+2}  , {6.66E+2, 1.E-2, 7.947723266E-2, 7.464496665E+2, 8.722718632E+2, 3.960786755E+0, 9.211575282E-1, 1.110129832E+0, 1.110070631E+0, -1.891434312E+3}
196  , {7.16E+2, 1.E-2, 7.392649286E-2, 9.283769181E+2, 4.04200591E+0, 1.133717699E+0, 1.133668593E+0, 4.028960679E+2}  , {7.16E+2, 1.E-2, 7.392649286E-2, 7.931074144E+2, 9.283769181E+2, 4.04200591E+0, 9.447532117E-1, 1.133717699E+0, 1.133668593E+0, -2.100968817E+3}
197  , {7.66E+2, 1.E-2, 6.910055494E-2, 9.856098142E+2, 4.119264426E+0, 1.155275584E+0, 1.155234203E+0, 4.159521691E+2}  , {7.66E+2, 1.E-2, 6.910055494E-2, 8.408931762E+2, 9.856098142E+2, 4.119264426E+0, 9.663171895E-1, 1.155275584E+0, 1.155234203E+0, -2.314463374E+3}
198  , {8.16E+2, 1.E-2, 6.486613978E-2, 1.043873678E+3, 4.192940888E+0, 1.174984966E+0, 1.174949628E+0, 4.286116931E+2}  , {8.16E+2, 1.E-2, 6.486613978E-2, 8.897100411E+2, 1.043873678E+3, 4.192940888E+0, 9.860314115E-1, 1.174984966E+0, 1.174949628E+0, -2.531729723E+3}
199  , {8.66E+2, 1.E-2, 6.112076176E-2, 1.103080237E+3, 4.263356512E+0, 1.193009454E+0, 1.192978932E+0, 4.40910879E+2}  , {8.66E+2, 1.E-2, 6.112076176E-2, 9.394697103E+2, 1.103080237E+3, 4.263356512E+0, 1.00405981E+0, 1.193009454E+0, 1.192978932E+0, -2.752597029E+3}
200  , {9.16E+2, 1.E-2, 5.778432004E-2, 1.16314905E+3, 4.330787182E+0, 1.209499163E+0, 1.209472541E+0, 4.528805691E+2}  , {9.16E+2, 1.E-2, 5.778432004E-2, 9.900917223E+2, 1.16314905E+3, 4.330787182E+0, 1.020552721E+0, 1.209499163E+0, 1.209472541E+0, -2.976909336E+3}
201  , {9.66E+2, 1.E-2, 5.479330284E-2, 1.22400689E+3, 4.395472433E+0, 1.22459275E+0, 1.224569332E+0, 4.645472351E+2}  , {9.66E+2, 1.E-2, 5.479330284E-2, 1.041502834E+3, 1.22400689E+3, 4.395472433E+0, 1.035648963E+0, 1.22459275E+0, 1.224569332E+0, -3.204523536E+3}
202  , {1.016E+3, 1.E-2, 5.209670393E-2, 1.285587197E+3, 4.457622249E+0, 1.238418302E+0, 1.238397548E+0, 4.759337887E+2}  , {1.016E+3, 1.E-2, 5.209670393E-2, 1.093636473E+3, 1.285587197E+3, 4.457622249E+0, 1.049476742E+0, 1.238418302E+0, 1.238397548E+0, -3.435307731E+3}
203  , {1.066E+3, 1.E-2, 4.965308888E-2, 1.347829557E+3, 4.517422242E+0, 1.251093748E+0, 1.251075233E+0, 4.87060222E+2}  , {1.066E+3, 1.E-2, 4.965308888E-2, 1.146432218E+3, 1.347829557E+3, 4.517422242E+0, 1.062154072E+0, 1.251093748E+0, 1.251075233E+0, -3.669139892E+3}
204  , {2.16E+2, 1.E-1, 2.486826186E+0, 4.391151258E+2, 2.478112397E+0, 7.781063746E-1, 7.53186339E-1, 2.314231315E+2}  , {2.16E+2, 1.E-1, 2.486826186E+0, 3.989032282E+2, 4.391151258E+2, 2.478112397E+0, 5.755764177E-1, 7.781063746E-1, 7.53186339E-1, -1.363690496E+2}
205  , {2.66E+2, 1.E-1, 2.004455136E+0, 4.789986937E+2, 2.644049351E+0, 8.199018006E-1, 8.094088439E-1, 2.549699686E+2}  , {2.66E+2, 1.E-1, 2.004455136E+0, 4.291098246E+2, 4.789986937E+2, 2.644049351E+0, 6.246026263E-1, 8.199018006E-1, 8.094088439E-1, -2.742073029E+2}
206  , {3.16E+2, 1.E-1, 1.68194314E+0, 5.211919289E+2, 2.789283073E+0, 8.677495114E-1, 8.621948769E-1, 2.758590947E+2}  , {3.16E+2, 1.E-1, 1.68194314E+0, 4.61736887E+2, 5.211919289E+2, 2.789283073E+0, 6.751779588E-1, 8.677495114E-1, 8.621948769E-1, -4.196765641E+2}
207  , {3.66E+2, 1.E-1, 1.449818149E+0, 5.657298133E+2, 2.920045807E+0, 9.131429027E-1, 9.097699996E-1, 2.950131586E+2}  , {3.66E+2, 1.E-1, 1.449818149E+0, 4.967556457E+2, 5.657298133E+2, 2.920045807E+0, 7.218710051E-1, 9.131429027E-1, 9.097699996E-1, -5.719811197E+2}
208  , {4.16E+2, 1.E-1, 1.274385234E+0, 6.124388048E+2, 3.03961277E+0, 9.545352243E-1, 9.522915251E-1, 3.128704572E+2}  , {4.16E+2, 1.E-1, 1.274385234E+0, 5.339695968E+2, 6.124388048E+2, 3.03961277E+0, 7.639780203E-1, 9.545352243E-1, 9.522915251E-1, -7.305093154E+2}
209  , {4.66E+2, 1.E-1, 1.137005248E+0, 6.611176673E+2, 3.150073938E+0, 9.919963349E-1, 9.904041983E-1, 3.296882772E+2}  , {4.66E+2, 1.E-1, 1.137005248E+0, 5.731673256E+2, 6.611176673E+2, 3.150073938E+0, 8.018724634E-1, 9.919963349E-1, 9.904041983E-1, -8.947671292E+2}
210  , {5.16E+2, 1.E-1, 1.026454403E+0, 7.115799429E+2, 3.252907994E+0, 1.025947113E+0, 1.024761888E+0, 3.456400516E+2}  , {5.16E+2, 1.E-1, 1.026454403E+0, 6.141572033E+2, 7.115799429E+2, 3.252907994E+0, 8.361052864E-1, 1.025947113E+0, 1.024761888E+0, -1.064343322E+3}
211  , {5.66E+2, 1.E-1, 9.355473622E-1, 7.636610538E+2, 3.349222831E+0, 1.056815769E+0, 1.055900437E+0, 3.608532909E+2}  , {5.66E+2, 1.E-1, 9.355473622E-1, 6.567717566E+2, 7.636610538E+2, 3.349222831E+0, 8.671676509E-1, 1.056815769E+0, 1.055900437E+0, -1.238888366E+3}
212  , {6.16E+2, 1.E-1, 8.594627307E-1, 8.17216138E+2, 3.439877873E+0, 1.084959095E+0, 1.084231447E+0, 3.754267419E+2}  , {6.16E+2, 1.E-1, 8.594627307E-1, 7.008643796E+2, 8.17216138E+2, 3.439877873E+0, 8.954497237E-1, 1.084959095E+0, 1.084231447E+0, -1.41810039E+3}
213  , {6.66E+2, 1.E-1, 7.948415686E-1, 8.721163302E+2, 3.525556019E+0, 1.110662697E+0, 1.110070631E+0, 3.894392505E+2}  , {6.66E+2, 1.E-1, 7.948415686E-1, 7.463050944E+2, 8.721163302E+2, 3.525556019E+0, 9.212561086E-1, 1.110662697E+0, 1.110070631E+0, -1.601715214E+3}
214  , {7.16E+2, 1.E-1, 7.392705665E-1, 9.282456486E+2, 3.606810343E+0, 1.134159609E+0, 1.133668593E+0, 4.029549639E+2}  , {7.16E+2, 1.E-1, 7.392705665E-1, 7.929771765E+2, 9.282456486E+2, 3.606810343E+0, 9.448312846E-1, 1.134159609E+0, 1.133668593E+0, -1.789499029E+3}
215  , {7.66E+2, 1.E-1, 6.909704835E-1, 9.854988267E+2, 3.684096267E+0, 1.155647918E+0, 1.155234203E+0, 4.160267505E+2}  , {7.66E+2, 1.E-1, 6.909704835E-1, 8.407748444E+2, 9.854988267E+2, 3.684096267E+0, 9.663805679E-1, 1.155647918E+0, 1.155234203E+0, -1.981242896E+3}
216  , {8.16E+2, 1.E-1, 6.486004016E-1, 1.043779893E+3, 3.757794502E+0, 1.175302894E+0, 1.174949628E+0, 4.286986388E+2}  , {8.16E+2, 1.E-1, 6.486004016E-1, 8.896017587E+2, 1.043779893E+3, 3.757794502E+0, 9.860839726E-1, 1.175302894E+0, 1.174949628E+0, -2.176758555E+3}
217  , {8.66E+2, 1.E-1, 6.111303977E-1, 1.103001225E+3, 3.828227711E+0, 1.193284034E+0, 1.192978932E+0, 4.410076452E+2}  , {8.66E+2, 1.E-1, 6.111303977E-1, 9.39370026E+2, 1.103001225E+3, 3.828227711E+0, 1.004104215E+0, 1.193284034E+0, 1.192978932E+0, -2.375875172E+3}
218  , {9.16E+2, 1.E-1, 5.777562235E-1, 1.16308286E+3, 3.895672784E+0, 1.209738641E+0, 1.209472541E+0, 4.529851787E+2}  , {9.16E+2, 1.E-1, 5.777562235E-1, 9.899994803E+2, 1.16308286E+3, 3.895672784E+0, 1.020590846E+0, 1.209738641E+0, 1.209472541E+0, -2.57843679E+3}
219  , {9.66E+2, 1.E-1, 5.478406447E-1, 1.223951931E+3, 3.960369979E+0, 1.2248034E+0, 1.224569332E+0, 4.646581315E+2}  , {9.66E+2, 1.E-1, 5.478406447E-1, 1.041417099E+3, 1.223951931E+3, 3.960369979E+0, 1.035682158E+0, 1.2248034E+0, 1.224569332E+0, -2.784300301E+3}
220  , {1.016E+3, 1.E-1, 5.208721919E-1, 1.285542153E+3, 4.022529808E+0, 1.238604985E+0, 1.238397548E+0, 4.760497327E+2}  , {1.016E+3, 1.E-1, 5.208721919E-1, 1.093556477E+3, 1.285542153E+3, 4.022529808E+0, 1.049506E+0, 1.238604985E+0, 1.238397548E+0, -2.993333808E+3}
221  , {1.066E+3, 1.E-1, 4.964355705E-1, 1.34779333E+3, 4.082338276E+0, 1.251260285E+0, 1.251075233E+0, 4.871802172E+2}  , {1.066E+3, 1.E-1, 4.964355705E-1, 1.146357321E+3, 1.34779333E+3, 4.082338276E+0, 1.062180136E+0, 1.251260285E+0, 1.251075233E+0, -3.205415281E+3}
222  , {2.20E+2, 1.E+0, 1.16703286E+3, 8.683239354E+1, 5.505734566E-1, 1.95885155E+0, 7.576942616E-1, 9.535503477E+2}  , {2.66E+2, 5.1795E-1, 1.072031066E+1, 4.254876367E+2, 4.738024679E+2, 2.319596193E+0, 6.427453842E-1, 8.683836586E-1, 8.094088439E-1, -1.915249506E+2}
223  , {2.330282499E+2, 1.E+0, 1.116903622E+3, 1.126571928E+2, 6.646002577E-1, 2.011052229E+0, 7.724408378E-1, 8.571842464E+2}  , {3.16E+2, 5.1795E-1, 8.866368952E+0, 4.59302791E+2, 5.17720166E+2, 2.470828925E+0, 6.833341373E-1, 8.922113512E-1, 8.621948769E-1, -3.214791494E+2}
224  , {2.330282499E+2, 1.E+0, 2.60056432E+1, 4.352962968E+2, 2.049149574E+0, 1.032158413E+0, 7.724408378E-1, 2.235011365E+2}  , {3.66E+2, 5.1795E-1, 7.588812072E+0, 4.949499667E+2, 5.632017556E+2, 2.60438178E+0, 7.262056716E-1, 9.276786804E-1, 9.097699996E-1, -4.582537647E+2}
225  , {2.70E+2, 1.E+0, 2.114700655E+1, 4.711107444E+2, 2.191977498E+0, 9.323476551E-1, 8.138035995E-1, 2.46913928E+2}  , {4.16E+2, 5.1795E-1, 6.6442523E+0, 5.325457223E+2, 6.105003178E+2, 2.725465421E+0, 7.665452876E-1, 9.640907818E-1, 9.522915251E-1, -6.012478927E+2}
226  , {3.20E+2, 1.E+0, 1.723430305E+1, 5.17280852E+2, 2.348878807E+0, 9.243272791E-1, 8.661970555E-1, 2.718047464E+2}  , {4.66E+2, 5.1795E-1, 5.9138178E+0, 5.719973032E+2, 6.595803194E+2, 2.836840196E+0, 8.035169363E-1, 9.987288679E-1, 9.904041983E-1, -7.499702279E+2}
227  , {3.70E+2, 1.E+0, 1.465943954E+1, 5.640248879E+2, 2.484565304E+0, 9.476337552E-1, 9.133502547E-1, 2.930513634E+2}  , {5.16E+2, 5.1795E-1, 5.330607466E+0, 6.131671242E+2, 7.103324104E+2, 2.940266523E+0, 8.372243184E-1, 1.030936206E+0, 1.024761888E+0, -9.040104016E+2}
228  , {4.20E+2, 1.E+0, 1.279586846E+1, 6.121542004E+2, 2.606534037E+0, 9.779757625E-1, 9.554935853E-1, 3.121629304E+2}  , {5.66E+2, 5.1795E-1, 4.853493743E+0, 6.559154804E+2, 7.626324203E+2, 3.036987084E+0, 8.679670625E-1, 1.060657083E+0, 1.055900437E+0, -1.063019209E+3}
229  , {4.70E+2, 1.E+0, 1.137110751E+1, 6.618337726E+2, 2.718258346E+0, 1.009113763E+0, 9.932830471E-1, 3.297928738E+2}  , {6.16E+2, 5.1795E-1, 4.455590046E+0, 7.001113458E+2, 8.163585756E+2, 3.127932232E+0, 8.96044226E-1, 1.088006386E+0, 1.084231447E+0, -1.226694909E+3}
230  , {5.20E+2, 1.E+0, 1.024105219E+1, 7.130457217E+2, 2.821779384E+0, 1.039098422E+0, 1.027365251E+0, 3.462974606E+2}  , {6.66E+2, 5.1795E-1, 4.11849327E+0, 7.456340944E+2, 8.713961064E+2, 3.213825034E+0, 9.2171345E-1, 1.113138541E+0, 1.110070631E+0, -1.394773378E+3}
231  , {5.70E+2, 1.E+0, 9.32034508E+0, 7.657137375E+2, 2.918465999E+0, 1.067303834E+0, 1.058265531E+0, 3.61900332E+2}  , {7.16E+2, 5.1795E-1, 3.829145224E+0, 7.923729248E+2, 9.276381001E+2, 3.295242675E+0, 9.451934164E-1, 1.136210722E+0, 1.133668593E+0, -1.567020831E+3}
232  , {6.20E+2, 1.E+0, 8.554490597E+0, 8.197435523E+2, 3.009310433E+0, 1.093561114E+0, 1.08638662E+0, 3.767553743E+2}  , {7.66E+2, 5.1795E-1, 3.578005866E+0, 8.402259753E+2, 9.849853802E+2, 3.372655764E+0, 9.66674516E-1, 1.1573748E+0, 1.155234203E+0, -1.74322834E+3}
233  , {6.70E+2, 1.E+0, 7.906773232E+0, 8.750373567E+2, 3.095068043E+0, 1.117870908E+0, 1.112037962E+0, 3.909751084E+2}  , {8.16E+2, 5.1795E-1, 3.357934289E+0, 8.890995969E+2, 1.043346213E+3, 3.446454959E+0, 9.863277451E-1, 1.176776641E+0, 1.174949628E+0, -1.923207649E+3}
234  , {7.20E+2, 1.E+0, 7.351424646E+0, 9.31499309E+2, 3.176333646E+0, 1.140301565E+0, 1.135466079E+0, 4.046452569E+2}  , {8.66E+2, 5.1795E-1, 3.163474379E+0, 9.38907801E+2, 1.102636012E+3, 3.516969662E+0, 1.004310168E+0, 1.194556335E+0, 1.192978932E+0, -2.106787927E+3}
235  , {7.70E+2, 1.E+0, 6.86976529E+0, 9.890378182E+2, 3.253587797E+0, 1.160950947E+0, 1.156877252E+0, 4.178330492E+2}  , {9.16E+2, 5.1795E-1, 2.990382931E+0, 9.895718066E+2, 1.162777049E+3, 3.584481464E+0, 1.020767676E+0, 1.210847972E+0, 1.209472541E+0, -2.293813215E+3}
236  , {8.20E+2, 1.E+0, 6.447892424E+0, 1.047566572E+3, 3.327226892E+0, 1.179930587E+0, 1.176451949E+0, 4.305923792E+2}  , {9.66E+2, 5.1795E-1, 2.835308799E+0, 1.041019622E+3, 1.223698136E+3, 3.649233983E+0, 1.03583613E+0, 1.225778994E+0, 1.224569332E+0, -2.484140406E+3}
237  , {8.70E+2, 1.E+0, 6.075228187E+0, 1.107005004E+3, 3.397583678E+0, 1.197357862E+0, 1.194353025E+0, 4.429672231E+2}  , {1.016E+3, 5.1795E-1, 2.695569229E+0, 1.093185622E+3, 1.285334277E+3, 3.71144018E+0, 1.049641719E+0, 1.239469451E+0, 1.238397548E+0, -2.6776376E+3}
238  , {9.20E+2, 1.E+0, 5.743572329E+0, 1.167278461E+3, 3.464941685E+0, 1.213351251E+0, 1.210729923E+0, 4.549940129E+2}  , {1.066E+3, 5.1795E-1, 2.568990553E+0, 1.146010121E+3, 1.347626277E+3, 3.771287889E+0, 1.06230105E+0, 1.252031389E+0, 1.251075233E+0, -2.874182768E+3}
239  , {9.70E+2, 1.E+0, 5.446464232E+0, 1.228318172E+3, 3.529545622E+0, 1.22802709E+0, 1.225720645E+0, 4.667033381E+2}  , {2.330282499E+2, 1.E+0, 1.116903622E+3, 1.117618604E+2, 1.126571928E+2, 6.646002577E-1, 9.53030319E-1, 2.011052229E+0, 7.724408378E-1, -4.310877453E+1}
240  , {1.020E+3, 1.E+0, 5.178741503E+0, 1.290061079E+3, 3.591609022E+0, 1.24149727E+0, 1.239452608E+0, 4.781211965E+2}  , {2.330282499E+2, 1.E+0, 2.60056432E+1, 3.968431044E+2, 4.352962968E+2, 2.049149574E+0, 6.802616759E-1, 1.032158413E+0, 7.724408378E-1, -8.066663447E+1}
241  , {1.070E+3, 1.E+0, 4.936227415E+0, 1.352449573E+3, 3.651319963E+0, 1.253867674E+0, 1.252043051E+0, 4.892699275E+2}  , {2.66E+2, 1.E+0, 2.155818097E+1, 4.209870371E+2, 4.6737314E+2, 2.178030874E+0, 6.662310299E-1, 9.366092546E-1, 8.094088439E-1, -1.583691755E+2}
242  , {2.20E+2, 1.E+1, 1.185633186E+3, 8.944320189E+1, 5.276696788E-1, 1.903856175E+0, 7.576942616E-1, 1.003148747E+3}  , {3.16E+2, 1.E+0, 1.748525303E+1, 4.563947588E+2, 5.135858099E+2, 2.337259039E+0, 6.931434181E-1, 9.232229315E-1, 8.621948769E-1, -2.821790975E+2}
243  , {2.70E+2, 1.0E+1, 9.894607858E+2, 1.887248227E+2, 9.335353209E-1, 2.142451131E+0, 8.138035995E-1, 6.641525173E+2}  , {3.66E+2, 1.E+0, 1.483421042E+1, 4.928271476E+2, 5.602388919E+2, 2.474277192E+0, 7.313069027E-1, 9.453710516E-1, 9.097699996E-1, -4.127583047E+2}
244  , {3.20E+2, 1.0E+1, 4.482772002E+2, 3.629005817E+2, 1.513661952E+0, 7.617495479E+0, 8.661970555E-1, 2.191440164E+2}  , {4.16E+2, 1.E+0, 1.292629419E+1, 5.308856036E+2, 6.082472964E+2, 2.597187325E+0, 7.695363473E-1, 9.754769274E-1, 9.522915251E-1, -5.495443234E+2}
245  , {3.70E+2, 1.E+1, 1.927369325E+2, 4.99290267E+2, 1.917146569E+0, 1.557419115E+0, 9.133502547E-1, 2.659746253E+2}  , {4.66E+2, 1.E+0, 1.147281341E+1, 5.706396611E+2, 6.578022396E+2, 2.709643922E+0, 8.054224988E-1, 1.006651426E+0, 9.904041983E-1, -6.920544065E+2}
246  , {4.20E+2, 1.E+1, 1.474750648E+2, 5.679882127E+2, 2.091678028E+0, 1.261643459E+0, 9.554935853E-1, 2.987681682E+2}  , {5.16E+2, 1.E+0, 1.032286536E+1, 6.120216711E+2, 7.088939992E+2, 2.813764447E+0, 8.38516999E-1, 1.036760927E+0, 1.024761888E+0, -8.398807836E+2}
247  , {4.70E+2, 1.0E+1, 1.233191172E+2, 6.286736778E+2, 2.228277185E+0, 1.180883351E+0, 9.932830471E-1, 3.241813901E+2}  , {5.66E+2, 1.E+0, 9.387718047E+0, 6.549267294E+2, 7.61448889E+2, 2.910957441E+0, 8.688888527E-1, 1.065118347E+0, 1.055900437E+0, -9.92675182E+2}
248  , {5.20E+2, 1.0E+1, 1.073305338E+2, 6.869694278E+2, 2.346167334E+0, 1.15623326E+0, 1.027365251E+0, 3.457082453E+2}  , {6.16E+2, 1.E+0, 8.611007831E+0, 6.992429352E+2, 8.153733603E+2, 3.002238907E+0, 8.967290163E-1, 1.091532789E+0, 1.084231447E+0, -1.150136231E+3}
249  , {5.70E+2, 1.0E+1, 9.563339786E+1, 7.446444524E+2, 2.452069543E+0, 1.153027376E+0, 1.058265531E+0, 3.647789332E+2}  , {6.66E+2, 1.E+0, 7.954902838E+0, 7.448609802E+2, 8.70569618E+2, 3.088379797E+0, 9.222399371E-1, 1.11599643E+0, 1.110070631E+0, -1.311999965E+3}
250  , {6.20E+2, 1.0E+1, 8.656077872E+1, 8.024278418E+2, 2.549237927E+0, 1.159358586E+0, 1.08638662E+0, 3.82133223E+2}  , {7.16E+2, 1.E+0, 7.392928421E+0, 7.916771576E+2, 9.269415539E+2, 3.169985785E+0, 9.456101719E-1, 1.138574074E+0, 1.133668593E+0, -1.478032665E+3}
251  , {6.70E+2, 1.0E+1, 7.924771529E+1, 8.606533504E+2, 2.639549855E+0, 1.170151035E+0, 1.112037962E+0, 3.982078546E+2}  , {7.66E+2, 1.E+0, 6.905937903E+0, 8.39594265E+2, 9.843971888E+2, 3.2475453E+0, 9.670127599E-1, 1.159361973E+0, 1.155234203E+0, -1.648025435E+3}
252  , {7.20E+2, 1.0E+1, 7.318938174E+1, 9.194757706E+2, 2.724217432E+0, 1.182949222E+0, 1.135466079E+0, 4.132839338E+2}  , {8.16E+2, 1.E+0, 6.479708386E+0, 8.885218325E+2, 1.042849765E+3, 3.321460616E+0, 9.866082474E-1, 1.178470915E+0, 1.174949628E+0, -1.82179003E+3}
253  , {7.70E+2, 1.0E+1, 6.80662969E+1, 9.789600008E+2, 2.804086944E+0, 1.196469438E+0, 1.156877252E+0, 4.275537729E+2}  , {8.66E+2, 1.E+0, 6.103436029E+0, 9.383761121E+2, 1.102218249E+3, 3.392068976E+0, 1.004547163E+0, 1.19601801E+0, 1.192978932E+0, -1.999155621E+3}
254  , {8.20E+2, 1.0E+1, 6.366395382E+1, 1.039122905E+3, 2.879784045E+0, 1.210011066E+0, 1.176451949E+0, 4.411550682E+2}  , {9.16E+2, 1.E+0, 5.768757078E+0, 9.890799457E+2, 1.162427512E+3, 3.45965742E+0, 1.020971177E+0, 1.212121783E+0, 1.209472541E+0, -2.179966251E+3}
255  , {8.70E+2, 1.0E+1, 5.98316849E+1, 1.099954916E+3, 2.951791759E+0, 1.223186904E+0, 1.194353025E+0, 4.54189981E+2}  , {9.66E+2, 1.E+0, 5.469090071E+0, 1.040562545E+3, 1.223408318E+3, 3.524473452E+0, 1.036013344E+0, 1.226898838E+0, 1.224569332E+0, -2.36407881E+3}
256  , {9.20E+2, 1.0E+1, 5.645984003E+1, 1.161431955E+3, 3.020495782E+0, 1.235787709E+0, 1.210729923E+0, 4.667365265E+2}  , {1.016E+3, 1.E+0, 5.199181651E+0, 1.092759198E+3, 1.285097159E+3, 3.586732867E+0, 1.049797941E+0, 1.240461488E+0, 1.238397548E+0, -2.551361396E+3}
257  , {9.70E+2, 1.0E+1, 5.346630143E+1, 1.223522328E+3, 3.086212432E+0, 1.247709024E+0, 1.225720645E+0, 4.788557353E+2}  , {1.066E+3, 1.E+0, 4.954785492E+0, 1.145610919E+3, 1.347436004E+3, 3.646625605E+0, 1.062440245E+0, 1.252916138E+0, 1.251075233E+0, -2.741691976E+3}
258  , {1.020E+3, 1.0E+1, 5.078810138E+1, 1.286190821E+3, 3.149206742E+0, 1.258909342E+0, 1.239452608E+0, 4.905963451E+2}  , {2.66E+2, 1.E+1, 1.008230103E+3, 1.703194875E+2, 1.802378583E+2, 9.018677036E-1, 9.301612759E-1, 2.102126092E+0, 8.094088439E-1, -6.957732164E+1}
259  , {1.070E+3, 1.0E+1, 4.837599107E+1, 1.349401174E+3, 3.209704663E+0, 1.269385128E+0, 1.252043051E+0, 5.019979809E+2}  , {3.16E+2, 1.0E+1, 5.597379223E+2, 3.135772661E+2, 3.3144277E+2, 1.414740357E+0, 1.057816023E+0, 7.326747935E+0, 8.621948769E-1, -1.334806868E+2}
260    , {3.66E+2, 1.0E+1, 1.98484141E+2, 4.425772486E+2, 4.929591076E+2, 1.899941649E+0, 8.452988087E-1, 1.609567063E+0, 9.097699996E-1, -2.528013949E+2}
261    , {4.16E+2, 1.0E+1, 1.500068446E+2, 4.96255415E+2, 5.629190397E+2, 2.07955064E+0, 8.285930224E-1, 1.273130713E+0, 9.522915251E-1, -3.688376511E+2}
262    , {4.66E+2, 1.0E+1, 1.248716971E+2, 5.438608125E+2, 6.239430107E+2, 2.218168838E+0, 8.412988511E-1, 1.184505847E+0, 9.904041983E-1, -4.898058661E+2}
263    , {5.16E+2, 1.0E+1, 1.08420603E+2, 5.90109262E+2, 6.823426542E+2, 2.33723529E+0, 8.623121952E-1, 1.157175578E+0, 1.024761888E+0, -6.159041474E+2}
264    , {5.66E+2, 1.0E+1, 9.645739419E+1, 6.36360029E+2, 7.400327446E+2, 2.443950306E+0, 8.856716458E-1, 1.152836473E+0, 1.055900437E+0, -7.469158442E+2}
265    , {6.16E+2, 1.0E+1, 8.721333487E+1, 6.83130483E+2, 7.977918478E+2, 2.541736295E+0, 9.091353898E-1, 1.158643119E+0, 1.084231447E+0, -8.825790747E+2}
266    , {6.66E+2, 1.0E+1, 7.978137012E+1, 7.306321252E+2, 8.559746706E+2, 2.632545825E+0, 9.317628034E-1, 1.169190991E+0, 1.110070631E+0, -1.022643394E+3}
267    , {7.16E+2, 1.0E+1, 7.363621112E+1, 7.789433503E+2, 9.147461007E+2, 2.717630134E+0, 9.531498283E-1, 1.181886486E+0, 1.133668593E+0, -1.166879825E+3}
268    , {7.66E+2, 1.0E+1, 6.844726264E+1, 8.28078421E+2, 9.741763015E+2, 2.797858158E+0, 9.731393965E-1, 1.195380226E+0, 1.155234203E+0, -1.315080928E+3}
269    , {8.16E+2, 1.0E+1, 6.399347241E+1, 8.780190691E+2, 1.034285007E+3, 2.873869735E+0, 9.916976438E-1, 1.208937453E+0, 1.174949628E+0, -1.467058634E+3}
270    , {8.66E+2, 1.0E+1, 6.012007272E+1, 9.287304392E+2, 1.095064237E+3, 2.946157327E+0, 1.008855348E+0, 1.222152034E+0, 1.192978932E+0, -1.622641806E+3}
271    , {9.16E+2, 1.0E+1, 5.671471735E+1, 9.801696947E+2, 1.156490771E+3, 3.015113221E+0, 1.024677616E+0, 1.234803749E+0, 1.209472541E+0, -1.781674016E+3}
272    , {9.66E+2, 1.0E+1, 5.369344469E+1, 1.032290871E+3, 1.218533345E+3, 3.081058518E+0, 1.039246983E+0, 1.246781594E+0, 1.224569332E+0, -1.944011658E+3}
273    , {1.016E+3, 1.0E+1, 5.099198481E+1, 1.085047668E+3, 1.28115692E+3, 3.144261845E+0, 1.052653463E+0, 1.258040108E+0, 1.238397548E+0, -2.109522366E+3}
274    , {1.066E+3, 1.0E+1, 4.856014402E+1, 1.138395059E+3, 1.344325256E+3, 3.204951926E+0, 1.06498854E+0, 1.268573381E+0, 1.251075233E+0, -2.278083694E+3}
275    , {2.66E+2, 1.00E+2, 1.201223185E+3, 1.219547998E+2, 2.052032762E+2, 6.942957691E-1, 9.718044513E-1, 1.650019263E+0, 8.094088439E-1, -6.272787475E+1}
276    , {3.16E+2, 1.E+2, 1.096322907E+3, 1.952242518E+2, 2.864382539E+2, 9.74260521E-1, 9.371579106E-1, 1.602434682E+0, 8.621948769E-1, -1.126420729E+2}
277    , {3.66E+2, 1.00E+2, 9.964902957E+2, 2.652422026E+2, 3.655944092E+2, 1.206875731E+0, 9.241881172E-1, 1.564208585E+0, 9.097699996E-1, -1.764743149E+2}
278    , {4.16E+2, 1.00E+2, 9.042838571E+2, 3.322798721E+2, 4.428646174E+2, 1.404819045E+0, 9.23833771E-1, 1.526675777E+0, 9.522915251E-1, -2.521248506E+2}
279    , {4.66E+2, 1.00E+2, 8.216796311E+2, 3.965925437E+2, 5.182944774E+2, 1.576083662E+0, 9.308032174E-1, 1.491059241E+0, 9.904041983E-1, -3.378624427E+2}
280    , {5.16E+2, 1.00E+2, 7.493262617E+2, 4.585756675E+2, 5.920288842E+2, 1.726413415E+0, 9.418566817E-1, 1.458855999E+0, 1.024761888E+0, -4.322536548E+2}
281    , {5.66E+2, 1.00E+2, 6.869405554E+2, 5.186601111E+2, 6.642331148E+2, 1.85999357E+0, 9.55105217E-1, 1.429920847E+0, 1.055900437E+0, -5.340962496E+2}
282    , {6.16E+2, 1.00E+2, 6.336268853E+2, 5.772738229E+2, 7.350954089E+2, 1.979982167E+0, 9.694646332E-1, 1.405425192E+0, 1.084231447E+0, -6.423951922E+2}
283    , {6.66E+2, 1.00E+2, 5.881088279E+2, 6.348278144E+2, 8.048643711E+2, 2.088891096E+0, 9.842611237E-1, 1.38621029E+0, 1.110070631E+0, -7.563736553E+2}
284    , {7.16E+2, 1.00E+2, 5.490721929E+2, 6.916740689E+2, 8.737994821E+2, 2.188702252E+0, 9.990470337E-1, 1.371968091E+0, 1.133668593E+0, -8.754367438E+2}
285    , {7.66E+2, 1.00E+2, 5.153528008E+2, 7.480905022E+2, 9.421323311E+2, 2.280958121E+0, 1.013517439E+0, 1.361993577E+0, 1.155234203E+0, -9.991234187E+2}
286    , {8.16E+2, 1.00E+2, 4.85979865E+2, 8.042873443E+2, 1.010057186E+3, 2.366861016E+0, 1.027466929E+0, 1.355520557E+0, 1.174949628E+0, -1.127071245E+3}
287    , {8.66E+2, 1.00E+2, 4.601677832E+2, 8.604182512E+2, 1.077730292E+3, 2.447353251E+0, 1.040764859E+0, 1.351804563E+0, 1.192978932E+0, -1.258989664E+3}
288    , {9.16E+2, 1.00E+2, 4.372911706E+2, 9.165918456E+2, 1.145272429E+3, 2.523178446E+0, 1.053337501E+0, 1.350179618E+0, 1.209472541E+0, -1.394639611E+3}
289    , {9.66E+2, 1.00E+2, 4.168552987E+2, 9.728824079E+2, 1.212773805E+3, 2.594929012E+0, 1.065153173E+0, 1.350091707E+0, 1.224569332E+0, -1.533819018E+3}
290    , {1.016E+3, 1.00E+2, 3.984683178E+2, 1.029338886E+3, 1.280299867E+3, 2.663082573E+0, 1.076209938E+0, 1.351103331E+0, 1.238397548E+0, -1.676353008E+3}
291    , {1.066E+3, 1.00E+2, 3.81818007E+2, 1.085991969E+3, 1.34789685E+3, 2.72802938E+0, 1.086525894E+0, 1.352880503E+0, 1.251075233E+0, -1.822087351E+3}
292  };  };
293    
294  const unsigned ntd = sizeof(td)/sizeof(TestData);  const unsigned ntd = sizeof(td)/sizeof(TestData);

Legend:
Removed from v.2110  
changed lines
  Added in v.2111

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