From: Wolfgang Bangerth Date: Thu, 10 May 2001 13:58:27 +0000 (+0000) Subject: Revert order in which matrix is stored (requires vmult instead of X-Git-Tag: v8.0.0~19205 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=c1ba99dd99fe98830e22646ff64d39b78a4c95ff;p=dealii.git Revert order in which matrix is stored (requires vmult instead of Tvmult then). git-svn-id: https://svn.dealii.org/trunk@4586 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/dofs/dof_tools.cc b/deal.II/deal.II/source/dofs/dof_tools.cc index 97590cf689..f505f86190 100644 --- a/deal.II/deal.II/source/dofs/dof_tools.cc +++ b/deal.II/deal.II/source/dofs/dof_tools.cc @@ -1391,8 +1391,6 @@ compute_intergrid_transfer_representation (const DoFHandler &c const InterGridMap &coarse_to_fine_grid_map, std::vector > &transfer_representation) { - - // store the weights with which a dof // on the parameter grid contributes // to a dof on the fine grid. see the @@ -1446,9 +1444,6 @@ compute_intergrid_transfer_representation (const DoFHandler &c // now compute the requested // representation - transfer_representation.clear (); - transfer_representation.resize (weights.size()); - const unsigned int n_global_parm_dofs = std::count_if (weight_mapping.begin(), weight_mapping.end(), std::bind2nd (std::not_equal_to (), -1)); @@ -1475,12 +1470,22 @@ compute_intergrid_transfer_representation (const DoFHandler &c // next copy over weights array // and replace respective // numbers + const unsigned int n_rows = weight_mapping.size(); + + transfer_representation.clear (); + transfer_representation.resize (n_rows); + const unsigned int n_coarse_dofs = coarse_grid.n_dofs(); for (unsigned int i=0; i::const_iterator j = weights[i].begin(); for (; j!=weights[i].end(); ++j) - transfer_representation[i][inverse_weight_mapping[j->first]] = j->second; + { + const unsigned int p = inverse_weight_mapping[j->first]; + Assert (psecond; + }; }; };