[1] | 1 | /*************************************************************************** |
---|
| 2 | module classe ta_c.h - description |
---|
| 3 | ***************************************************************************/ |
---|
| 4 | // mohamed.berrada@upmc.fr |
---|
| 5 | // locean-ipsl.upmc, Paris, March 18, 2010 |
---|
| 6 | //=========================================================================== |
---|
| 7 | // methode forward |
---|
[59] | 8 | forward (YREAL x1) |
---|
[58] | 9 | { |
---|
[67] | 10 | YS1 = ta_eb(Yi,Yj,Yk) +norm_fil*x1; |
---|
[62] | 11 | if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1 || Yk==NZ-1){ |
---|
| 12 | YS1=0.; |
---|
| 13 | } |
---|
| 14 | |
---|
| 15 | |
---|
[58] | 16 | ///* 1..NPCA from pca_ta 1..NPCA 1*/ |
---|
| 17 | //double sum=moy_ta(Yi,Yj,Yk); |
---|
| 18 | //for(int n=0;n<NPCA;n++) sum+=Yting[n]*stdev_ta[n]*shfs_ta(Yi,Yj,Yk,n); |
---|
| 19 | //YS1=sum; |
---|
| 20 | //// |
---|
[67] | 21 | // printf("norm_fac= %24.16e\n",norm_fac); |
---|
[68] | 22 | if(fabs(norm_fil*x1)>5) |
---|
[67] | 23 | printf("increm(%d,%d,%d) = %24.16e (%24.16e) \n",Yi+1,Yj+1,Yk+1,norm_fil*x1,YS1); |
---|
[1] | 24 | } |
---|
| 25 | //=========================================================================== |
---|
| 26 | // methode backward |
---|
| 27 | |
---|
[59] | 28 | backward (YREAL x1) |
---|
[1] | 29 | { |
---|
[67] | 30 | YJ1I1=norm_fil; |
---|
[62] | 31 | if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1 || Yk==NZ-1){ |
---|
| 32 | YJ1I1=0.; |
---|
| 33 | } |
---|
| 34 | |
---|
[58] | 35 | //for(int n=0;n<NPCA;n++) |
---|
| 36 | // Yjac[0][n] = stdev_ta[n]*shfs_ta(Yi,Yj,Yk,n); |
---|
| 37 | //// |
---|
[1] | 38 | } |
---|
| 39 | |
---|
| 40 | //=========================================================================== |
---|
| 41 | //*************************** FIN DU MODULE ta_c *************************** |
---|
| 42 | |
---|