Merge pull request #333 from arcaneframework/mab/petsc-improvements

BSR for PETSc

Author: mohd-afeef-badri

modules/testlab/CMakeLists.txt

@@ -85,6 +85,8 @@ if(FEMUTILS_HAS_SOLVER_BACKEND_PETSC)
   arcanefem_add_gpu_test(NAME [testlab]petsc_direct.bjacobi.gmres COMMAND  Testlab ARGS inputs/Test.petsc_direct.bjacobi.gmres.arc)
   arcanefem_add_gpu_test(NAME [testlab]petsc_direct.gamg COMMAND  Testlab ARGS inputs/Test.petsc_direct.gamg.arc)
   arcanefem_add_gpu_test(NAME [testlab]petsc_direct_3D COMMAND  Testlab ARGS inputs/Test.petsc_direct-3D.arc)
+  arcanefem_add_gpu_test(NAME [testlab]sphere_3D_petsc_direct_bsr COMMAND Testlab ARGS inputs/Test.sphere.3D.bsr.petsc_direct.arc)
+  arcanefem_add_gpu_test(NAME [testlab]sphere_3D_petsc_direct_af-bsr COMMAND Testlab ARGS inputs/Test.sphere.3D.bsr.atomic-free.petsc_direct.arc)
 
   arcanefem_add_gpu_test(NAME [testlab]L-shape_2D_petsc COMMAND ./Testlab ARGS inputs/Test.petsc.arc)
   arcanefem_add_gpu_test(NAME [testlab]L-shape_2D_petsc_coo-gpu COMMAND ./Testlab ARGS inputs/Test.L-shape.2D.coo-gpu.arc)

modules/testlab/FemModule.cc

@@ -138,7 +138,9 @@ startInit()
     }
 
     if (options()->bsr || options()->bsrAtomicFree()) {
-      bool use_csr_in_linear_system = options()->linearSystem.serviceName() == "HypreLinearSystem";
+      bool use_csr_in_linear_system = false;
+      if (options()->linearSystem.serviceName() == "HypreLinearSystem" || options()->linearSystem.serviceName() == "PETScLinearSystem")
+        use_csr_in_linear_system = true;
       m_bsr_format.initialize(mesh, 1, use_csr_in_linear_system, options()->bsrAtomicFree);
     }
   }

modules/testlab/inputs/Test.sphere.3D.bsr.atomic-free.petsc_direct.arc

@@ -0,0 +1,37 @@
+<?xml version="1.0"?>
+<case codename="Testlab" xml:lang="en" codeversion="1.0">
+  <arcane>
+    <title>Testlab: 3D Sphere with Horizontal Dirichlet Boundary (Atomic Free BSR-GPU Format, PETSc Linear System Solver, Penalty Method)</title>
+    <timeloop>TestlabLoop</timeloop>
+  </arcane>
+
+  <arcane-post-processing>
+   <output-period>1</output-period>
+   <format name="VtkHdfV2PostProcessor" />
+   <output>
+     <variable>U</variable>
+   </output>
+  </arcane-post-processing>
+
+  <meshes>
+    <mesh>
+      <filename>sphere_cut.msh</filename>
+    </mesh>
+  </meshes>
+
+  <fem>
+    <solution-comparison-file>poisson_test_ref_sphere_3D.txt</solution-comparison-file>
+    <f>5.5</f>
+    <enforce-Dirichlet-method>Penalty</enforce-Dirichlet-method>
+    <penalty>1.e31</penalty>
+    <dirichlet-boundary-condition>
+      <surface>horizontal</surface>
+      <value>0.5</value>
+    </dirichlet-boundary-condition>
+    <bsr-atomic-free>true</bsr-atomic-free>
+    <linear-system name="PETScLinearSystem">
+      <rtol>0.</rtol>
+      <atol>1e-5</atol>
+    </linear-system>
+  </fem>
+</case>

modules/testlab/inputs/Test.sphere.3D.bsr.petsc_direct.arc

@@ -0,0 +1,37 @@
+<?xml version="1.0"?>
+<case codename="Testlab" xml:lang="en" codeversion="1.0">
+  <arcane>
+    <title>Testlab: 3D Sphere with Horizontal Dirichlet Boundary (BSR Format, PETSc Linear System Solver, Penalty Method)</title>
+    <timeloop>TestlabLoop</timeloop>
+  </arcane>
+
+  <arcane-post-processing>
+   <output-period>1</output-period>
+   <format name="VtkHdfV2PostProcessor" />
+   <output>
+     <variable>U</variable>
+   </output>
+  </arcane-post-processing>
+
+  <meshes>
+    <mesh>
+      <filename>sphere_cut.msh</filename>
+    </mesh>
+  </meshes>
+
+  <fem>
+    <solution-comparison-file>poisson_test_ref_sphere_3D.txt</solution-comparison-file>
+    <f>5.5</f>
+    <enforce-Dirichlet-method>Penalty</enforce-Dirichlet-method>
+    <penalty>1.e31</penalty>
+    <dirichlet-boundary-condition>
+      <surface>horizontal</surface>
+      <value>0.5</value>
+    </dirichlet-boundary-condition>
+    <bsr>true</bsr>
+    <linear-system name="PETScLinearSystem">
+      <rtol>0.</rtol>
+      <atol>1e-5</atol>
+    </linear-system>
+  </fem>
+</case>