Source
1067
1067
{
1068
1068
1069
1069
1070
1070
Matrix<Float> work = ( vecWork_p[scale] )(blc,trc);
1071
1071
1072
1072
work = 0.0;
1073
1073
for(Int taylor1=0;taylor1<ntaylor;taylor1++)
1074
1074
{
1075
1075
Matrix<Float> coeffs1 = (matCoeffs_p[IND2(taylor1,scale)])(blc,trc);
1076
1076
Matrix<Float> resid = (matR_p[IND2(taylor1,scale)])(blc,trc);
1077
+
1078
+
/*
1077
1079
work = work + (Float)2.0 * coeffs1 * resid;
1078
1080
for(Int taylor2=0;taylor2<ntaylor;taylor2++)
1079
1081
{
1080
1082
Matrix<Float> coeffs2 = (matCoeffs_p[IND2(taylor2,scale)])(blc,trc);
1081
1083
work = work - (Float)((matA_p[scale])(taylor1,taylor2)) * coeffs1 * coeffs2;
1082
1084
}
1085
+
*/
1086
+
1087
+
work = work + coeffs1 * resid;
1088
+
1083
1089
}
1084
1090
findMaxAbsMask(vecWork_p[scale],vecScaleMasks_p[scale],maxScaleVal_p[scale],maxScalePos_p[scale]);
1085
1091
1086
1092
/*
1087
1093
1088
1094
vecWork_p[scale] = 0.0;
1089
1095
for(Int taylor1=0;taylor1<ntaylor;taylor1++)
1090
1096
{
1091
1097
vecWork_p[scale] = vecWork_p[scale] + (Float)2.0 * ( (matCoeffs_p[IND2(taylor1,scale)]) * (matR_p[IND2(taylor1,scale)]) );
1092
1098
for(Int taylor2=0;taylor2<ntaylor;taylor2++)