Changeset 2700
- Timestamp:
- 01/30/12 03:38:07 (4 months ago)
- Location:
- simo/trunk/src/models/Finland/c
- Files:
-
- 2 edited
-
DistributionmodelLibrary.c (modified) (30 diffs)
-
make_dll.py (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
simo/trunk/src/models/Finland/c/DistributionmodelLibrary.c
r2697 r2700 458 458 int i; 459 459 double low_boundary, high_boundary, low_mass, high_mass, stem_number, left_out; 460 //double dist_2[class_count][2]; 460 #ifdef __cplusplus 461 461 double** dist_2 = NULL; 462 462 dist_2 = new double*[class_count]; … … 467 467 } 468 468 } 469 #else 470 double dist_2[class_count][2]; 471 #endif 469 472 470 473 left_out = 0; … … 539 542 class_count_new = 3; 540 543 } 541 544 #ifdef __cplusplus 542 545 for (int i=0; i < class_count; i++) { 543 546 delete[] dist_2[i]; … … 546 549 delete[] dist_2; 547 550 dist_2 = NULL; 548 551 #endif 549 552 return class_count_new; 550 553 } … … 562 565 int i; 563 566 double low_boundary, high_boundary, low_mass, high_mass, stem_number, left_out; 564 //double dist_2[class_count][2]; 567 #ifdef __cplusplus 565 568 double** dist_2 = NULL; 566 569 dist_2 = new double*[class_count]; … … 571 574 } 572 575 } 576 #else 577 double dist_2[class_count][2]; 578 #endif 573 579 574 580 left_out = 0; … … 660 666 class_count_new = 3; 661 667 } 668 #ifdef __cplusplus 662 669 for (int i=0; i < class_count; i++) { 663 670 delete[] dist_2[i]; … … 666 673 delete[] dist_2; 667 674 dist_2 = NULL; 675 #endif 668 676 return class_count_new; 669 677 } … … 680 688 int i; 681 689 double low_boundary, high_boundary, low_mass, high_mass, stem_number, N_calc, left_out; 682 //double dist_2[class_count][2]; 690 691 #ifdef __cplusplus 683 692 double** dist_2 = NULL; 684 693 dist_2 = new double*[class_count]; … … 689 698 } 690 699 } 700 #else 701 double dist_2[class_count][2]; 702 #endif 691 703 692 704 N_calc = 0; … … 786 798 class_count_new = 3; 787 799 } 800 #ifdef __cplusplus 788 801 for (int i=0; i < class_count; i++) { 789 802 delete[] dist_2[i]; … … 792 805 delete[] dist_2; 793 806 dist_2 = NULL; 807 #endif 794 808 return class_count_new; 795 809 … … 835 849 int i; 836 850 double cur_dbh, cur_N, low_boundary, high_boundary, low_mass, high_mass, stem_number, left_out; 837 //double dist_2[class_count][2]; 851 852 #ifdef __cplusplus 838 853 double** dist_2 = NULL; 839 854 dist_2 = new double*[class_count]; … … 844 859 } 845 860 } 861 #else 862 double dist_2[class_count][2]; 863 #endif 846 864 847 865 left_out = 0; … … 922 940 class_count_new = 3; 923 941 } 942 #ifdef __cplusplus 924 943 for (int i=0; i < class_count; i++) { 925 944 delete[] dist_2[i]; … … 928 947 delete[] dist_2; 929 948 dist_2 = NULL; 949 #endif 930 950 return class_count_new; 931 951 } … … 3225 3245 double low_boundary, high_boundary, low_mass, high_mass, stem_number, left_out; 3226 3246 double sum_ppa, sum_ppa_d, ppa, ppa_d, D_g, step, counter1, counter2, counter3; 3227 //double dist_2[class_count][3]; 3228 3247 #ifdef __cplusplus 3229 3248 double** dist_2 = NULL; 3230 3249 dist_2 = new double*[class_count]; … … 3235 3254 } 3236 3255 } 3237 3256 #else 3257 double dist_2[class_count][3]; 3258 #endif 3238 3259 left_out = 0; 3239 3260 counter1 = 0; … … 3431 3452 } 3432 3453 } 3433 3454 #ifdef __cplusplus 3434 3455 for (int i=0; i < class_count; i++) { 3435 3456 delete[] dist_2[i]; … … 3438 3459 delete[] dist_2; 3439 3460 dist_2 = NULL; 3440 3461 #endif 3441 3462 return class_count_new; 3442 3463 … … 3555 3576 double low_boundary, high_boundary, low_mass, high_mass, stem_number, left_out; 3556 3577 double sum_ppa, sum_ppa_d, ppa, ppa_d, D_g, step; 3557 // double dist_2[class_count][3];3558 3578 // 3579 #ifdef __cplusplus 3559 3580 double** dist_2 = NULL; 3560 3581 dist_2 = new double*[class_count]; … … 3565 3586 } 3566 3587 } 3588 #else 3589 double dist_2[class_count][3]; 3590 #endif 3567 3591 3568 3592 left_out = 0; … … 3700 3724 } 3701 3725 } 3702 3726 3727 #ifdef __cplusplus 3703 3728 for (int i=0; i < class_count; i++) { 3704 3729 delete[] dist_2[i]; … … 3707 3732 delete[] dist_2; 3708 3733 dist_2 = NULL; 3734 #endif 3709 3735 3710 3736 return class_count_new; … … 3793 3819 double c = exp(-2.3205+0.0744*H_gM-0.0469*H_dom+0.2296*log(N)+1.2869/log(H_dom/H_gM+0.4)); 3794 3820 double h0 = b * pow( -log( 1.0 - 0.000000000000001 ), (1 / c)); 3795 double h100 = b * pow( -log( 1.0 - 0.999999999999999 ), (1 / c)); // HUOM! Weibull-jakauman prosenttipistejakauman yht älö?3821 double h100 = b * pow( -log( 1.0 - 0.999999999999999 ), (1 / c)); // HUOM! Weibull-jakauman prosenttipistejakauman yhtᅵlᅵ? 3796 3822 3797 3823 // how many height classes (??????) … … 3865 3891 double c = exp(-2.3205+0.0744*H_gM-0.0469*H_dom+0.2296*log(N)+1.2869/log(H_dom/H_gM+0.4)); 3866 3892 double h0 = b * pow( -log( 1.0 - 0.000000000000001 ), (1 / c)); 3867 double h100 = b * pow( -log( 1.0 - 0.999999999999999 ), (1 / c)); // HUOM! Weibull-jakauman prosenttipistejakauman yht älö?3893 double h100 = b * pow( -log( 1.0 - 0.999999999999999 ), (1 / c)); // HUOM! Weibull-jakauman prosenttipistejakauman yhtᅵlᅵ? 3868 3894 3869 3895 // how many height classes (??????) … … 4006 4032 h0 = b * pow( -log( 0 + 0.000000000000001 ), (1 / c)); 4007 4033 } 4008 double h100 = a + b * pow( -log( 1.0 - 0.999999999999999 ), (1 / c)); // HUOM! Weibull-jakauman prosenttipistejakauman yht älö?4034 double h100 = a + b * pow( -log( 1.0 - 0.999999999999999 ), (1 / c)); // HUOM! Weibull-jakauman prosenttipistejakauman yhtᅵlᅵ? 4009 4035 4010 4036 // how many height classes (??????) … … 4027 4053 return ret; 4028 4054 4029 //Skaalaus siten ett ämitattu H_ari toteutuu: kerroin=mitattu H_ari/jakaumasta laskettu H_ari. Siis jokaista pituutta4055 //Skaalaus siten ettᅵ mitattu H_ari toteutuu: kerroin=mitattu H_ari/jakaumasta laskettu H_ari. Siis jokaista pituutta 4030 4056 //korjataan kertoimella! 4031 4057 } … … 4058 4084 h0 = b * pow( -log( 0 + 0.000000000000001 ), (1 / c)); 4059 4085 } 4060 double h100 = a + b * pow( -log( 1.0 - 0.999999999999999 ), (1 / c)); // HUOM! Weibull-jakauman prosenttipistejakauman yht älö?4086 double h100 = a + b * pow( -log( 1.0 - 0.999999999999999 ), (1 / c)); // HUOM! Weibull-jakauman prosenttipistejakauman yhtᅵlᅵ? 4061 4087 4062 4088 // how many height classes (??????) … … 4573 4599 double c = exp(-2.4203+0.0895*H_gM-0.0637*H_dom+0.251*log(N)+1.2707/log(H_dom/H_gM+0.4)); 4574 4600 double h0 = b * pow( -log( 0 + 0.000000000000001 ), (1 / c)); 4575 double h100 = b * pow( -log( 1.0 - 0.999999999999999 ), (1 / c)); // HUOM! Weibull-jakauman prosenttipistejakauman yht älö?4601 double h100 = b * pow( -log( 1.0 - 0.999999999999999 ), (1 / c)); // HUOM! Weibull-jakauman prosenttipistejakauman yhtᅵlᅵ? 4576 4602 4577 4603 // how many height classes (??????) -
simo/trunk/src/models/Finland/c/make_dll.py
r2697 r2700 36 36 for filename in filenames[1:]: 37 37 print "Compiling %s..." % filename 38 os.system(' '.join(('gcc', fpic, '-c', dmode, posix_source, 39 '%s.h' % filename, '%s.c' % filename))) 38 to_run = ' '.join(('gcc', fpic, '-c', dmode, posix_source, 39 '%s.h' % filename, '%s.c' % filename)) 40 print to_run 41 os.system(to_run) 40 42 input=' '.join((input, '%s.o' % filename)) 41 43 print "output =", dllname
Note: See TracChangeset
for help on using the changeset viewer.
