source: altifloat/src/locate_d.h @ 84

forward(YREAL x, YREAL y)
{
        //Yj = longitude
        //Yk = latitude
        int j_ll=YS1_xinter(Yi,YTemps);
        int k_ll=YS5_xinter(Yi,YTemps);
        if (j_ll==Yj & k_ll==Yk)
        {
                //printf("j_ll=%d ; k_ll=%d\n",j_ll,k_ll);
                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
                YS2_uinter_d(Yi,YTemps)=YS1_u_d(Yj+1,Yk);
                YS3_uinter_d(Yi,YTemps)=YS1_u_d(Yj,Yk+1);
                YS4_uinter_d(Yi,YTemps)=YS1_u_d(Yj+1,Yk+1);
                YS5_uinter_d(Yi,YTemps)=YS1_v_d(Yj,Yk);
                YS6_uinter_d(Yi,YTemps)=YS1_v_d(Yj+1,Yk);
                YS7_uinter_d(Yi,YTemps)=YS1_v_d(Yj,Yk+1);
                YS8_uinter_d(Yi,YTemps)=YS1_v_d(Yj+1,Yk+1);
                
                
        }
}

backward(YREAL x, YREAL y)
{
        int j_ll=YS1_xinter(Yi,YTemps);
        int k_ll=YS5_xinter(Yi,YTemps);
        if (j_ll==Yj & k_ll==Yk)
        {
                YG1_u_d(Yj,Yk)+=YG1_uinter_d(Yi,YTemps);
                YG1_u_d(Yj+1,Yk)+=YG2_uinter_d(Yi,YTemps);
                YG1_u_d(Yj,Yk+1)+=YG3_uinter_d(Yi,YTemps);
                YG1_u_d(Yj+1,Yk+1)+=YG4_uinter_d(Yi,YTemps);
                YG1_v_d(Yj,Yk)+=YG5_uinter_d(Yi,YTemps);
                YG1_v_d(Yj+1,Yk)+=YG6_uinter_d(Yi,YTemps);
                YG1_v_d(Yj,Yk+1)+=YG7_uinter_d(Yi,YTemps);
                YG1_v_d(Yj+1,Yk+1)+=YG8_uinter_d(Yi,YTemps);

        }


//YG1_uinter(Yi,Ytemps)=


YG1_u(Yj,Yk)


        
        
}