]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Two fixes.
authorhartmann <hartmann@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 19 Mar 2001 09:54:06 +0000 (09:54 +0000)
committerhartmann <hartmann@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 19 Mar 2001 09:54:06 +0000 (09:54 +0000)
git-svn-id: https://svn.dealii.org/trunk@4241 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/fe/mapping_q.cc

index 5a887dd40982dd928d0b04946019275d53d76ab1..223effc7063042826354649c090b23f1c24a14a9 100644 (file)
@@ -442,11 +442,18 @@ MappingQ<dim>::set_laplace_on_quad_vector(std::vector<std::vector<double> > &loq
     }
   else
     {
-                                  // not precomputed, then do so now
+                                      // not precomputed, then do so now
       if (dim==2)
        compute_laplace_vector(loqvs);
-      else
-       Assert(false, ExcNotImplemented());
+      
+                                      // for dim==3 don't throw an
+                                      // ExcNotImplemented here to
+                                      // allow the creating of that
+                                      // MappingQ<3> object. But an
+                                      // ExcLaplaceVectorNotSet
+                                      // assertion is thrown when the
+                                      // apply_laplace_vector
+                                      // function is called.
     }
 
                                   // the sum of weights of the points
@@ -986,6 +993,16 @@ MappingQ<3>::add_quad_support_points(const Triangulation<3>::cell_iterator &cell
                                           // separately
          if (lines_at_boundary>0)
            {
+                                              // call of function
+                                              // apply_laplace_vector
+                                              // increases size of b
+                                              // about 1. There
+                                              // resize b for the
+                                              // case the mentioned
+                                              // function was already
+                                              // called.
+             b.resize(4*degree);
+             
                                               // b is of size
                                               // 4*degree, make sure
                                               // that this is the

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.