void ScalarFiniteElement<D> :: CalcMappedDShape (const MappedIntegrationPoint<D,D> & mip, SliceMatrix<> dshape) const { CalcDShape (mip.IP(), dshape); for (int i = 0; i < dshape.Height(); i++) { Vec<D> hv = dshape.Row(i); FlatVec<D> (&dshape(i,0)) = Trans (mip.GetJacobianInverse ()) * hv; } }
inline void MultMatMat (SliceMatrix<> a, SliceMatrix<SIMD<double>> b, SliceMatrix<SIMD<double>> c) { MultMatMat_intern (a.Height(), a.Width(), b.Width(), a, b, c); }