| 686 |
|
|
| 687 |
|
char *foo; // the pointer to the current string token |
| 688 |
|
int j; |
| 689 |
+ |
double chi, integralOfChidt; |
| 690 |
+ |
double eta[9]; |
| 691 |
|
|
| 692 |
|
// set the string tokenizer |
| 693 |
|
|
| 787 |
|
} |
| 788 |
|
boxMat[8] = atof( foo ); |
| 789 |
|
|
| 790 |
+ |
return NULL; |
| 791 |
+ |
|
| 792 |
+ |
//get chi and integralOfChidt, they should appear by pair |
| 793 |
+ |
foo = strtok(NULL, " ,;\t\n"); |
| 794 |
+ |
if(foo != NULL){ |
| 795 |
+ |
chi = atof(foo); |
| 796 |
+ |
|
| 797 |
+ |
foo = strtok(NULL, " ,;\t\n"); |
| 798 |
+ |
if(foo == NULL){ |
| 799 |
+ |
sprintf( painCave.errMsg, |
| 800 |
+ |
"chi and integralOfChidt should appear by pair in %s\n", inName ); |
| 801 |
+ |
return strdup( painCave.errMsg ); |
| 802 |
+ |
} |
| 803 |
+ |
integralOfChidt = atof( foo ); |
| 804 |
+ |
|
| 805 |
+ |
//push chi and integralOfChidt into SimInfo::properties which can be |
| 806 |
+ |
//retrieved by integrator later |
| 807 |
+ |
DoubleData* chiValue = new DoubleData(); |
| 808 |
+ |
chiValue->setID(CHIVALUE_ID); |
| 809 |
+ |
chiValue->setData(chi); |
| 810 |
+ |
simnfo->addProperty(chiValue); |
| 811 |
+ |
|
| 812 |
+ |
DoubleData* integralOfChidtValue = new DoubleData(); |
| 813 |
+ |
integralOfChidtValue->setID(INTEGRALOFCHIDT_ID); |
| 814 |
+ |
integralOfChidtValue->setData(integralOfChidt); |
| 815 |
+ |
simnfo->addProperty(integralOfChidtValue); |
| 816 |
+ |
|
| 817 |
+ |
} |
| 818 |
+ |
else |
| 819 |
+ |
return NULL; |
| 820 |
+ |
|
| 821 |
+ |
//get eta |
| 822 |
+ |
for(int i = 0 ; i < 9; i++){ |
| 823 |
+ |
foo = strtok(NULL, " ,;\t"); |
| 824 |
+ |
if(foo == NULL){ |
| 825 |
+ |
sprintf( painCave.errMsg, |
| 826 |
+ |
"error in reading eta[%d] from %s\n", i, inName ); |
| 827 |
+ |
return strdup( painCave.errMsg ); |
| 828 |
+ |
} |
| 829 |
+ |
eta[i] = atof( foo ); |
| 830 |
+ |
} |
| 831 |
+ |
|
| 832 |
+ |
//push eta into SimInfo::properties which can be |
| 833 |
+ |
//retrieved by integrator later |
| 834 |
+ |
//simnfo->setBoxM( theBoxMat3 ); |
| 835 |
+ |
DoubleArrayData* etaValue = new DoubleArrayData(); |
| 836 |
+ |
etaValue->setID(ETAVALUE_ID); |
| 837 |
+ |
etaValue->setData(eta, 9); |
| 838 |
+ |
simnfo->addProperty(etaValue); |
| 839 |
+ |
|
| 840 |
+ |
|
| 841 |
|
return NULL; |
| 842 |
+ |
|
| 843 |
+ |
|
| 844 |
+ |
|
| 845 |
|
} |
| 846 |
|
|
| 847 |
|
|