| 6 | { |
| 7 | template<typename VOLFE> |
| 8 | void GlobalInterfaceSpace::VolDiffOp<VOLFE>::CalcMatrix |
| 9 | (const FiniteElement & bfel, |
| 10 | const BaseMappedIntegrationPoint & mip, |
| 11 | BareSliceMatrix<double,ColMajor> mat, |
| 12 | LocalHeap & lh) const |
| 13 | { |
| 14 | auto & ip = mip.IP(); |
| 15 | auto & fel = dynamic_cast<const VOLFE&> (bfel); |
| 16 | mat.AddSize(Dim(), bfel.GetNDof()) = 0; |
| 17 | |
| 18 | int fnr = ip.FacetNr(); |
| 19 | |
| 20 | if (fnr != -1 && fel.GetFacet(fnr)) |
| 21 | { |
| 22 | fel.GetFacet(fnr)->CalcShape(mip, mat.Row(0)); |
| 23 | return; |
| 24 | } |
| 25 | if(fnr == -1) |
| 26 | { |
| 27 | fel.CalcShape(mip, mat.Row(0)); |
| 28 | } |
| 29 | } |
| 30 | |
| 31 | template<typename VOLFE> |
| 32 | void GlobalInterfaceSpace::ParameterGradDiffOp<VOLFE>::CalcMatrix |