Update for backsubstitution bug fix.

Author: chrispbradley

src/c/expected_results/LaplaceEquation.part0.exnode

@@ -13,94 +13,94 @@
   0.0000000000000000E+00
   0.0000000000000000E+00
   0.0000000000000000E+00
-  -3.2707901429383768E-01
+  -2.6809351305219919E-01
  Node:            2
   2.0000000000000000E+00
   0.0000000000000000E+00
   0.0000000000000000E+00
   5.1441697033516798E-01
-  1.0143463773488193E-01
+  -2.2204460492503131E-16
  Node:            3
   0.0000000000000000E+00
   3.3333333333333331E-01
   0.0000000000000000E+00
   1.7891898861256758E-01
-  7.1059119943811921E-01
+  -2.2204460492503131E-16
  Node:            4
   2.0000000000000000E+00
   3.3333333333333331E-01
   0.0000000000000000E+00
   4.2342691170625335E-01
-  6.2093061073913491E-01
+  -3.3306690738754696E-16
  Node:            5
   0.0000000000000000E+00
   6.6666666666666663E-01
   0.0000000000000000E+00
   2.7840556692655594E-01
-  7.1059119943811921E-01
+  6.6613381477509392E-16
  Node:            6
   2.0000000000000000E+00
   6.6666666666666663E-01
   0.0000000000000000E+00
   3.5372850150917595E-01
-  6.2093061073913502E-01
+  -1.1102230246251565E-16
  Node:            7
   0.0000000000000000E+00
   1.0000000000000000E+00
   0.0000000000000000E+00
   3.4274305408244493E-01
-  2.5370520868685476E-01
+  -3.8857805861880479E-16
  Node:            8
   2.0000000000000000E+00
   1.0000000000000000E+00
   0.0000000000000000E+00
   2.6103842060345522E-01
-  -3.4680425042767582E-02
+  -2.2204460492503131E-16
  Node:            9
   0.0000000000000000E+00
   0.0000000000000000E+00
   3.0000000000000000E+00
   7.3896157939653973E-01
-  -1.8309866618667348E-01
+  1.1102230246251565E-16
  Node:            10
   2.0000000000000000E+00
   0.0000000000000000E+00
   3.0000000000000000E+00
   6.5725694591755268E-01
-  1.4064952969343414E-01
+  0.0000000000000000E+00
  Node:            11
   0.0000000000000000E+00
   3.3333333333333331E-01
   3.0000000000000000E+00
   6.4627149849081977E-01
-  1.4088754611866103E+00
+  -1.1102230246251565E-16
  Node:            12
   2.0000000000000000E+00
   3.3333333333333331E-01
   3.0000000000000000E+00
   7.2159443307344184E-01
-  1.2147895567375782E+00
+  8.8817841970012523E-16
  Node:            13
   0.0000000000000000E+00
   6.6666666666666663E-01
   3.0000000000000000E+00
   5.7657308829374232E-01
-  1.4088754611866103E+00
+  0.0000000000000000E+00
  Node:            14
   2.0000000000000000E+00
   6.6666666666666663E-01
   3.0000000000000000E+00
   8.2108101138743061E-01
-  1.2147895567375782E+00
+  -4.4408920985006262E-16
  Node:            15
   0.0000000000000000E+00
   1.0000000000000000E+00
   3.0000000000000000E+00
   4.8558302966482708E-01
-  1.6240457247709605E-01
+  0.0000000000000000E+00
  Node:            16
   2.0000000000000000E+00
   1.0000000000000000E+00
   3.0000000000000000E+00
   1.0000000000000000E+00
-  -1.1333584306898714E-01
+  2.6809351305220197E-01

src/fortran/expected_results/LaplaceEquation.part0.exnode

@@ -13,94 +13,94 @@
   0.0000000000000000E+00
   0.0000000000000000E+00
   0.0000000000000000E+00
-  -3.2707901429383768E-01
+  -2.6809351305219997E-01
  Node:            2
   2.0000000000000000E+00
   0.0000000000000000E+00
   0.0000000000000000E+00
-  5.1441697033516798E-01
-  1.0143463773488193E-01
+  5.1441697033516887E-01
+  8.8817841970012523E-16
  Node:            3
   0.0000000000000000E+00
   3.3333333333333331E-01
   0.0000000000000000E+00
-  1.7891898861256758E-01
-  7.1059119943811921E-01
+  1.7891898861256827E-01
+  3.2196467714129540E-15
  Node:            4
   2.0000000000000000E+00
   3.3333333333333331E-01
   0.0000000000000000E+00
   4.2342691170625335E-01
-  6.2093061073913491E-01
+  -1.1102230246251565E-15
  Node:            5
   0.0000000000000000E+00
   6.6666666666666663E-01
   0.0000000000000000E+00
-  2.7840556692655594E-01
-  7.1059119943811921E-01
+  2.7840556692655516E-01
+  -1.5543122344752192E-15
  Node:            6
   2.0000000000000000E+00
   6.6666666666666663E-01
   0.0000000000000000E+00
-  3.5372850150917595E-01
-  6.2093061073913502E-01
+  3.5372850150917712E-01
+  1.9984014443252818E-15
  Node:            7
   0.0000000000000000E+00
   1.0000000000000000E+00
   0.0000000000000000E+00
-  3.4274305408244493E-01
-  2.5370520868685476E-01
+  3.4274305408244438E-01
+  -9.4368957093138306E-16
  Node:            8
   2.0000000000000000E+00
   1.0000000000000000E+00
   0.0000000000000000E+00
-  2.6103842060345522E-01
-  -3.4680425042767582E-02
+  2.6103842060345495E-01
+  -2.2204460492503131E-15
  Node:            9
   0.0000000000000000E+00
   0.0000000000000000E+00
   3.0000000000000000E+00
-  7.3896157939653973E-01
-  -1.8309866618667348E-01
+  7.3896157939654039E-01
+  -3.3306690738754696E-16
  Node:            10
   2.0000000000000000E+00
   0.0000000000000000E+00
   3.0000000000000000E+00
-  6.5725694591755268E-01
-  1.4064952969343414E-01
+  6.5725694591755346E-01
+  4.4408920985006262E-16
  Node:            11
   0.0000000000000000E+00
   3.3333333333333331E-01
   3.0000000000000000E+00
-  6.4627149849081977E-01
-  1.4088754611866103E+00
+  6.4627149849082066E-01
+  2.9976021664879227E-15
  Node:            12
   2.0000000000000000E+00
   3.3333333333333331E-01
   3.0000000000000000E+00
-  7.2159443307344184E-01
-  1.2147895567375782E+00
+  7.2159443307344251E-01
+  1.7763568394002505E-15
  Node:            13
   0.0000000000000000E+00
   6.6666666666666663E-01
   3.0000000000000000E+00
-  5.7657308829374232E-01
-  1.4088754611866103E+00
+  5.7657308829374310E-01
+  -2.6645352591003757E-15
  Node:            14
   2.0000000000000000E+00
   6.6666666666666663E-01
   3.0000000000000000E+00
-  8.2108101138743061E-01
-  1.2147895567375782E+00
+  8.2108101138743017E-01
+  -1.9984014443252818E-15
  Node:            15
   0.0000000000000000E+00
   1.0000000000000000E+00
   3.0000000000000000E+00
-  4.8558302966482708E-01
-  1.6240457247709605E-01
+  4.8558302966482819E-01
+  8.8817841970012523E-16
  Node:            16
   2.0000000000000000E+00
   1.0000000000000000E+00
   3.0000000000000000E+00
   1.0000000000000000E+00
-  -1.1333584306898714E-01
+  2.6809351305220153E-01

src/fortran/laplace_equation.F90

@@ -215,14 +215,14 @@ PROGRAM LaplaceEquation
   !Create a decomposition
   CALL cmfe_Decomposition_Initialise(decomposition,err)
   CALL cmfe_Decomposition_CreateStart(DECOMPOSITION_USER_NUMBER,mesh,decomposition,err)
-! Set the decomposition to be a general decomposition with the specified number of domains
+  !Set the decomposition to be a general decomposition with the specified number of domains
   CALL cmfe_Decomposition_TypeSet(decomposition,CMFE_DECOMPOSITION_CALCULATED_TYPE,err)
   CALL cmfe_Decomposition_NumberOfDomainsSet(decomposition,numberOfComputationalNodes,err)
   !Finish the decomposition
   CALL cmfe_Decomposition_CreateFinish(decomposition,err)
 
   !Destory the mesh now that we have decomposed it
-! CALL cmfe_Mesh_Destroy(mesh,err)
+  !CALL cmfe_Mesh_Destroy(mesh,err)
 
   !-----------------------------------------------------------------------------------------------------------
   ! GEOMETRIC FIELD
@@ -284,12 +284,12 @@ PROGRAM LaplaceEquation
   CALL cmfe_EquationsSet_EquationsCreateStart(equationsSet,equations,err)
   !Set the equations matrices sparsity type
   CALL cmfe_Equations_SparsityTypeSet(equations,CMFE_EQUATIONS_SPARSE_MATRICES,err)
-! CALL cmfe_Equations_SparsityTypeSet(equations,CMFE_EQUATIONS_FULL_MATRICES,err)
+  !CALL cmfe_Equations_SparsityTypeSet(equations,CMFE_EQUATIONS_FULL_MATRICES,err)
   !Set the equations set output
-  CALL cmfe_Equations_OutputTypeSet(equations,CMFE_EQUATIONS_NO_OUTPUT,err)
-! CALL cmfe_Equations_OutputTypeSet(equations,CMFE_EQUATIONS_TIMING_OUTPUT,err)
-! CALL cmfe_Equations_OutputTypeSet(equations,CMFE_EQUATIONS_MATRIX_OUTPUT,err)
-! CALL cmfe_Equations_OutputTypeSet(equations,CMFE_EQUATIONS_ELEMENT_MATRIX_OUTPUT,err)
+  !CALL cmfe_Equations_OutputTypeSet(equations,CMFE_EQUATIONS_NO_OUTPUT,err)
+  !CALL cmfe_Equations_OutputTypeSet(equations,CMFE_EQUATIONS_TIMING_OUTPUT,err)
+  !CALL cmfe_Equations_OutputTypeSet(equations,CMFE_EQUATIONS_MATRIX_OUTPUT,err)
+  CALL cmfe_Equations_OutputTypeSet(equations,CMFE_EQUATIONS_ELEMENT_MATRIX_OUTPUT,err)
   !Finish the equations set equations
   CALL cmfe_EquationsSet_EquationsCreateFinish(equationsSet,err)
 
@@ -317,23 +317,22 @@ PROGRAM LaplaceEquation
   CALL cmfe_Solver_Initialise(solver,err)
   CALL cmfe_Problem_SolversCreateStart(problem,err)
   CALL cmfe_Problem_SolverGet(problem,CMFE_CONTROL_LOOP_NODE,1,solver,err)
-  CALL cmfe_Solver_OutputTypeSet(solver,CMFE_SOLVER_NO_OUTPUT,err)
-! CALL cmfe_Solver_OutputTypeSet(solver,CMFE_SOLVER_PROGRESS_OUTPUT,err)
-! CALL cmfe_Solver_OutputTypeSet(solver,CMFE_SOLVER_TIMING_OUTPUT,err)
-! CALL cmfe_Solver_OutputTypeSet(solver,CMFE_SOLVER_SOLVER_OUTPUT,err)
-! CALL cmfe_Solver_OutputTypeSet(solver,CMFE_SOLVER_MATRIX_OUTPUT,err)
-
-! CALL cmfe_Solver_LinearTypeSet(solver,CMFE_SOLVER_LINEAR_ITERATIVE_SOLVE_TYPE,err)
-! CALL cmfe_Solver_LinearIterativeAbsoluteToleranceSet(solver,1.0E-12_CMISSRP,err)
-! CALL cmfe_Solver_LinearIterativeRelativeToleranceSet(solver,1.0E-12_CMISSRP,err)
-
-  CALL cmfe_Solver_LinearTypeSet(solver,CMFE_SOLVER_LINEAR_DIRECT_SOLVE_TYPE,err)
-
-! CALL cmfe_Solver_LinearTypeSet(solver,CMFE_SOLVER_LINEAR_DIRECT_SOLVE_TYPE,err)
-! CALL cmfe_Solver_LibraryTypeSet(solver,CMFE_SOLVER_MUMPS_LIBRARY,err)
-! CALL cmfe_Solver_LibraryTypeSet(solver,CMFE_SOLVER_LAPACK_LIBRARY,err)
-! CALL cmfe_Solver_LibraryTypeSet(solver,CMFE_SOLVER_SUPERLU_LIBRARY,err)
-! CALL cmfe_Solver_LibraryTypeSet(solver,CMFE_SOLVER_PASTIX_LIBRARY,err)
+  !CALL cmfe_Solver_OutputTypeSet(solver,CMFE_SOLVER_NO_OUTPUT,err)
+  !CALL cmfe_Solver_OutputTypeSet(solver,CMFE_SOLVER_PROGRESS_OUTPUT,err)
+  !CALL cmfe_Solver_OutputTypeSet(solver,CMFE_SOLVER_TIMING_OUTPUT,err)
+  !CALL cmfe_Solver_OutputTypeSet(solver,CMFE_SOLVER_SOLVER_OUTPUT,err)
+  CALL cmfe_Solver_OutputTypeSet(solver,CMFE_SOLVER_MATRIX_OUTPUT,err)
+  
+  CALL cmfe_Solver_LinearTypeSet(solver,CMFE_SOLVER_LINEAR_ITERATIVE_SOLVE_TYPE,err)
+  CALL cmfe_Solver_LinearIterativeAbsoluteToleranceSet(solver,1.0E-12_CMISSRP,err)
+  CALL cmfe_Solver_LinearIterativeRelativeToleranceSet(solver,1.0E-12_CMISSRP,err)
+  
+  !CALL cmfe_Solver_LinearTypeSet(solver,CMFE_SOLVER_LINEAR_DIRECT_SOLVE_TYPE,err)
+  !CALL cmfe_Solver_LinearTypeSet(solver,CMFE_SOLVER_LINEAR_DIRECT_SOLVE_TYPE,err)
+  !CALL cmfe_Solver_LibraryTypeSet(solver,CMFE_SOLVER_MUMPS_LIBRARY,err)
+  !CALL cmfe_Solver_LibraryTypeSet(solver,CMFE_SOLVER_LAPACK_LIBRARY,err)
+  !CALL cmfe_Solver_LibraryTypeSet(solver,CMFE_SOLVER_SUPERLU_LIBRARY,err)
+  !CALL cmfe_Solver_LibraryTypeSet(solver,CMFE_SOLVER_PASTIX_LIBRARY,err)
   !Finish the creation of the problem solver
   CALL cmfe_Problem_SolversCreateFinish(problem,err)