42 #define SIMPLE_UOP(OPER) \
43 G4ErrorMatrixIter a=m.begin(); \
44 G4ErrorMatrixIter e=m.end(); \
45 for(;a!=e; a++) (*a) OPER t;
47 #define SIMPLE_BOP(OPER) \
48 G4ErrorMatrixIter a=m.begin(); \
49 G4ErrorMatrixConstIter b=mat2.m.begin(); \
50 G4ErrorMatrixIter e=m.end(); \
51 for(;a!=e; a++, b++) (*a) OPER (*b);
53 #define SIMPLE_TOP(OPER) \
54 G4ErrorMatrixConstIter a=mat1.m.begin(); \
55 G4ErrorMatrixConstIter b=mat2.m.begin(); \
56 G4ErrorMatrixIter t=mret.m.begin(); \
57 G4ErrorMatrixConstIter e=mat1.m.end(); \
58 for(;a!=e; a++, b++, t++) (*t) = (*a) OPER (*b);
62 #define CHK_DIM_2(r1,r2,c1,c2,fun) \
63 if (r1!=r2 || c1!=c2) { \
64 G4ErrorMatrix::error("Range error in Matrix function " #fun "(1)."); \
67 #define CHK_DIM_1(c1,r2,fun) \
69 G4ErrorMatrix::error("Range error in Matrix function " #fun "(2)."); \
75 :
m(p*q), nrow(p), ncol(q)
81 :
m(p*q), nrow(p), ncol(q)
98 for( ; a<b; a+=(ncol+1)) *a = 1.0;
100 error(
"Invalid dimension in G4ErrorMatrix(G4int,G4int,1).");
105 error(
"G4ErrorMatrix: initialization must be either 0 or 1.");
118 :
m(mat1.size), nrow(mat1.nrow), ncol(mat1.ncol), size(mat1.size)
124 :
m(mat1.nrow*mat1.nrow), nrow(mat1.nrow), ncol(mat1.nrow)
137 for(
G4int k=1;k<=j;k++)
140 if(j!=k) *mkj = *sjk;
160 {
error(
"G4ErrorMatrix::sub: Index out of range"); }
181 {
error(
"G4ErrorMatrix::sub: Index out of range"); }
221 for(;a<
e; a++, b++) (*b) = -(*a);
283 for (
G4int j=0; j<m1cols; j++)
293 (*pt) += temp * (*pb);
334 if (&mat1 ==
this) {
return *
this; }
358 if(os.flags() & std::ios::fixed)
359 { width = os.precision()+3; }
361 { width = os.precision()+7; }
367 os << q(irow,icol) <<
" ";
381 for (; pme < pl; pme++, pt+=
nrow)
383 if (pt >= ptl) { pt -= (
size-1); }
399 *(b++) = (*
f)(*(a++), ir, ic);
407 {
error(
"dfinv_matrix: G4ErrorMatrix is not NxN"); }
409 if (n==1) {
return 0; }
418 *m21 = -(*m22) * (*m11) * (*m21);
431 for (
G4int j=1;j<=im2;j++)
439 for (
G4int k=j;k<=im2;k++)
441 s31 += (*mkj) * (*(mik++));
442 s32 += (*(mjkp++)) * (*mkpi);
446 *mij = -(*mii) * (((*(mij-
n)))*( (*(mii-1)))+(s31));
452 *(mii-1) = -(*mii) * (*mimim) * (*(mii-1));
453 *(mimim+1) = -(*(mimim+1));
474 s33 += (*mikj) * (*(miik++));
484 for (
G4int k=j;k<=ni;k++)
486 s34 += *mikij * (*(miik++));
495 if (nxch==0)
return 0;
496 for (
G4int mq=1;mq<=nxch;mq++)
498 G4int k = nxch - mq + 1;
521 error(
"dfact_matrix: G4ErrorMatrix is not NxN");
538 G4int jrange = 0, jover = 1, junder = -1;
548 p = (std::fabs(*mjj));
551 for (
G4int i=j+1;i<=
n;i++)
553 q = (std::fabs(*(mij)));
582 ir[nxch] = (((j)<<12)+(k));
596 t = (std::fabs(det));
600 if (jfail == jrange) jfail = junder;
605 if (jfail==jrange) jfail = jover;
622 for (
G4int i=1;i<j;i++)
624 s11 += (*mik) * (*(mji++));
625 s12 += (*mijp) * (*(mki++));
630 *(mjk++) = -s11 * (*mjj);
631 *(mkjp) = -(((*(mjj+1)))*((*(mkjp-1)))+(s12));
639 if (nxch%2==1) det = -det;
640 if (jfail !=jrange) det = 0.0;
648 {
error(
"G4ErrorMatrix::invert: G4ErrorMatrix is not NxN"); }
653 if (
ncol > max_array)
657 ir =
new G4int [max_array+1];
665 G4double c11,c12,c13,c21,c22,c23,c31,c32,c33;
667 c11 = (*(
m.begin()+4)) * (*(
m.begin()+8))
668 - (*(
m.begin()+5)) * (*(
m.begin()+7));
669 c12 = (*(
m.begin()+5)) * (*(
m.begin()+6))
670 - (*(
m.begin()+3)) * (*(
m.begin()+8));
671 c13 = (*(
m.begin()+3)) * (*(
m.begin()+7))
672 - (*(
m.begin()+4)) * (*(
m.begin()+6));
673 c21 = (*(
m.begin()+7)) * (*(
m.begin()+2))
674 - (*(
m.begin()+8)) * (*(
m.begin()+1));
675 c22 = (*(
m.begin()+8)) * (*
m.begin())
676 - (*(
m.begin()+6)) * (*(
m.begin()+2));
677 c23 = (*(
m.begin()+6)) * (*(
m.begin()+1))
678 - (*(
m.begin()+7)) * (*
m.begin());
679 c31 = (*(
m.begin()+1)) * (*(
m.begin()+5))
680 - (*(
m.begin()+2)) * (*(
m.begin()+4));
681 c32 = (*(
m.begin()+2)) * (*(
m.begin()+3))
682 - (*
m.begin()) * (*(
m.begin()+5));
683 c33 = (*
m.begin()) * (*(
m.begin()+4))
684 - (*(
m.begin()+1)) * (*(
m.begin()+3));
685 t1 = std::fabs(*
m.begin());
686 t2 = std::fabs(*(
m.begin()+3));
687 t3 = std::fabs(*(
m.begin()+6));
692 temp = *(
m.begin()+6);
693 det = c23*c12-c22*c13;
698 det = c22*c33-c23*c32;
703 temp = *(
m.begin()+6);
704 det = c23*c12-c22*c13;
708 temp = *(
m.begin()+3);
709 det = c13*c32-c12*c33;
732 det = (*
m.begin())*(*(
m.begin()+3)) - (*(
m.begin()+1))*(*(
m.begin()+2));
739 temp = ss*(*(
m.begin()+3));
740 *(
m.begin()+1) *= -ss;
741 *(
m.begin()+2) *= -ss;
742 *(
m.begin()+3) = ss*(*
m.begin());
747 if ((*(
m.begin()))==0)
752 *(
m.begin()) = 1.0/(*(
m.begin()));
781 {
error(
"G4ErrorMatrix::determinant: G4ErrorMatrix is not NxN"); }
782 if (
ncol > max_array)
786 ir =
new G4int [max_array+1];
809 G4Exception(
"G4ErrorMatrix::error()",
"GEANT4e-Error",
974 - m[
F01]*Det3_123_023
975 + m[
F02]*Det3_123_013
976 - m[
F03]*Det3_123_012;
987 m[
F00] = Det3_123_123 * oneOverDet;
988 m[
F01] = Det3_023_123 * mn1OverDet;
989 m[
F02] = Det3_013_123 * oneOverDet;
990 m[
F03] = Det3_012_123 * mn1OverDet;
992 m[
F10] = Det3_123_023 * mn1OverDet;
993 m[
F11] = Det3_023_023 * oneOverDet;
994 m[
F12] = Det3_013_023 * mn1OverDet;
995 m[
F13] = Det3_012_023 * oneOverDet;
997 m[
F20] = Det3_123_013 * oneOverDet;
998 m[
F21] = Det3_023_013 * mn1OverDet;
999 m[
F22] = Det3_013_013 * oneOverDet;
1000 m[
F23] = Det3_012_013 * mn1OverDet;
1002 m[
F30] = Det3_123_012 * mn1OverDet;
1003 m[
F31] = Det3_023_012 * oneOverDet;
1004 m[
F32] = Det3_013_012 * mn1OverDet;
1005 m[
F33] = Det3_012_012 * oneOverDet;
1050 + m[
M12]*Det2_23_01;
1052 + m[
M13]*Det2_23_01;
1054 + m[
M14]*Det2_23_01;
1056 + m[
M13]*Det2_23_02;
1058 + m[
M14]*Det2_23_02;
1060 + m[
M14]*Det2_23_03;
1062 + m[
M13]*Det2_23_12;
1064 + m[
M14]*Det2_23_12;
1066 + m[
M14]*Det2_23_13;
1068 + m[
M14]*Det2_23_23;
1070 + m[
M12]*Det2_24_01;
1072 + m[
M13]*Det2_24_01;
1074 + m[
M14]*Det2_24_01;
1076 + m[
M13]*Det2_24_02;
1078 + m[
M14]*Det2_24_02;
1080 + m[
M14]*Det2_24_03;
1082 + m[
M13]*Det2_24_12;
1084 + m[
M14]*Det2_24_12;
1086 + m[
M14]*Det2_24_13;
1088 + m[
M14]*Det2_24_23;
1090 + m[
M12]*Det2_34_01;
1092 + m[
M13]*Det2_34_01;
1094 + m[
M14]*Det2_34_01;
1096 + m[
M13]*Det2_34_02;
1098 + m[
M14]*Det2_34_02;
1100 + m[
M14]*Det2_34_03;
1102 + m[
M13]*Det2_34_12;
1104 + m[
M14]*Det2_34_12;
1106 + m[
M14]*Det2_34_13;
1108 + m[
M14]*Det2_34_23;
1110 + m[
M22]*Det2_34_01;
1112 + m[
M23]*Det2_34_01;
1114 + m[
M24]*Det2_34_01;
1116 + m[
M23]*Det2_34_02;
1118 + m[
M24]*Det2_34_02;
1120 + m[
M24]*Det2_34_03;
1122 + m[
M23]*Det2_34_12;
1124 + m[
M24]*Det2_34_12;
1126 + m[
M24]*Det2_34_13;
1128 + m[
M24]*Det2_34_23;
1132 G4double Det4_0123_0123 = m[
M00]*Det3_123_123 - m[
M01]*Det3_123_023
1133 + m[
M02]*Det3_123_013 - m[
M03]*Det3_123_012;
1134 G4double Det4_0123_0124 = m[
M00]*Det3_123_124 - m[
M01]*Det3_123_024
1135 + m[
M02]*Det3_123_014 - m[
M04]*Det3_123_012;
1136 G4double Det4_0123_0134 = m[
M00]*Det3_123_134 - m[
M01]*Det3_123_034
1137 + m[
M03]*Det3_123_014 - m[
M04]*Det3_123_013;
1138 G4double Det4_0123_0234 = m[
M00]*Det3_123_234 - m[
M02]*Det3_123_034
1139 + m[
M03]*Det3_123_024 - m[
M04]*Det3_123_023;
1140 G4double Det4_0123_1234 = m[
M01]*Det3_123_234 - m[
M02]*Det3_123_134
1141 + m[
M03]*Det3_123_124 - m[
M04]*Det3_123_123;
1142 G4double Det4_0124_0123 = m[
M00]*Det3_124_123 - m[
M01]*Det3_124_023
1143 + m[
M02]*Det3_124_013 - m[
M03]*Det3_124_012;
1144 G4double Det4_0124_0124 = m[
M00]*Det3_124_124 - m[
M01]*Det3_124_024
1145 + m[
M02]*Det3_124_014 - m[
M04]*Det3_124_012;
1146 G4double Det4_0124_0134 = m[
M00]*Det3_124_134 - m[
M01]*Det3_124_034
1147 + m[
M03]*Det3_124_014 - m[
M04]*Det3_124_013;
1148 G4double Det4_0124_0234 = m[
M00]*Det3_124_234 - m[
M02]*Det3_124_034
1149 + m[
M03]*Det3_124_024 - m[
M04]*Det3_124_023;
1150 G4double Det4_0124_1234 = m[
M01]*Det3_124_234 - m[
M02]*Det3_124_134
1151 + m[
M03]*Det3_124_124 - m[
M04]*Det3_124_123;
1152 G4double Det4_0134_0123 = m[
M00]*Det3_134_123 - m[
M01]*Det3_134_023
1153 + m[
M02]*Det3_134_013 - m[
M03]*Det3_134_012;
1154 G4double Det4_0134_0124 = m[
M00]*Det3_134_124 - m[
M01]*Det3_134_024
1155 + m[
M02]*Det3_134_014 - m[
M04]*Det3_134_012;
1156 G4double Det4_0134_0134 = m[
M00]*Det3_134_134 - m[
M01]*Det3_134_034
1157 + m[
M03]*Det3_134_014 - m[
M04]*Det3_134_013;
1158 G4double Det4_0134_0234 = m[
M00]*Det3_134_234 - m[
M02]*Det3_134_034
1159 + m[
M03]*Det3_134_024 - m[
M04]*Det3_134_023;
1160 G4double Det4_0134_1234 = m[
M01]*Det3_134_234 - m[
M02]*Det3_134_134
1161 + m[
M03]*Det3_134_124 - m[
M04]*Det3_134_123;
1162 G4double Det4_0234_0123 = m[
M00]*Det3_234_123 - m[
M01]*Det3_234_023
1163 + m[
M02]*Det3_234_013 - m[
M03]*Det3_234_012;
1164 G4double Det4_0234_0124 = m[
M00]*Det3_234_124 - m[
M01]*Det3_234_024
1165 + m[
M02]*Det3_234_014 - m[
M04]*Det3_234_012;
1166 G4double Det4_0234_0134 = m[
M00]*Det3_234_134 - m[
M01]*Det3_234_034
1167 + m[
M03]*Det3_234_014 - m[
M04]*Det3_234_013;
1168 G4double Det4_0234_0234 = m[
M00]*Det3_234_234 - m[
M02]*Det3_234_034
1169 + m[
M03]*Det3_234_024 - m[
M04]*Det3_234_023;
1170 G4double Det4_0234_1234 = m[
M01]*Det3_234_234 - m[
M02]*Det3_234_134
1171 + m[
M03]*Det3_234_124 - m[
M04]*Det3_234_123;
1172 G4double Det4_1234_0123 = m[
M10]*Det3_234_123 - m[
M11]*Det3_234_023
1173 + m[
M12]*Det3_234_013 - m[
M13]*Det3_234_012;
1174 G4double Det4_1234_0124 = m[
M10]*Det3_234_124 - m[
M11]*Det3_234_024
1175 + m[
M12]*Det3_234_014 - m[
M14]*Det3_234_012;
1176 G4double Det4_1234_0134 = m[
M10]*Det3_234_134 - m[
M11]*Det3_234_034
1177 + m[
M13]*Det3_234_014 - m[
M14]*Det3_234_013;
1178 G4double Det4_1234_0234 = m[
M10]*Det3_234_234 - m[
M12]*Det3_234_034
1179 + m[
M13]*Det3_234_024 - m[
M14]*Det3_234_023;
1180 G4double Det4_1234_1234 = m[
M11]*Det3_234_234 - m[
M12]*Det3_234_134
1181 + m[
M13]*Det3_234_124 - m[
M14]*Det3_234_123;
1186 - m[
M01]*Det4_1234_0234
1187 + m[
M02]*Det4_1234_0134
1188 - m[
M03]*Det4_1234_0124
1189 + m[
M04]*Det4_1234_0123;
1198 G4double mn1OverDet = - oneOverDet;
1200 m[
M00] = Det4_1234_1234 * oneOverDet;
1201 m[
M01] = Det4_0234_1234 * mn1OverDet;
1202 m[
M02] = Det4_0134_1234 * oneOverDet;
1203 m[
M03] = Det4_0124_1234 * mn1OverDet;
1204 m[
M04] = Det4_0123_1234 * oneOverDet;
1206 m[
M10] = Det4_1234_0234 * mn1OverDet;
1207 m[
M11] = Det4_0234_0234 * oneOverDet;
1208 m[
M12] = Det4_0134_0234 * mn1OverDet;
1209 m[
M13] = Det4_0124_0234 * oneOverDet;
1210 m[
M14] = Det4_0123_0234 * mn1OverDet;
1212 m[
M20] = Det4_1234_0134 * oneOverDet;
1213 m[
M21] = Det4_0234_0134 * mn1OverDet;
1214 m[
M22] = Det4_0134_0134 * oneOverDet;
1215 m[
M23] = Det4_0124_0134 * mn1OverDet;
1216 m[
M24] = Det4_0123_0134 * oneOverDet;
1218 m[
M30] = Det4_1234_0124 * mn1OverDet;
1219 m[
M31] = Det4_0234_0124 * oneOverDet;
1220 m[
M32] = Det4_0134_0124 * mn1OverDet;
1221 m[
M33] = Det4_0124_0124 * oneOverDet;
1222 m[
M34] = Det4_0123_0124 * mn1OverDet;
1224 m[
M40] = Det4_1234_0123 * oneOverDet;
1225 m[
M41] = Det4_0234_0123 * mn1OverDet;
1226 m[
M42] = Det4_0134_0123 * oneOverDet;
1227 m[
M43] = Det4_0124_0123 * mn1OverDet;
1228 m[
M44] = Det4_0123_0123 * oneOverDet;
1288 + m[
A22]*Det2_34_01;
1290 + m[
A23]*Det2_34_01;
1292 + m[
A24]*Det2_34_01;
1294 + m[
A25]*Det2_34_01;
1296 + m[
A23]*Det2_34_02;
1298 + m[
A24]*Det2_34_02;
1300 + m[
A25]*Det2_34_02;
1302 + m[
A24]*Det2_34_03;
1304 + m[
A25]*Det2_34_03;
1306 + m[
A25]*Det2_34_04;
1308 + m[
A23]*Det2_34_12;
1310 + m[
A24]*Det2_34_12;
1312 + m[
A25]*Det2_34_12;
1314 + m[
A24]*Det2_34_13;
1316 + m[
A25]*Det2_34_13;
1318 + m[
A25]*Det2_34_14;
1320 + m[
A24]*Det2_34_23;
1322 + m[
A25]*Det2_34_23;
1324 + m[
A25]*Det2_34_24;
1326 + m[
A25]*Det2_34_34;
1328 + m[
A22]*Det2_35_01;
1330 + m[
A23]*Det2_35_01;
1332 + m[
A24]*Det2_35_01;
1334 + m[
A25]*Det2_35_01;
1336 + m[
A23]*Det2_35_02;
1338 + m[
A24]*Det2_35_02;
1340 + m[
A25]*Det2_35_02;
1342 + m[
A24]*Det2_35_03;
1344 + m[
A25]*Det2_35_03;
1346 + m[
A25]*Det2_35_04;
1348 + m[
A23]*Det2_35_12;
1350 + m[
A24]*Det2_35_12;
1352 + m[
A25]*Det2_35_12;
1354 + m[
A24]*Det2_35_13;
1356 + m[
A25]*Det2_35_13;
1358 + m[
A25]*Det2_35_14;
1360 + m[
A24]*Det2_35_23;
1362 + m[
A25]*Det2_35_23;
1364 + m[
A25]*Det2_35_24;
1366 + m[
A25]*Det2_35_34;
1368 + m[
A22]*Det2_45_01;
1370 + m[
A23]*Det2_45_01;
1372 + m[
A24]*Det2_45_01;
1374 + m[
A25]*Det2_45_01;
1376 + m[
A23]*Det2_45_02;
1378 + m[
A24]*Det2_45_02;
1380 + m[
A25]*Det2_45_02;
1382 + m[
A24]*Det2_45_03;
1384 + m[
A25]*Det2_45_03;
1386 + m[
A25]*Det2_45_04;
1388 + m[
A23]*Det2_45_12;
1390 + m[
A24]*Det2_45_12;
1392 + m[
A25]*Det2_45_12;
1394 + m[
A24]*Det2_45_13;
1396 + m[
A25]*Det2_45_13;
1398 + m[
A25]*Det2_45_14;
1400 + m[
A24]*Det2_45_23;
1402 + m[
A25]*Det2_45_23;
1404 + m[
A25]*Det2_45_24;
1406 + m[
A25]*Det2_45_34;
1408 + m[
A32]*Det2_45_01;
1410 + m[
A33]*Det2_45_01;
1412 + m[
A34]*Det2_45_01;
1414 + m[
A35]*Det2_45_01;
1416 + m[
A33]*Det2_45_02;
1418 + m[
A34]*Det2_45_02;
1420 + m[
A35]*Det2_45_02;
1422 + m[
A34]*Det2_45_03;
1424 + m[
A35]*Det2_45_03;
1426 + m[
A35]*Det2_45_04;
1428 + m[
A33]*Det2_45_12;
1430 + m[
A34]*Det2_45_12;
1432 + m[
A35]*Det2_45_12;
1434 + m[
A34]*Det2_45_13;
1436 + m[
A35]*Det2_45_13;
1438 + m[
A35]*Det2_45_14;
1440 + m[
A34]*Det2_45_23;
1442 + m[
A35]*Det2_45_23;
1444 + m[
A35]*Det2_45_24;
1446 + m[
A35]*Det2_45_34;
1450 G4double Det4_1234_0123 = m[
A10]*Det3_234_123 - m[
A11]*Det3_234_023
1451 + m[
A12]*Det3_234_013 - m[
A13]*Det3_234_012;
1452 G4double Det4_1234_0124 = m[
A10]*Det3_234_124 - m[
A11]*Det3_234_024
1453 + m[
A12]*Det3_234_014 - m[
A14]*Det3_234_012;
1454 G4double Det4_1234_0125 = m[
A10]*Det3_234_125 - m[
A11]*Det3_234_025
1455 + m[
A12]*Det3_234_015 - m[
A15]*Det3_234_012;
1456 G4double Det4_1234_0134 = m[
A10]*Det3_234_134 - m[
A11]*Det3_234_034
1457 + m[
A13]*Det3_234_014 - m[
A14]*Det3_234_013;
1458 G4double Det4_1234_0135 = m[
A10]*Det3_234_135 - m[
A11]*Det3_234_035
1459 + m[
A13]*Det3_234_015 - m[
A15]*Det3_234_013;
1460 G4double Det4_1234_0145 = m[
A10]*Det3_234_145 - m[
A11]*Det3_234_045
1461 + m[
A14]*Det3_234_015 - m[
A15]*Det3_234_014;
1462 G4double Det4_1234_0234 = m[
A10]*Det3_234_234 - m[
A12]*Det3_234_034
1463 + m[
A13]*Det3_234_024 - m[
A14]*Det3_234_023;
1464 G4double Det4_1234_0235 = m[
A10]*Det3_234_235 - m[
A12]*Det3_234_035
1465 + m[
A13]*Det3_234_025 - m[
A15]*Det3_234_023;
1466 G4double Det4_1234_0245 = m[
A10]*Det3_234_245 - m[
A12]*Det3_234_045
1467 + m[
A14]*Det3_234_025 - m[
A15]*Det3_234_024;
1468 G4double Det4_1234_0345 = m[
A10]*Det3_234_345 - m[
A13]*Det3_234_045
1469 + m[
A14]*Det3_234_035 - m[
A15]*Det3_234_034;
1470 G4double Det4_1234_1234 = m[
A11]*Det3_234_234 - m[
A12]*Det3_234_134
1471 + m[
A13]*Det3_234_124 - m[
A14]*Det3_234_123;
1472 G4double Det4_1234_1235 = m[
A11]*Det3_234_235 - m[
A12]*Det3_234_135
1473 + m[
A13]*Det3_234_125 - m[
A15]*Det3_234_123;
1474 G4double Det4_1234_1245 = m[
A11]*Det3_234_245 - m[
A12]*Det3_234_145
1475 + m[
A14]*Det3_234_125 - m[
A15]*Det3_234_124;
1476 G4double Det4_1234_1345 = m[
A11]*Det3_234_345 - m[
A13]*Det3_234_145
1477 + m[
A14]*Det3_234_135 - m[
A15]*Det3_234_134;
1478 G4double Det4_1234_2345 = m[
A12]*Det3_234_345 - m[
A13]*Det3_234_245
1479 + m[
A14]*Det3_234_235 - m[
A15]*Det3_234_234;
1480 G4double Det4_1235_0123 = m[
A10]*Det3_235_123 - m[
A11]*Det3_235_023
1481 + m[
A12]*Det3_235_013 - m[
A13]*Det3_235_012;
1482 G4double Det4_1235_0124 = m[
A10]*Det3_235_124 - m[
A11]*Det3_235_024
1483 + m[
A12]*Det3_235_014 - m[
A14]*Det3_235_012;
1484 G4double Det4_1235_0125 = m[
A10]*Det3_235_125 - m[
A11]*Det3_235_025
1485 + m[
A12]*Det3_235_015 - m[
A15]*Det3_235_012;
1486 G4double Det4_1235_0134 = m[
A10]*Det3_235_134 - m[
A11]*Det3_235_034
1487 + m[
A13]*Det3_235_014 - m[
A14]*Det3_235_013;
1488 G4double Det4_1235_0135 = m[
A10]*Det3_235_135 - m[
A11]*Det3_235_035
1489 + m[
A13]*Det3_235_015 - m[
A15]*Det3_235_013;
1490 G4double Det4_1235_0145 = m[
A10]*Det3_235_145 - m[
A11]*Det3_235_045
1491 + m[
A14]*Det3_235_015 - m[
A15]*Det3_235_014;
1492 G4double Det4_1235_0234 = m[
A10]*Det3_235_234 - m[
A12]*Det3_235_034
1493 + m[
A13]*Det3_235_024 - m[
A14]*Det3_235_023;
1494 G4double Det4_1235_0235 = m[
A10]*Det3_235_235 - m[
A12]*Det3_235_035
1495 + m[
A13]*Det3_235_025 - m[
A15]*Det3_235_023;
1496 G4double Det4_1235_0245 = m[
A10]*Det3_235_245 - m[
A12]*Det3_235_045
1497 + m[
A14]*Det3_235_025 - m[
A15]*Det3_235_024;
1498 G4double Det4_1235_0345 = m[
A10]*Det3_235_345 - m[
A13]*Det3_235_045
1499 + m[
A14]*Det3_235_035 - m[
A15]*Det3_235_034;
1500 G4double Det4_1235_1234 = m[
A11]*Det3_235_234 - m[
A12]*Det3_235_134
1501 + m[
A13]*Det3_235_124 - m[
A14]*Det3_235_123;
1502 G4double Det4_1235_1235 = m[
A11]*Det3_235_235 - m[
A12]*Det3_235_135
1503 + m[
A13]*Det3_235_125 - m[
A15]*Det3_235_123;
1504 G4double Det4_1235_1245 = m[
A11]*Det3_235_245 - m[
A12]*Det3_235_145
1505 + m[
A14]*Det3_235_125 - m[
A15]*Det3_235_124;
1506 G4double Det4_1235_1345 = m[
A11]*Det3_235_345 - m[
A13]*Det3_235_145
1507 + m[
A14]*Det3_235_135 - m[
A15]*Det3_235_134;
1508 G4double Det4_1235_2345 = m[
A12]*Det3_235_345 - m[
A13]*Det3_235_245
1509 + m[
A14]*Det3_235_235 - m[
A15]*Det3_235_234;
1510 G4double Det4_1245_0123 = m[
A10]*Det3_245_123 - m[
A11]*Det3_245_023
1511 + m[
A12]*Det3_245_013 - m[
A13]*Det3_245_012;
1512 G4double Det4_1245_0124 = m[
A10]*Det3_245_124 - m[
A11]*Det3_245_024
1513 + m[
A12]*Det3_245_014 - m[
A14]*Det3_245_012;
1514 G4double Det4_1245_0125 = m[
A10]*Det3_245_125 - m[
A11]*Det3_245_025
1515 + m[
A12]*Det3_245_015 - m[
A15]*Det3_245_012;
1516 G4double Det4_1245_0134 = m[
A10]*Det3_245_134 - m[
A11]*Det3_245_034
1517 + m[
A13]*Det3_245_014 - m[
A14]*Det3_245_013;
1518 G4double Det4_1245_0135 = m[
A10]*Det3_245_135 - m[
A11]*Det3_245_035
1519 + m[
A13]*Det3_245_015 - m[
A15]*Det3_245_013;
1520 G4double Det4_1245_0145 = m[
A10]*Det3_245_145 - m[
A11]*Det3_245_045
1521 + m[
A14]*Det3_245_015 - m[
A15]*Det3_245_014;
1522 G4double Det4_1245_0234 = m[
A10]*Det3_245_234 - m[
A12]*Det3_245_034
1523 + m[
A13]*Det3_245_024 - m[
A14]*Det3_245_023;
1524 G4double Det4_1245_0235 = m[
A10]*Det3_245_235 - m[
A12]*Det3_245_035
1525 + m[
A13]*Det3_245_025 - m[
A15]*Det3_245_023;
1526 G4double Det4_1245_0245 = m[
A10]*Det3_245_245 - m[
A12]*Det3_245_045
1527 + m[
A14]*Det3_245_025 - m[
A15]*Det3_245_024;
1528 G4double Det4_1245_0345 = m[
A10]*Det3_245_345 - m[
A13]*Det3_245_045
1529 + m[
A14]*Det3_245_035 - m[
A15]*Det3_245_034;
1530 G4double Det4_1245_1234 = m[
A11]*Det3_245_234 - m[
A12]*Det3_245_134
1531 + m[
A13]*Det3_245_124 - m[
A14]*Det3_245_123;
1532 G4double Det4_1245_1235 = m[
A11]*Det3_245_235 - m[
A12]*Det3_245_135
1533 + m[
A13]*Det3_245_125 - m[
A15]*Det3_245_123;
1534 G4double Det4_1245_1245 = m[
A11]*Det3_245_245 - m[
A12]*Det3_245_145
1535 + m[
A14]*Det3_245_125 - m[
A15]*Det3_245_124;
1536 G4double Det4_1245_1345 = m[
A11]*Det3_245_345 - m[
A13]*Det3_245_145
1537 + m[
A14]*Det3_245_135 - m[
A15]*Det3_245_134;
1538 G4double Det4_1245_2345 = m[
A12]*Det3_245_345 - m[
A13]*Det3_245_245
1539 + m[
A14]*Det3_245_235 - m[
A15]*Det3_245_234;
1540 G4double Det4_1345_0123 = m[
A10]*Det3_345_123 - m[
A11]*Det3_345_023
1541 + m[
A12]*Det3_345_013 - m[
A13]*Det3_345_012;
1542 G4double Det4_1345_0124 = m[
A10]*Det3_345_124 - m[
A11]*Det3_345_024
1543 + m[
A12]*Det3_345_014 - m[
A14]*Det3_345_012;
1544 G4double Det4_1345_0125 = m[
A10]*Det3_345_125 - m[
A11]*Det3_345_025
1545 + m[
A12]*Det3_345_015 - m[
A15]*Det3_345_012;
1546 G4double Det4_1345_0134 = m[
A10]*Det3_345_134 - m[
A11]*Det3_345_034
1547 + m[
A13]*Det3_345_014 - m[
A14]*Det3_345_013;
1548 G4double Det4_1345_0135 = m[
A10]*Det3_345_135 - m[
A11]*Det3_345_035
1549 + m[
A13]*Det3_345_015 - m[
A15]*Det3_345_013;
1550 G4double Det4_1345_0145 = m[
A10]*Det3_345_145 - m[
A11]*Det3_345_045
1551 + m[
A14]*Det3_345_015 - m[
A15]*Det3_345_014;
1552 G4double Det4_1345_0234 = m[
A10]*Det3_345_234 - m[
A12]*Det3_345_034
1553 + m[
A13]*Det3_345_024 - m[
A14]*Det3_345_023;
1554 G4double Det4_1345_0235 = m[
A10]*Det3_345_235 - m[
A12]*Det3_345_035
1555 + m[
A13]*Det3_345_025 - m[
A15]*Det3_345_023;
1556 G4double Det4_1345_0245 = m[
A10]*Det3_345_245 - m[
A12]*Det3_345_045
1557 + m[
A14]*Det3_345_025 - m[
A15]*Det3_345_024;
1558 G4double Det4_1345_0345 = m[
A10]*Det3_345_345 - m[
A13]*Det3_345_045
1559 + m[
A14]*Det3_345_035 - m[
A15]*Det3_345_034;
1560 G4double Det4_1345_1234 = m[
A11]*Det3_345_234 - m[
A12]*Det3_345_134
1561 + m[
A13]*Det3_345_124 - m[
A14]*Det3_345_123;
1562 G4double Det4_1345_1235 = m[
A11]*Det3_345_235 - m[
A12]*Det3_345_135
1563 + m[
A13]*Det3_345_125 - m[
A15]*Det3_345_123;
1564 G4double Det4_1345_1245 = m[
A11]*Det3_345_245 - m[
A12]*Det3_345_145
1565 + m[
A14]*Det3_345_125 - m[
A15]*Det3_345_124;
1566 G4double Det4_1345_1345 = m[
A11]*Det3_345_345 - m[
A13]*Det3_345_145
1567 + m[
A14]*Det3_345_135 - m[
A15]*Det3_345_134;
1568 G4double Det4_1345_2345 = m[
A12]*Det3_345_345 - m[
A13]*Det3_345_245
1569 + m[
A14]*Det3_345_235 - m[
A15]*Det3_345_234;
1570 G4double Det4_2345_0123 = m[
A20]*Det3_345_123 - m[
A21]*Det3_345_023
1571 + m[
A22]*Det3_345_013 - m[
A23]*Det3_345_012;
1572 G4double Det4_2345_0124 = m[
A20]*Det3_345_124 - m[
A21]*Det3_345_024
1573 + m[
A22]*Det3_345_014 - m[
A24]*Det3_345_012;
1574 G4double Det4_2345_0125 = m[
A20]*Det3_345_125 - m[
A21]*Det3_345_025
1575 + m[
A22]*Det3_345_015 - m[
A25]*Det3_345_012;
1576 G4double Det4_2345_0134 = m[
A20]*Det3_345_134 - m[
A21]*Det3_345_034
1577 + m[
A23]*Det3_345_014 - m[
A24]*Det3_345_013;
1578 G4double Det4_2345_0135 = m[
A20]*Det3_345_135 - m[
A21]*Det3_345_035
1579 + m[
A23]*Det3_345_015 - m[
A25]*Det3_345_013;
1580 G4double Det4_2345_0145 = m[
A20]*Det3_345_145 - m[
A21]*Det3_345_045
1581 + m[
A24]*Det3_345_015 - m[
A25]*Det3_345_014;
1582 G4double Det4_2345_0234 = m[
A20]*Det3_345_234 - m[
A22]*Det3_345_034
1583 + m[
A23]*Det3_345_024 - m[
A24]*Det3_345_023;
1584 G4double Det4_2345_0235 = m[
A20]*Det3_345_235 - m[
A22]*Det3_345_035
1585 + m[
A23]*Det3_345_025 - m[
A25]*Det3_345_023;
1586 G4double Det4_2345_0245 = m[
A20]*Det3_345_245 - m[
A22]*Det3_345_045
1587 + m[
A24]*Det3_345_025 - m[
A25]*Det3_345_024;
1588 G4double Det4_2345_0345 = m[
A20]*Det3_345_345 - m[
A23]*Det3_345_045
1589 + m[
A24]*Det3_345_035 - m[
A25]*Det3_345_034;
1590 G4double Det4_2345_1234 = m[
A21]*Det3_345_234 - m[
A22]*Det3_345_134
1591 + m[
A23]*Det3_345_124 - m[
A24]*Det3_345_123;
1592 G4double Det4_2345_1235 = m[
A21]*Det3_345_235 - m[
A22]*Det3_345_135
1593 + m[
A23]*Det3_345_125 - m[
A25]*Det3_345_123;
1594 G4double Det4_2345_1245 = m[
A21]*Det3_345_245 - m[
A22]*Det3_345_145
1595 + m[
A24]*Det3_345_125 - m[
A25]*Det3_345_124;
1596 G4double Det4_2345_1345 = m[
A21]*Det3_345_345 - m[
A23]*Det3_345_145
1597 + m[
A24]*Det3_345_135 - m[
A25]*Det3_345_134;
1598 G4double Det4_2345_2345 = m[
A22]*Det3_345_345 - m[
A23]*Det3_345_245
1599 + m[
A24]*Det3_345_235 - m[
A25]*Det3_345_234;
1603 G4double Det5_01234_01234 = m[
A00]*Det4_1234_1234 - m[
A01]*Det4_1234_0234
1604 + m[
A02]*Det4_1234_0134 - m[
A03]*Det4_1234_0124 + m[
A04]*Det4_1234_0123;
1605 G4double Det5_01234_01235 = m[
A00]*Det4_1234_1235 - m[
A01]*Det4_1234_0235
1606 + m[
A02]*Det4_1234_0135 - m[
A03]*Det4_1234_0125 + m[
A05]*Det4_1234_0123;
1607 G4double Det5_01234_01245 = m[
A00]*Det4_1234_1245 - m[
A01]*Det4_1234_0245
1608 + m[
A02]*Det4_1234_0145 - m[
A04]*Det4_1234_0125 + m[
A05]*Det4_1234_0124;
1609 G4double Det5_01234_01345 = m[
A00]*Det4_1234_1345 - m[
A01]*Det4_1234_0345
1610 + m[
A03]*Det4_1234_0145 - m[
A04]*Det4_1234_0135 + m[
A05]*Det4_1234_0134;
1611 G4double Det5_01234_02345 = m[
A00]*Det4_1234_2345 - m[
A02]*Det4_1234_0345
1612 + m[
A03]*Det4_1234_0245 - m[
A04]*Det4_1234_0235 + m[
A05]*Det4_1234_0234;
1613 G4double Det5_01234_12345 = m[
A01]*Det4_1234_2345 - m[
A02]*Det4_1234_1345
1614 + m[
A03]*Det4_1234_1245 - m[
A04]*Det4_1234_1235 + m[
A05]*Det4_1234_1234;
1615 G4double Det5_01235_01234 = m[
A00]*Det4_1235_1234 - m[
A01]*Det4_1235_0234
1616 + m[
A02]*Det4_1235_0134 - m[
A03]*Det4_1235_0124 + m[
A04]*Det4_1235_0123;
1617 G4double Det5_01235_01235 = m[
A00]*Det4_1235_1235 - m[
A01]*Det4_1235_0235
1618 + m[
A02]*Det4_1235_0135 - m[
A03]*Det4_1235_0125 + m[
A05]*Det4_1235_0123;
1619 G4double Det5_01235_01245 = m[
A00]*Det4_1235_1245 - m[
A01]*Det4_1235_0245
1620 + m[
A02]*Det4_1235_0145 - m[
A04]*Det4_1235_0125 + m[
A05]*Det4_1235_0124;
1621 G4double Det5_01235_01345 = m[
A00]*Det4_1235_1345 - m[
A01]*Det4_1235_0345
1622 + m[
A03]*Det4_1235_0145 - m[
A04]*Det4_1235_0135 + m[
A05]*Det4_1235_0134;
1623 G4double Det5_01235_02345 = m[
A00]*Det4_1235_2345 - m[
A02]*Det4_1235_0345
1624 + m[
A03]*Det4_1235_0245 - m[
A04]*Det4_1235_0235 + m[
A05]*Det4_1235_0234;
1625 G4double Det5_01235_12345 = m[
A01]*Det4_1235_2345 - m[
A02]*Det4_1235_1345
1626 + m[
A03]*Det4_1235_1245 - m[
A04]*Det4_1235_1235 + m[
A05]*Det4_1235_1234;
1627 G4double Det5_01245_01234 = m[
A00]*Det4_1245_1234 - m[
A01]*Det4_1245_0234
1628 + m[
A02]*Det4_1245_0134 - m[
A03]*Det4_1245_0124 + m[
A04]*Det4_1245_0123;
1629 G4double Det5_01245_01235 = m[
A00]*Det4_1245_1235 - m[
A01]*Det4_1245_0235
1630 + m[
A02]*Det4_1245_0135 - m[
A03]*Det4_1245_0125 + m[
A05]*Det4_1245_0123;
1631 G4double Det5_01245_01245 = m[
A00]*Det4_1245_1245 - m[
A01]*Det4_1245_0245
1632 + m[
A02]*Det4_1245_0145 - m[
A04]*Det4_1245_0125 + m[
A05]*Det4_1245_0124;
1633 G4double Det5_01245_01345 = m[
A00]*Det4_1245_1345 - m[
A01]*Det4_1245_0345
1634 + m[
A03]*Det4_1245_0145 - m[
A04]*Det4_1245_0135 + m[
A05]*Det4_1245_0134;
1635 G4double Det5_01245_02345 = m[
A00]*Det4_1245_2345 - m[
A02]*Det4_1245_0345
1636 + m[
A03]*Det4_1245_0245 - m[
A04]*Det4_1245_0235 + m[
A05]*Det4_1245_0234;
1637 G4double Det5_01245_12345 = m[
A01]*Det4_1245_2345 - m[
A02]*Det4_1245_1345
1638 + m[
A03]*Det4_1245_1245 - m[
A04]*Det4_1245_1235 + m[
A05]*Det4_1245_1234;
1639 G4double Det5_01345_01234 = m[
A00]*Det4_1345_1234 - m[
A01]*Det4_1345_0234
1640 + m[
A02]*Det4_1345_0134 - m[
A03]*Det4_1345_0124 + m[
A04]*Det4_1345_0123;
1641 G4double Det5_01345_01235 = m[
A00]*Det4_1345_1235 - m[
A01]*Det4_1345_0235
1642 + m[
A02]*Det4_1345_0135 - m[
A03]*Det4_1345_0125 + m[
A05]*Det4_1345_0123;
1643 G4double Det5_01345_01245 = m[
A00]*Det4_1345_1245 - m[
A01]*Det4_1345_0245
1644 + m[
A02]*Det4_1345_0145 - m[
A04]*Det4_1345_0125 + m[
A05]*Det4_1345_0124;
1645 G4double Det5_01345_01345 = m[
A00]*Det4_1345_1345 - m[
A01]*Det4_1345_0345
1646 + m[
A03]*Det4_1345_0145 - m[
A04]*Det4_1345_0135 + m[
A05]*Det4_1345_0134;
1647 G4double Det5_01345_02345 = m[
A00]*Det4_1345_2345 - m[
A02]*Det4_1345_0345
1648 + m[
A03]*Det4_1345_0245 - m[
A04]*Det4_1345_0235 + m[
A05]*Det4_1345_0234;
1649 G4double Det5_01345_12345 = m[
A01]*Det4_1345_2345 - m[
A02]*Det4_1345_1345
1650 + m[
A03]*Det4_1345_1245 - m[
A04]*Det4_1345_1235 + m[
A05]*Det4_1345_1234;
1651 G4double Det5_02345_01234 = m[
A00]*Det4_2345_1234 - m[
A01]*Det4_2345_0234
1652 + m[
A02]*Det4_2345_0134 - m[
A03]*Det4_2345_0124 + m[
A04]*Det4_2345_0123;
1653 G4double Det5_02345_01235 = m[
A00]*Det4_2345_1235 - m[
A01]*Det4_2345_0235
1654 + m[
A02]*Det4_2345_0135 - m[
A03]*Det4_2345_0125 + m[
A05]*Det4_2345_0123;
1655 G4double Det5_02345_01245 = m[
A00]*Det4_2345_1245 - m[
A01]*Det4_2345_0245
1656 + m[
A02]*Det4_2345_0145 - m[
A04]*Det4_2345_0125 + m[
A05]*Det4_2345_0124;
1657 G4double Det5_02345_01345 = m[
A00]*Det4_2345_1345 - m[
A01]*Det4_2345_0345
1658 + m[
A03]*Det4_2345_0145 - m[
A04]*Det4_2345_0135 + m[
A05]*Det4_2345_0134;
1659 G4double Det5_02345_02345 = m[
A00]*Det4_2345_2345 - m[
A02]*Det4_2345_0345
1660 + m[
A03]*Det4_2345_0245 - m[
A04]*Det4_2345_0235 + m[
A05]*Det4_2345_0234;
1661 G4double Det5_02345_12345 = m[
A01]*Det4_2345_2345 - m[
A02]*Det4_2345_1345
1662 + m[
A03]*Det4_2345_1245 - m[
A04]*Det4_2345_1235 + m[
A05]*Det4_2345_1234;
1663 G4double Det5_12345_01234 = m[
A10]*Det4_2345_1234 - m[
A11]*Det4_2345_0234
1664 + m[
A12]*Det4_2345_0134 - m[
A13]*Det4_2345_0124 + m[
A14]*Det4_2345_0123;
1665 G4double Det5_12345_01235 = m[
A10]*Det4_2345_1235 - m[
A11]*Det4_2345_0235
1666 + m[
A12]*Det4_2345_0135 - m[
A13]*Det4_2345_0125 + m[
A15]*Det4_2345_0123;
1667 G4double Det5_12345_01245 = m[
A10]*Det4_2345_1245 - m[
A11]*Det4_2345_0245
1668 + m[
A12]*Det4_2345_0145 - m[
A14]*Det4_2345_0125 + m[
A15]*Det4_2345_0124;
1669 G4double Det5_12345_01345 = m[
A10]*Det4_2345_1345 - m[
A11]*Det4_2345_0345
1670 + m[
A13]*Det4_2345_0145 - m[
A14]*Det4_2345_0135 + m[
A15]*Det4_2345_0134;
1671 G4double Det5_12345_02345 = m[
A10]*Det4_2345_2345 - m[
A12]*Det4_2345_0345
1672 + m[
A13]*Det4_2345_0245 - m[
A14]*Det4_2345_0235 + m[
A15]*Det4_2345_0234;
1673 G4double Det5_12345_12345 = m[
A11]*Det4_2345_2345 - m[
A12]*Det4_2345_1345
1674 + m[
A13]*Det4_2345_1245 - m[
A14]*Det4_2345_1235 + m[
A15]*Det4_2345_1234;
1679 - m[
A01]*Det5_12345_02345
1680 + m[
A02]*Det5_12345_01345
1681 - m[
A03]*Det5_12345_01245
1682 + m[
A04]*Det5_12345_01235
1683 - m[
A05]*Det5_12345_01234;
1692 G4double mn1OverDet = - oneOverDet;
1694 m[
A00] = Det5_12345_12345*oneOverDet;
1695 m[
A01] = Det5_02345_12345*mn1OverDet;
1696 m[
A02] = Det5_01345_12345*oneOverDet;
1697 m[
A03] = Det5_01245_12345*mn1OverDet;
1698 m[
A04] = Det5_01235_12345*oneOverDet;
1699 m[
A05] = Det5_01234_12345*mn1OverDet;
1701 m[
A10] = Det5_12345_02345*mn1OverDet;
1702 m[
A11] = Det5_02345_02345*oneOverDet;
1703 m[
A12] = Det5_01345_02345*mn1OverDet;
1704 m[
A13] = Det5_01245_02345*oneOverDet;
1705 m[
A14] = Det5_01235_02345*mn1OverDet;
1706 m[
A15] = Det5_01234_02345*oneOverDet;
1708 m[
A20] = Det5_12345_01345*oneOverDet;
1709 m[
A21] = Det5_02345_01345*mn1OverDet;
1710 m[
A22] = Det5_01345_01345*oneOverDet;
1711 m[
A23] = Det5_01245_01345*mn1OverDet;
1712 m[
A24] = Det5_01235_01345*oneOverDet;
1713 m[
A25] = Det5_01234_01345*mn1OverDet;
1715 m[
A30] = Det5_12345_01245*mn1OverDet;
1716 m[
A31] = Det5_02345_01245*oneOverDet;
1717 m[
A32] = Det5_01345_01245*mn1OverDet;
1718 m[
A33] = Det5_01245_01245*oneOverDet;
1719 m[
A34] = Det5_01235_01245*mn1OverDet;
1720 m[
A35] = Det5_01234_01245*oneOverDet;
1722 m[
A40] = Det5_12345_01235*oneOverDet;
1723 m[
A41] = Det5_02345_01235*mn1OverDet;
1724 m[
A42] = Det5_01345_01235*oneOverDet;
1725 m[
A43] = Det5_01245_01235*mn1OverDet;
1726 m[
A44] = Det5_01235_01235*oneOverDet;
1727 m[
A45] = Det5_01234_01235*mn1OverDet;
1729 m[
A50] = Det5_12345_01234*mn1OverDet;
1730 m[
A51] = Det5_02345_01234*oneOverDet;
1731 m[
A52] = Det5_01345_01234*mn1OverDet;
1732 m[
A53] = Det5_01245_01234*oneOverDet;
1733 m[
A54] = Det5_01235_01234*mn1OverDet;
1734 m[
A55] = Det5_01234_01234*oneOverDet;
std::vector< ExP01TrackerHit * > a
virtual G4int num_row() const
std::ostream & operator<<(std::ostream &, const BasicVector3D< float > &)
G4ErrorMatrix & operator=(const G4ErrorMatrix &m2)
G4ErrorMatrix & operator*=(G4double t)
void message(RunManager *runmanager)
virtual void invertHaywood6(G4int &ierr)
G4ErrorMatrix sub(G4int min_row, G4int max_row, G4int min_col, G4int max_col) const
static double normal(HepRandomEngine *eptr)
virtual void invertHaywood5(G4int &ierr)
BasicVector3D< float > operator*(const BasicVector3D< float > &v, double a)
static void error(const char *s)
static constexpr double m
std::vector< G4double >::const_iterator G4ErrorMatrixConstIter
G4int dfinv_matrix(G4int *ir)
virtual void invertHaywood4(G4int &ierr)
G4ErrorMatrix dsum(const G4ErrorMatrix &, const G4ErrorMatrix &)
G4ErrorMatrix apply(G4double(*f)(G4double, G4int, G4int)) const
BasicVector3D< float > operator+(const BasicVector3D< float > &v)
std::vector< G4double >::iterator G4ErrorMatrixIter
std::vector< G4double > m
double epsilon(double density, double temperature)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
G4double determinant() const
#define CHK_DIM_2(r1, r2, c1, c2, fun)
std::vector< G4double > m
G4ErrorMatrix & operator/=(G4double t)
G4ErrorMatrix operator-() const
G4ErrorMatrix & operator+=(const G4ErrorMatrix &m2)
BasicVector3D< float > operator/(const BasicVector3D< float > &v, double a)
virtual G4int num_col() const
G4ErrorMatrix & operator-=(const G4ErrorMatrix &m2)
#define CHK_DIM_1(c1, r2, fun)
G4int dfact_matrix(G4double &det, G4int *ir)
virtual void invert(G4int &ierr)
BasicVector3D< float > operator-(const BasicVector3D< float > &v)