From: wolf Date: Mon, 12 Apr 1999 17:08:00 +0000 (+0000) Subject: Reinsert computation of prolongation and restriction matrices for X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=40d98d8b50a21b8058602602038caabc32ecc5d9;p=dealii-svn.git 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 --- 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); }; -*/}; +};