Rev | Line | |
---|
[84] | 1 | forward(YREAL xd, YREAL yd, YREAL ud, YREAL vd) |
---|
| 2 | { |
---|
[87] | 3 | |
---|
| 4 | //what we had up to 01/04 |
---|
| 5 | /*YS1=xd+rdtflo*ud+xd*rdtflo*(ur_dx+ur_dy);//mot cle YS1: yao sortie 1 |
---|
[85] | 6 | printf("xd=%f 2term=%f 3term=%f\n",xd,rdtflo*ud,xd*rdtflo*(ur_dx+ur_dy)); |
---|
[87] | 7 | YS2=yd+rdtflo*vd+yd*rdtflo*(vr_dx+vr_dy);*/ |
---|
| 8 | |
---|
| 9 | //what we should have corrected on 04/04 |
---|
| 10 | YS1=xd+rdtflo*ud+rdtflo*(xd*ur_dx+yd*ur_dy); |
---|
| 11 | YS2=yd+rdtflo*vd+rdtflo*(xd*vr_dx+yd*vr_dy); |
---|
[108] | 12 | printf("ud=%f\n",ud); |
---|
[84] | 13 | } |
---|
| 14 | |
---|
| 15 | backward(YREAL xd, YREAL yd, YREAL ud, YREAL vd) |
---|
| 16 | { |
---|
[87] | 17 | /*YJ1I1=1+rdtflo*(ur_dx+ur_dy); //mot cle, les entrees du gradient du out par rapport au in, dYs1/dx,y,u,v |
---|
[84] | 18 | YJ1I2=0; |
---|
| 19 | YJ1I3=rdtflo; |
---|
| 20 | YJ1I4=0; |
---|
| 21 | |
---|
| 22 | YJ2I1=0; //dYs2/dx,y,u,v |
---|
| 23 | YJ2I2=1+rdtflo*(vr_dx+vr_dy); |
---|
| 24 | YJ2I3=0; |
---|
[87] | 25 | YJ2I4=rdtflo;*/ |
---|
| 26 | |
---|
| 27 | |
---|
| 28 | //what we sould have corrected on 04/04 |
---|
| 29 | YJ1I1=1+rdtflo*ur_dx; |
---|
| 30 | YJ1I2=rdtflo*ur_dy; |
---|
| 31 | YJ1I3=rdtflo; |
---|
| 32 | YJ1I4=0; |
---|
| 33 | |
---|
| 34 | YJ2I1=rdtflo*vr_dx; |
---|
| 35 | YJ2I2=1+rdtflo*vr_dy; |
---|
| 36 | YJ2I3=0; |
---|
[84] | 37 | YJ2I4=rdtflo; |
---|
| 38 | |
---|
| 39 | } |
---|
Note: See
TracBrowser
for help on using the repository browser.