source: altifloat/src/locate_d.h

forward(YREAL x, YREAL y)
{
        //Yj = longitude
        //Yk = latitude
        int j_ll=YS1_xinter_u(Yi,YTemps); //iIndexLowerLeft
        int k_ll=YS5_xinter_u(Yi,YTemps); //jIndexLowerLeft     
        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);             
                //printf("u_d=%f\n",YS1_u_d(Yj,Yk));    
        }
        
        j_ll=YS1_xinter_v(Yi,YTemps); //iIndexLowerLeft
        k_ll=YS5_xinter_v(Yi,YTemps); //jIndexLowerLeft 
        if (j_ll==Yj && k_ll==Yk)
        {
                //printf("j_ll=%d ; k_ll=%d\n",j_ll,k_ll);
                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);             
                //printf("u_d=%f\n",YS1_u_d(Yj,Yk));    
        }
}

backward(YREAL x, YREAL y)
{
        int j_ll=YS1_xinter_u(Yi,YTemps);
        int k_ll=YS5_xinter_u(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);
        }
        j_ll=YS1_xinter_v(Yi,YTemps);
        k_ll=YS5_xinter_v(Yi,YTemps);
        if (j_ll==Yj && k_ll==Yk)
        {
                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)


        
        
}