From 40d98d8b50a21b8058602602038caabc32ecc5d9 Mon Sep 17 00:00:00 2001 From: wolf Date: Mon, 12 Apr 1999 17:08:00 +0000 Subject: [PATCH] Reinsert computation of prolongation and restriction matrices for FESystem objects from the base classes' respective matrices. git-svn-id: https://svn.dealii.org/trunk@1125 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/fe/fe_system.cc | 32 ++++++++++++++------------ 1 file changed, 17 insertions(+), 15 deletions(-) diff --git a/deal.II/deal.II/source/fe/fe_system.cc b/deal.II/deal.II/source/fe/fe_system.cc index cd83811a90..329b252dc8 100644 --- a/deal.II/deal.II/source/fe/fe_system.cc +++ b/deal.II/deal.II/source/fe/fe_system.cc @@ -262,28 +262,30 @@ void FESystem::initialize () build_cell_table(); build_face_table(); // distribute the matrices of the base - // finite element to the matrices of + // finite elements to the matrices of // this object -/* for (unsigned int i=0; itotal_dofs; ++i) - for (unsigned int j=0; jtotal_dofs; ++j) - for (unsigned int n=0; n::children_per_cell; ++child) { - restriction[child] (i*n_sub_elements + n, - j*n_sub_elements + n) - = base_element->restrict(child)(i,j); - prolongation[child] (i*n_sub_elements + n, - j*n_sub_elements + n) - = base_element->prolongate(child)(i,j); + restriction[child] (component_to_system_index (component,i), + component_to_system_index (component, j)) + = base_element(component_to_base_table[component]).restrict(child)(i,j); + prolongation[child] (component_to_system_index (component,i), + component_to_system_index (component, j)) + = base_element(component_to_base_table[component]).prolongate(child)(i,j); }; - interface_constraints (i*n_sub_elements + n, - j*n_sub_elements + n) - = base_element->constraints()(i,j); +// still to do: +// interface_constraints (// component_to_system_index (component,i), +// component_to_system_index (component, j)) +// = base_element(component_to_base_table[component]).constraints()(i,j); }; -*/}; +}; -- 2.39.5