source: altifloat/src/locate_d.h @ 108

Last change on this file since 108 was 108, checked in by leila_ocean, 11 years ago

with filter
File size: 1.2 KB
Line 
1forward(YREAL x, YREAL y)
2{
3        //Yj = longitude
4        //Yk = latitude
5        int j_ll=YS1_xinter(Yi,YTemps);
6        int k_ll=YS5_xinter(Yi,YTemps);
7        if (j_ll==Yj && k_ll==Yk)
8        {
9                //printf("j_ll=%d ; k_ll=%d\n",j_ll,k_ll);
10                YS1_uinter_d(Yi,YTemps)=YS1_u_d(Yj,Yk);//la 1ere sortie du module uinter pour le floater i est la valeur du module u en Yj Yk
11                YS2_uinter_d(Yi,YTemps)=YS1_u_d(Yj+1,Yk);
12                YS3_uinter_d(Yi,YTemps)=YS1_u_d(Yj,Yk+1);
13                YS4_uinter_d(Yi,YTemps)=YS1_u_d(Yj+1,Yk+1);
14                YS5_uinter_d(Yi,YTemps)=YS1_v_d(Yj,Yk);
15                YS6_uinter_d(Yi,YTemps)=YS1_v_d(Yj+1,Yk);
16                YS7_uinter_d(Yi,YTemps)=YS1_v_d(Yj,Yk+1);
17                YS8_uinter_d(Yi,YTemps)=YS1_v_d(Yj+1,Yk+1);
18               
19                //printf("u_d=%f\n",YS1_u_d(Yj,Yk));
20
21               
22               
23        }
24}
25
26backward(YREAL x, YREAL y)
27{
28        int j_ll=YS1_xinter(Yi,YTemps);
29        int k_ll=YS5_xinter(Yi,YTemps);
30        if (j_ll==Yj && k_ll==Yk)
31        {
32                YG1_u_d(Yj,Yk)+=YG1_uinter_d(Yi,YTemps);
33                YG1_u_d(Yj+1,Yk)+=YG2_uinter_d(Yi,YTemps);
34                YG1_u_d(Yj,Yk+1)+=YG3_uinter_d(Yi,YTemps);
35                YG1_u_d(Yj+1,Yk+1)+=YG4_uinter_d(Yi,YTemps);
36                YG1_v_d(Yj,Yk)+=YG5_uinter_d(Yi,YTemps);
37                YG1_v_d(Yj+1,Yk)+=YG6_uinter_d(Yi,YTemps);
38                YG1_v_d(Yj,Yk+1)+=YG7_uinter_d(Yi,YTemps);
39                YG1_v_d(Yj+1,Yk+1)+=YG8_uinter_d(Yi,YTemps);
40
41        }
42
43
44//YG1_uinter(Yi,Ytemps)=


YG1_u(Yj,Yk)


45       
46       
47}
Note: See TracBrowser for help on using the repository browser.