Changeset 86 for altifloat/src/Yworkdir/Y2floater_delta.h
- Timestamp:
- 04/01/13 19:31:21 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
altifloat/src/Yworkdir/Y2floater_delta.h
r84 r86 2816 2816 Ytps = Ytttt_pdt(0, YTemps, YidTraj); 2817 2817 Yting[15]=YS8_xinter( Yi, Ytps); 2818 Yting[16]=YS1_r_float_d( Yi, YTemps-1); 2819 Yting[17]=YS2_r_float_d( Yi, YTemps-1); 2818 Ytps = Ytttt_pdt(0, YTemps-1, YidTraj); 2819 Yting[16]=YS1_r_float( Yi, Ytps); 2820 Ytps = Ytttt_pdt(0, YTemps-1, YidTraj); 2821 Yting[17]=YS2_r_float( Yi, Ytps); 2820 2822 Yur_d(Yi)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] ); 2821 2823 Yting[0]=YS1_r_float_d( Yi, YTemps-1); … … 2863 2865 Ytps = Ytttt_pdt(0, YTemps, YidTraj); 2864 2866 Yting[15]=YS8_xinter( Yi, Ytps); 2865 Yting[16]=YS1_r_float_d( Yi, YTemps-1); 2866 Yting[17]=YS2_r_float_d( Yi, YTemps-1); 2867 Ytps = Ytttt_pdt(0, YTemps-1, YidTraj); 2868 Yting[16]=YS1_r_float( Yi, Ytps); 2869 Ytps = Ytttt_pdt(0, YTemps-1, YidTraj); 2870 Yting[17]=YS2_r_float( Yi, Ytps); 2867 2871 Ytbeta[0]=YG1_uinter_d( Yi, YTemps); 2868 2872 Ytbeta[1]=YG2_uinter_d( Yi, YTemps); … … 2889 2893 Ytps = Ytttt_pdt(0, YTemps, YidTraj); 2890 2894 Ytbeta[15]=YG8_xinter( Yi, Ytps); 2891 Ytbeta[16]=YG1_r_float_d( Yi, YTemps-1); 2892 Ytbeta[17]=YG2_r_float_d( Yi, YTemps-1); 2895 Ytps = Ytttt_pdt(0, YTemps-1, YidTraj); 2896 Ytbeta[16]=YG1_r_float( Yi, Ytps); 2897 Ytps = Ytttt_pdt(0, YTemps-1, YidTraj); 2898 Ytbeta[17]=YG2_r_float( Yi, Ytps); 2893 2899 memset(Yjac, 0, 36*sizeof(double)); 2894 2900 Yur_d(Yi)->backward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] ); … … 2958 2964 Ytps = Ytttt_pdt(0, YTemps, YidTraj); 2959 2965 Yting[15]=YS8_xinter( Yi, Ytps); 2960 Yting[16]=YS1_r_float_d( Yi, YTemps-1); 2961 Yting[17]=YS2_r_float_d( Yi, YTemps-1); 2966 Ytps = Ytttt_pdt(0, YTemps-1, YidTraj); 2967 Yting[16]=YS1_r_float( Yi, Ytps); 2968 Ytps = Ytttt_pdt(0, YTemps-1, YidTraj); 2969 Yting[17]=YS2_r_float( Yi, Ytps); 2962 2970 if (!strcmp(nmmod, "ur_d") || All) 2963 2971 { … … 3023 3031 Ytps = Ytttt_pdt(0, YTemps, YidTraj); 3024 3032 Yting[15]=YS8_xinter( Yi, Ytps); 3025 Yting[16]=YS1_r_float_d( Yi, YTemps-1); 3026 Yting[17]=YS2_r_float_d( Yi, YTemps-1); 3033 Ytps = Ytttt_pdt(0, YTemps-1, YidTraj); 3034 Yting[16]=YS1_r_float( Yi, Ytps); 3035 Ytps = Ytttt_pdt(0, YTemps-1, YidTraj); 3036 Yting[17]=YS2_r_float( Yi, Ytps); 3027 3037 if (!strcmp(nmmod, "ur_d") || All) 3028 3038 { … … 3113 3123 Ytps = Ytttt_pdt(0, YTemps, YidTraj); 3114 3124 Yting[15]=YS8_xinter( Yi, Ytps); 3115 Yting[16]=YS1_r_float_d( Yi, YTemps-1); 3116 Yting[17]=YS2_r_float_d( Yi, YTemps-1); 3125 Ytps = Ytttt_pdt(0, YTemps-1, YidTraj); 3126 Yting[16]=YS1_r_float( Yi, Ytps); 3127 Ytps = Ytttt_pdt(0, YTemps-1, YidTraj); 3128 Yting[17]=YS2_r_float( Yi, Ytps); 3117 3129 memset(Yjac, 0, 36*sizeof(double)); 3118 3130 Yur_d(Yi)->backward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] ); … … 3142 3154 Ytps = Ytttt_pdt(0, YTemps, YidTraj); 3143 3155 YG8_xinter( Yi, Ytps) += Ytbeta[15]; 3144 YG1_r_float_d( Yi, YTemps-1) += Ytbeta[16]; 3145 YG2_r_float_d( Yi, YTemps-1) += Ytbeta[17]; 3156 Ytps = Ytttt_pdt(0, YTemps-1, YidTraj); 3157 YG1_r_float( Yi, Ytps) += Ytbeta[16]; 3158 Ytps = Ytttt_pdt(0, YTemps-1, YidTraj); 3159 YG2_r_float( Yi, Ytps) += Ytbeta[17]; 3146 3160 if (Ytestad_module) { 3147 3161 Yting[0]=YD1_uinter_d( Yi, YTemps); … … 3169 3183 Ytps = Ytttt_pdt(0, YTemps, YidTraj); 3170 3184 Yting[15]=YD8_xinter( Yi, Ytps); 3171 Yting[16]=YD1_r_float_d( Yi, YTemps-1); 3172 Yting[17]=YD2_r_float_d( Yi, YTemps-1); 3185 Ytps = Ytttt_pdt(0, YTemps-1, YidTraj); 3186 Yting[16]=YD1_r_float( Yi, Ytps); 3187 Ytps = Ytttt_pdt(0, YTemps-1, YidTraj); 3188 Yting[17]=YD2_r_float( Yi, Ytps); 3173 3189 YLTRes = Yprosca(&YD1Y_ur_d, &YG1Y_ur_d, YNBS_ur_d); 3174 3190 YAdRes = Yprosca(Yting, Ytbeta, YNBI_ur_d);
Note: See TracChangeset
for help on using the changeset viewer.