test

Author: ESeNonFossiIo

include/interfaces/navier_stokes-chorin-temam.h

@@ -296,12 +296,12 @@ energies_and_residuals(const typename DoFHandler<dim,spacedim>::active_cell_iter
           res += ((v-ue)/delta_t)*v_test;
           res +=  (ue * grad_v )*v_test;
           res += nu * scalar_product( sym_grad_v,
-                    sym_grad_v_test);
+                                      sym_grad_v_test);
 
           res +=  grad_p*grad_p_test + (1/delta_t)*(div_v * p_test);
 
 // Updating of the solution:
-          res += (u - v + delta_t * grad_p)*u_test;
+          res += (u - v + delta_t *grad_p)*u_test;
 
           residual[0][i] += res * JxW[quad];
         }
@@ -352,18 +352,20 @@ NavierStokes<dim,spacedim,LAC>::compute_system_operators(
 
   auto Au = linear_operator<BVEC>(matrices[0]->block(2,2));
   auto Au_inv  =  inverse_operator(Au, solver_GMRES, AMG_u);
-  
+
   // Preconditioner
   //////////////////////////////
 
   BlockLinearOperator<VEC> diag_inv
   = block_diagonal_operator<num_blocks, VEC>(
-    {{ 
+  {
+    {
       Av_inv, Ap_inv, Au_inv
-    }}
-    );
+    }
+  }
+  );
   prec_op = block_forward_substitution(
-              BlockLinearOperator< VEC >(S), 
+              BlockLinearOperator< VEC >(S),
               diag_inv);
 }
 

tests/navier_stokes-chorin-temam_01.cc

@@ -0,0 +1,42 @@
+#include "pidomus.h"
+#include "interfaces/navier_stokes-chorin-temam.h"
+#include "tests.h"
+
+/**
+ * Test:     Navier Stokes interface.
+ * Method:   Direct
+ * Problem:  Non time depending Navier Stokes Equations
+ * Exact solution:
+ * \f[
+ *    u=\big( 2*(x^2)*y, -2*x*(y^2) \big)
+ *    \textrm{ and }p=xy;
+ * \f]
+ */
+
+using namespace dealii;
+int main (int argc, char *argv[])
+{
+
+  Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv,
+                                                      numbers::invalid_unsigned_int);
+
+  initlog();
+  deallog.depth_file(1);
+  deallog.threshold_double(1.0e-3);
+
+  NavierStokes<2,2,LATrilinos> energy;
+  piDoMUS<2,2,LATrilinos> navier_stokes ("",energy);
+  ParameterAcceptor::initialize(
+    SOURCE_DIR "/parameters/navier_stokes_chorin-temam_01.prm",
+    "used_parameters.prm");
+
+  navier_stokes.run ();
+
+  auto sol = navier_stokes.get_solution();
+  for (unsigned int i = 0 ; i<sol.size(); ++i)
+    {
+      deallog << sol[i] << std::endl ;
+    }
+
+  return 0;
+}