| 1137 | |
| 1138 | template <typename MIR, typename T, ORDERING ORD> |
| 1139 | void T_Evaluate (const MIR & mir, |
| 1140 | BareSliceMatrix<T,ORD> result) const |
| 1141 | { |
| 1142 | int hd = Dimensions()[0]; |
| 1143 | c1->Evaluate (mir, result); |
| 1144 | STACK_ARRAY(T, hmem, hd*hd); |
| 1145 | FlatMatrix<T,ORD> tmp (hd, hd, &hmem[0]); |
| 1146 | |
| 1147 | for (size_t i = 0; i < mir.Size(); i++) |
| 1148 | { |
| 1149 | for (int j = 0; j < hd; j++) |
| 1150 | for (int k = 0; k < hd; k++) |
| 1151 | tmp(j,k) = result(j*hd+k, i); |
| 1152 | for (int j = 0; j < hd; j++) |
| 1153 | for (int k = 0; k < hd; k++) |
| 1154 | result(j*hd+k, i) = 0.5*(tmp(j,k)-tmp(k,j)); |
| 1155 | } |
| 1156 | } |
| 1157 | |
| 1158 | template <typename MIR, typename T, ORDERING ORD> |
| 1159 | void T_Evaluate (const MIR & ir, |
nothing calls this directly
no test coverage detected