]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Use the correct potential. Do not scale the eigenvalue in the output.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Mon, 6 Jul 2009 17:05:49 +0000 (17:05 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Mon, 6 Jul 2009 17:05:49 +0000 (17:05 +0000)
git-svn-id: https://svn.dealii.org/trunk@19021 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-36/step-36.cc

index 980060b5f17bcf08b4111b688f27eb5c16dc89cf..d7b3d40bc13df8757622f9339d25ef59678708b7 100644 (file)
@@ -229,9 +229,8 @@ void EigenvalueProblem<dim>::assemble_system ()
                += (fe_values.shape_grad (i, q_point) *
                    fe_values.shape_grad (j, q_point) 
                    + 
+                   potential_values[q_point] *
                    fe_values.shape_value (i, q_point) *
-                   //               potential_values[q_point] *
-                   0 * // infinite potential well
                    fe_values.shape_value (j, q_point)
                    ) * fe_values.JxW (q_point);
              
@@ -377,7 +376,7 @@ void EigenvalueProblem<dim>::run ()
   for (unsigned int i=0; i<eigenvalues.size(); ++i)
     std::cout << std::endl 
              << "      eigenvalue " << i 
-             << " : " << eigenvalues[i] * 4 / 3.1415926 / 3.1415926;
+             << " : " << eigenvalues[i];
 }
 
 

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.