diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/circle_1_axi.gml b/Tests/lfs-data/Elliptic/circle_radius_1/circle_1_axi.gml
new file mode 100644
index 0000000000000000000000000000000000000000..3d3f95d480476fdd81f7b0e8d6c7b3c355f12aeb
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/circle_1_axi.gml
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:80301740d9f9a3428b0801f2439a481f430a94298edd1f40866fb969cd473825
+size 480
diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/circle_1e1_axi.prj b/Tests/lfs-data/Elliptic/circle_radius_1/circle_1e1_axi.prj
new file mode 100644
index 0000000000000000000000000000000000000000..d34028d67a67a0a45760dac0256bdb9180bcebf1
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/circle_1e1_axi.prj
@@ -0,0 +1,112 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+    <mesh axially_symmetric="true">line_1_lines_1e1.vtu</mesh>
+    <geometry>circle_1_axi.gml</geometry>
+    <processes>
+        <process>
+            <name>GW23</name>
+            <type>GROUNDWATER_FLOW</type>
+            <integration_order>2</integration_order>
+            <hydraulic_conductivity>kappa</hydraulic_conductivity>
+            <process_variables>
+                <process_variable>pressure</process_variable>
+            </process_variables>
+            <secondary_variables>
+                <secondary_variable type="static" internal_name="darcy_velocity" output_name="v"/>
+            </secondary_variables>
+        </process>
+    </processes>
+    <time_loop>
+        <processes>
+            <process ref="GW23">
+                <nonlinear_solver>basic_picard</nonlinear_solver>
+                <convergence_criterion>
+                    <type>DeltaX</type>
+                    <norm_type>NORM2</norm_type>
+                    <abstol>1.e-6</abstol>
+                </convergence_criterion>
+                <time_discretization>
+                    <type>BackwardEuler</type>
+                </time_discretization>
+                <output>
+                    <variables>
+                        <variable> pressure </variable>
+                        <variable> v      </variable>
+                    </variables>
+                </output>
+                <time_stepping>
+                    <type>SingleStep</type>
+                </time_stepping>
+            </process>
+        </processes>
+        <output>
+            <type>VTK</type>
+            <prefix>circle_1e1_axi</prefix>
+        </output>
+    </time_loop>
+    <nonlinear_solvers>
+        <nonlinear_solver>
+            <name>basic_picard</name>
+            <type>Picard</type>
+            <max_iter>10</max_iter>
+            <linear_solver>general_linear_solver</linear_solver>
+        </nonlinear_solver>
+    </nonlinear_solvers>
+    <linear_solvers>
+        <linear_solver>
+            <name>general_linear_solver</name>
+            <lis>-i cg -p jacobi -tol 1e-16 -maxiter 10000</lis>
+            <eigen>
+                <solver_type>CG</solver_type>
+                <precon_type>DIAGONAL</precon_type>
+                <max_iteration_step>10000</max_iteration_step>
+                <error_tolerance>1e-16</error_tolerance>
+            </eigen>
+            <petsc>
+                <prefix>gw</prefix>
+                <parameters>-gw_ksp_type cg -gw_pc_type bjacobi -gw_ksp_rtol 1e-16 -gw_ksp_max_it 10000</parameters>
+            </petsc>
+        </linear_solver>
+    </linear_solvers>
+    <parameters>
+        <parameter>
+            <name>kappa</name>
+            <type>Constant</type>
+            <value>1</value>
+        </parameter>
+        <parameter>
+            <name>p0</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+        <parameter>
+            <name>pressure_bc</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+    </parameters>
+    <process_variables>
+        <process_variable>
+            <name>pressure</name>
+            <components>1</components>
+            <order>1</order>
+            <initial_condition>p0</initial_condition>
+            <boundary_conditions>
+                <boundary_condition>
+                    <geometrical_set>geometry</geometrical_set>
+                    <geometry>outer</geometry>
+                    <type>Dirichlet</type>
+                    <parameter>pressure_bc</parameter>
+                </boundary_condition>
+            </boundary_conditions>
+            <source_terms>
+                <source_term>
+                    <geometrical_set>geometry</geometrical_set>
+                    <geometry>inner</geometry>
+                    <type>Nodal</type>
+                    <value>1</value>
+                </source_term>
+            </source_terms>
+        </process_variable>
+    </process_variables>
+</OpenGeoSysProject>
diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/circle_1e2_axi.prj b/Tests/lfs-data/Elliptic/circle_radius_1/circle_1e2_axi.prj
new file mode 100644
index 0000000000000000000000000000000000000000..77b9d2c4e7f97f408695b52e66da16dcc663e5d4
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/circle_1e2_axi.prj
@@ -0,0 +1,112 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+    <mesh axially_symmetric="true">line_1_lines_1e2.vtu</mesh>
+    <geometry>circle_1_axi.gml</geometry>
+    <processes>
+        <process>
+            <name>GW23</name>
+            <type>GROUNDWATER_FLOW</type>
+            <integration_order>2</integration_order>
+            <hydraulic_conductivity>kappa</hydraulic_conductivity>
+            <process_variables>
+                <process_variable>pressure</process_variable>
+            </process_variables>
+            <secondary_variables>
+                <secondary_variable type="static" internal_name="darcy_velocity" output_name="v"/>
+            </secondary_variables>
+        </process>
+    </processes>
+    <time_loop>
+        <processes>
+            <process ref="GW23">
+                <nonlinear_solver>basic_picard</nonlinear_solver>
+                <convergence_criterion>
+                    <type>DeltaX</type>
+                    <norm_type>NORM2</norm_type>
+                    <abstol>1.e-6</abstol>
+                </convergence_criterion>
+                <time_discretization>
+                    <type>BackwardEuler</type>
+                </time_discretization>
+                <output>
+                    <variables>
+                        <variable> pressure </variable>
+                        <variable> v      </variable>
+                    </variables>
+                </output>
+                <time_stepping>
+                    <type>SingleStep</type>
+                </time_stepping>
+            </process>
+        </processes>
+        <output>
+            <type>VTK</type>
+            <prefix>circle_1e2_axi</prefix>
+        </output>
+    </time_loop>
+    <nonlinear_solvers>
+        <nonlinear_solver>
+            <name>basic_picard</name>
+            <type>Picard</type>
+            <max_iter>10</max_iter>
+            <linear_solver>general_linear_solver</linear_solver>
+        </nonlinear_solver>
+    </nonlinear_solvers>
+    <linear_solvers>
+        <linear_solver>
+            <name>general_linear_solver</name>
+            <lis>-i cg -p jacobi -tol 1e-16 -maxiter 10000</lis>
+            <eigen>
+                <solver_type>CG</solver_type>
+                <precon_type>DIAGONAL</precon_type>
+                <max_iteration_step>10000</max_iteration_step>
+                <error_tolerance>1e-16</error_tolerance>
+            </eigen>
+            <petsc>
+                <prefix>gw</prefix>
+                <parameters>-gw_ksp_type cg -gw_pc_type bjacobi -gw_ksp_rtol 1e-16 -gw_ksp_max_it 10000</parameters>
+            </petsc>
+        </linear_solver>
+    </linear_solvers>
+    <parameters>
+        <parameter>
+            <name>kappa</name>
+            <type>Constant</type>
+            <value>1</value>
+        </parameter>
+        <parameter>
+            <name>p0</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+        <parameter>
+            <name>pressure_bc</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+    </parameters>
+    <process_variables>
+        <process_variable>
+            <name>pressure</name>
+            <components>1</components>
+            <order>1</order>
+            <initial_condition>p0</initial_condition>
+            <boundary_conditions>
+                <boundary_condition>
+                    <geometrical_set>geometry</geometrical_set>
+                    <geometry>outer</geometry>
+                    <type>Dirichlet</type>
+                    <parameter>pressure_bc</parameter>
+                </boundary_condition>
+            </boundary_conditions>
+            <source_terms>
+                <source_term>
+                    <geometrical_set>geometry</geometrical_set>
+                    <geometry>inner</geometry>
+                    <type>Nodal</type>
+                    <value>1</value>
+                </source_term>
+            </source_terms>
+        </process_variable>
+    </process_variables>
+</OpenGeoSysProject>
diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/circle_1e3_axi.prj b/Tests/lfs-data/Elliptic/circle_radius_1/circle_1e3_axi.prj
new file mode 100644
index 0000000000000000000000000000000000000000..89c7ddc6b473046ea101b5e55758948b9194e16e
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/circle_1e3_axi.prj
@@ -0,0 +1,112 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+    <mesh axially_symmetric="true">line_1_lines_1e3.vtu</mesh>
+    <geometry>circle_1_axi.gml</geometry>
+    <processes>
+        <process>
+            <name>GW23</name>
+            <type>GROUNDWATER_FLOW</type>
+            <integration_order>2</integration_order>
+            <hydraulic_conductivity>kappa</hydraulic_conductivity>
+            <process_variables>
+                <process_variable>pressure</process_variable>
+            </process_variables>
+            <secondary_variables>
+                <secondary_variable type="static" internal_name="darcy_velocity" output_name="v"/>
+            </secondary_variables>
+        </process>
+    </processes>
+    <time_loop>
+        <processes>
+            <process ref="GW23">
+                <nonlinear_solver>basic_picard</nonlinear_solver>
+                <convergence_criterion>
+                    <type>DeltaX</type>
+                    <norm_type>NORM2</norm_type>
+                    <abstol>1.e-6</abstol>
+                </convergence_criterion>
+                <time_discretization>
+                    <type>BackwardEuler</type>
+                </time_discretization>
+                <output>
+                    <variables>
+                        <variable> pressure </variable>
+                        <variable> v      </variable>
+                    </variables>
+                </output>
+                <time_stepping>
+                    <type>SingleStep</type>
+                </time_stepping>
+            </process>
+        </processes>
+        <output>
+            <type>VTK</type>
+            <prefix>circle_1e3_axi</prefix>
+        </output>
+    </time_loop>
+    <nonlinear_solvers>
+        <nonlinear_solver>
+            <name>basic_picard</name>
+            <type>Picard</type>
+            <max_iter>10</max_iter>
+            <linear_solver>general_linear_solver</linear_solver>
+        </nonlinear_solver>
+    </nonlinear_solvers>
+    <linear_solvers>
+        <linear_solver>
+            <name>general_linear_solver</name>
+            <lis>-i cg -p jacobi -tol 1e-16 -maxiter 10000</lis>
+            <eigen>
+                <solver_type>CG</solver_type>
+                <precon_type>DIAGONAL</precon_type>
+                <max_iteration_step>50000</max_iteration_step>
+                <error_tolerance>1e-13</error_tolerance>
+            </eigen>
+            <petsc>
+                <prefix>gw</prefix>
+                <parameters>-gw_ksp_type cg -gw_pc_type bjacobi -gw_ksp_rtol 1e-16 -gw_ksp_max_it 10000</parameters>
+            </petsc>
+        </linear_solver>
+    </linear_solvers>
+    <parameters>
+        <parameter>
+            <name>kappa</name>
+            <type>Constant</type>
+            <value>1</value>
+        </parameter>
+        <parameter>
+            <name>p0</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+        <parameter>
+            <name>pressure_bc</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+    </parameters>
+    <process_variables>
+        <process_variable>
+            <name>pressure</name>
+            <components>1</components>
+            <order>1</order>
+            <initial_condition>p0</initial_condition>
+            <boundary_conditions>
+                <boundary_condition>
+                    <geometrical_set>geometry</geometrical_set>
+                    <geometry>outer</geometry>
+                    <type>Dirichlet</type>
+                    <parameter>pressure_bc</parameter>
+                </boundary_condition>
+            </boundary_conditions>
+            <source_terms>
+                <source_term>
+                    <geometrical_set>geometry</geometrical_set>
+                    <geometry>inner</geometry>
+                    <type>Nodal</type>
+                    <value>1</value>
+                </source_term>
+            </source_terms>
+        </process_variable>
+    </process_variables>
+</OpenGeoSysProject>
diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/circle_1e4_axi.prj b/Tests/lfs-data/Elliptic/circle_radius_1/circle_1e4_axi.prj
new file mode 100644
index 0000000000000000000000000000000000000000..2ef9e67c9d9a6ba8876d82abd519f694ed68e942
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/circle_1e4_axi.prj
@@ -0,0 +1,112 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+    <mesh axially_symmetric="true">line_1_lines_1e4.vtu</mesh>
+    <geometry>circle_1_axi.gml</geometry>
+    <processes>
+        <process>
+            <name>GW23</name>
+            <type>GROUNDWATER_FLOW</type>
+            <integration_order>2</integration_order>
+            <hydraulic_conductivity>kappa</hydraulic_conductivity>
+            <process_variables>
+                <process_variable>pressure</process_variable>
+            </process_variables>
+            <secondary_variables>
+                <secondary_variable type="static" internal_name="darcy_velocity" output_name="v"/>
+            </secondary_variables>
+        </process>
+    </processes>
+    <time_loop>
+        <processes>
+            <process ref="GW23">
+                <nonlinear_solver>basic_picard</nonlinear_solver>
+                <convergence_criterion>
+                    <type>DeltaX</type>
+                    <norm_type>NORM2</norm_type>
+                    <abstol>1.e-6</abstol>
+                </convergence_criterion>
+                <time_discretization>
+                    <type>BackwardEuler</type>
+                </time_discretization>
+                <output>
+                    <variables>
+                        <variable> pressure </variable>
+                        <variable> v      </variable>
+                    </variables>
+                </output>
+                <time_stepping>
+                    <type>SingleStep</type>
+                </time_stepping>
+            </process>
+        </processes>
+        <output>
+            <type>VTK</type>
+            <prefix>circle_1e4_axi</prefix>
+        </output>
+    </time_loop>
+    <nonlinear_solvers>
+        <nonlinear_solver>
+            <name>basic_picard</name>
+            <type>Picard</type>
+            <max_iter>10</max_iter>
+            <linear_solver>general_linear_solver</linear_solver>
+        </nonlinear_solver>
+    </nonlinear_solvers>
+    <linear_solvers>
+        <linear_solver>
+            <name>general_linear_solver</name>
+            <lis>-i cg -p jacobi -tol 1e-16 -maxiter 10000</lis>
+            <eigen>
+                <solver_type>CG</solver_type>
+                <precon_type>DIAGONAL</precon_type>
+                <max_iteration_step>20000</max_iteration_step>
+                <error_tolerance>1e-14</error_tolerance>
+            </eigen>
+            <petsc>
+                <prefix>gw</prefix>
+                <parameters>-gw_ksp_type cg -gw_pc_type bjacobi -gw_ksp_rtol 1e-16 -gw_ksp_max_it 10000</parameters>
+            </petsc>
+        </linear_solver>
+    </linear_solvers>
+    <parameters>
+        <parameter>
+            <name>kappa</name>
+            <type>Constant</type>
+            <value>1</value>
+        </parameter>
+        <parameter>
+            <name>p0</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+        <parameter>
+            <name>pressure_bc</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+    </parameters>
+    <process_variables>
+        <process_variable>
+            <name>pressure</name>
+            <components>1</components>
+            <order>1</order>
+            <initial_condition>p0</initial_condition>
+            <boundary_conditions>
+                <boundary_condition>
+                    <geometrical_set>geometry</geometrical_set>
+                    <geometry>outer</geometry>
+                    <type>Dirichlet</type>
+                    <parameter>pressure_bc</parameter>
+                </boundary_condition>
+            </boundary_conditions>
+            <source_terms>
+                <source_term>
+                    <geometrical_set>geometry</geometrical_set>
+                    <geometry>inner</geometry>
+                    <type>Nodal</type>
+                    <value>1</value>
+                </source_term>
+            </source_terms>
+        </process_variable>
+    </process_variables>
+</OpenGeoSysProject>
diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/circle_1e5_axi.prj b/Tests/lfs-data/Elliptic/circle_radius_1/circle_1e5_axi.prj
new file mode 100644
index 0000000000000000000000000000000000000000..436391d00c192f63916919435c04998faca7326f
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/circle_1e5_axi.prj
@@ -0,0 +1,112 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+    <mesh axially_symmetric="true">line_1_lines_1e5.vtu</mesh>
+    <geometry>circle_1_axi.gml</geometry>
+    <processes>
+        <process>
+            <name>GW23</name>
+            <type>GROUNDWATER_FLOW</type>
+            <integration_order>2</integration_order>
+            <hydraulic_conductivity>kappa</hydraulic_conductivity>
+            <process_variables>
+                <process_variable>pressure</process_variable>
+            </process_variables>
+            <secondary_variables>
+                <secondary_variable type="static" internal_name="darcy_velocity" output_name="v"/>
+            </secondary_variables>
+        </process>
+    </processes>
+    <time_loop>
+        <processes>
+            <process ref="GW23">
+                <nonlinear_solver>basic_picard</nonlinear_solver>
+                <convergence_criterion>
+                    <type>DeltaX</type>
+                    <norm_type>NORM2</norm_type>
+                    <abstol>1.e-6</abstol>
+                </convergence_criterion>
+                <time_discretization>
+                    <type>BackwardEuler</type>
+                </time_discretization>
+                <output>
+                    <variables>
+                        <variable> pressure </variable>
+                        <variable> v      </variable>
+                    </variables>
+                </output>
+                <time_stepping>
+                    <type>SingleStep</type>
+                </time_stepping>
+            </process>
+        </processes>
+        <output>
+            <type>VTK</type>
+            <prefix>circle_1e5_axi</prefix>
+        </output>
+    </time_loop>
+    <nonlinear_solvers>
+        <nonlinear_solver>
+            <name>basic_picard</name>
+            <type>Picard</type>
+            <max_iter>10</max_iter>
+            <linear_solver>general_linear_solver</linear_solver>
+        </nonlinear_solver>
+    </nonlinear_solvers>
+    <linear_solvers>
+        <linear_solver>
+            <name>general_linear_solver</name>
+            <lis>-i cg -p jacobi -tol 1e-16 -maxiter 10000</lis>
+            <eigen>
+                <solver_type>SparseLU</solver_type>
+                <!--precon_type>DIAGONAL</precon_type>
+                <max_iteration_step>60000</max_iteration_step>
+                <error_tolerance>1e-6</error_tolerance-->
+            </eigen>
+            <petsc>
+                <prefix>gw</prefix>
+                <parameters>-gw_ksp_type cg -gw_pc_type bjacobi -gw_ksp_rtol 1e-16 -gw_ksp_max_it 10000</parameters>
+            </petsc>
+        </linear_solver>
+    </linear_solvers>
+    <parameters>
+        <parameter>
+            <name>kappa</name>
+            <type>Constant</type>
+            <value>1</value>
+        </parameter>
+        <parameter>
+            <name>p0</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+        <parameter>
+            <name>pressure_bc</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+    </parameters>
+    <process_variables>
+        <process_variable>
+            <name>pressure</name>
+            <components>1</components>
+            <order>1</order>
+            <initial_condition>p0</initial_condition>
+            <boundary_conditions>
+                <boundary_condition>
+                    <geometrical_set>geometry</geometrical_set>
+                    <geometry>outer</geometry>
+                    <type>Dirichlet</type>
+                    <parameter>pressure_bc</parameter>
+                </boundary_condition>
+            </boundary_conditions>
+            <source_terms>
+                <source_term>
+                    <geometrical_set>geometry</geometrical_set>
+                    <geometry>inner</geometry>
+                    <type>Nodal</type>
+                    <value>1</value>
+                </source_term>
+            </source_terms>
+        </process_variable>
+    </process_variables>
+</OpenGeoSysProject>
diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/circle_1e6_axi.prj b/Tests/lfs-data/Elliptic/circle_radius_1/circle_1e6_axi.prj
new file mode 100644
index 0000000000000000000000000000000000000000..86470dd0ef918fdfcf28e6ccf9ee4c75eda7afb9
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/circle_1e6_axi.prj
@@ -0,0 +1,109 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+    <mesh axially_symmetric="true">line_1_lines_1e6.vtu</mesh>
+    <geometry>circle_1_axi.gml</geometry>
+    <processes>
+        <process>
+            <name>GW23</name>
+            <type>GROUNDWATER_FLOW</type>
+            <integration_order>2</integration_order>
+            <hydraulic_conductivity>kappa</hydraulic_conductivity>
+            <process_variables>
+                <process_variable>pressure</process_variable>
+            </process_variables>
+            <secondary_variables>
+                <secondary_variable type="static" internal_name="darcy_velocity" output_name="v"/>
+            </secondary_variables>
+        </process>
+    </processes>
+    <time_loop>
+        <processes>
+            <process ref="GW23">
+                <nonlinear_solver>basic_picard</nonlinear_solver>
+                <convergence_criterion>
+                    <type>DeltaX</type>
+                    <norm_type>NORM2</norm_type>
+                    <abstol>1.e-6</abstol>
+                </convergence_criterion>
+                <time_discretization>
+                    <type>BackwardEuler</type>
+                </time_discretization>
+                <output>
+                    <variables>
+                        <variable> pressure </variable>
+                        <variable> v      </variable>
+                    </variables>
+                </output>
+                <time_stepping>
+                    <type>SingleStep</type>
+                </time_stepping>
+            </process>
+        </processes>
+        <output>
+            <type>VTK</type>
+            <prefix>circle_1e6_axi</prefix>
+        </output>
+    </time_loop>
+    <nonlinear_solvers>
+        <nonlinear_solver>
+            <name>basic_picard</name>
+            <type>Picard</type>
+            <max_iter>10</max_iter>
+            <linear_solver>general_linear_solver</linear_solver>
+        </nonlinear_solver>
+    </nonlinear_solvers>
+    <linear_solvers>
+        <linear_solver>
+            <name>general_linear_solver</name>
+            <lis>-i cg -p jacobi -tol 1e-16 -maxiter 10000</lis>
+            <eigen>
+                <solver_type>SparseLU</solver_type>
+            </eigen>
+            <petsc>
+                <prefix>gw</prefix>
+                <parameters>-gw_ksp_type cg -gw_pc_type bjacobi -gw_ksp_rtol 1e-16 -gw_ksp_max_it 10000</parameters>
+            </petsc>
+        </linear_solver>
+    </linear_solvers>
+    <parameters>
+        <parameter>
+            <name>kappa</name>
+            <type>Constant</type>
+            <value>1</value>
+        </parameter>
+        <parameter>
+            <name>p0</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+        <parameter>
+            <name>pressure_bc</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+    </parameters>
+    <process_variables>
+        <process_variable>
+            <name>pressure</name>
+            <components>1</components>
+            <order>1</order>
+            <initial_condition>p0</initial_condition>
+            <boundary_conditions>
+                <boundary_condition>
+                    <geometrical_set>geometry</geometrical_set>
+                    <geometry>outer</geometry>
+                    <type>Dirichlet</type>
+                    <parameter>pressure_bc</parameter>
+                </boundary_condition>
+            </boundary_conditions>
+            <source_terms>
+                <source_term>
+                    <geometrical_set>geometry</geometrical_set>
+                    <geometry>inner</geometry>
+                    <type>Nodal</type>
+                    <value>1</value>
+                </source_term>
+            </source_terms>
+        </process_variable>
+    </process_variables>
+</OpenGeoSysProject>
diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e1.vtu b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e1.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..9f40daa18d3d9b3e6faf1d6ca3538175abac4be2
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e1.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:53d1a2370aa28aeba2a319af9f291478b0aa5b1d028e658151e9255904119607
+size 1463
diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e1_expected.vtu b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e1_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..3600c1c212a2375e2c8ae96e04b4c78fb78a8600
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e1_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:a66c26f401a08aa2ec20457681ec2951b8ddcd8bab90ee23afe87c7f5b33522d
+size 1816
diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e2.vtu b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e2.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..313464d2ea316276f42b50f434da333319770d44
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e2.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:7592da0c9492e0b0330710ab89409ed4ba6b9c0d72e63f0d4bb6e580301719ca
+size 2359
diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e2_expected.vtu b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e2_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..858269185341004ff887ef58f5fd327a8fa9d981
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e2_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:fb099137f82d787c9bab2c471aa84f10ad4cb9d1b1053da82610ae74954057f7
+size 3632
diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e3.vtu b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e3.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..ba5493965eb02c424f46f52b86fc3e3223a58d0b
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e3.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:272081dc5b19adfc70a4bdea8b93c819ee5b7cf7a1000959521654e7e1d33535
+size 9771
diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e3_expected.vtu b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e3_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..3920f8119dcdbfa970694af3f83ecace363441de
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e3_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:e338763ddaa9084ce536972dd8c0e20202384b4be26424281e2efef440328fe8
+size 20068
diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e4.vtu b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e4.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..be6c40c4f1301f94353b9b2f8dc26e5045c061d8
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e4.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:84728aed338359f36053dcc6ab57da5e0caae70655f7752a30a5eea2fc1acc5d
+size 102595
diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e4_expected.vtu b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e4_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..b075e6ac333b72c7b695a10d079b901095b474aa
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e4_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:7fd038ae3e9bdca15f409719aad6847aa669c5807ae4b3501437ca410d079ee2
+size 202241
diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e5.vtu b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e5.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..98ebcdd388bd7faf8129f74cbe4e4b64389f173f
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e5.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:57f1f07c38a4e8de5041b6602e93235c85f6e49ab2f046ddcaec2476d6f90e99
+size 1230851
diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e5_expected.vtu b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e5_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..7051f756ec29e90b30ab6e779edbf2ef59229ece
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e5_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:600f5a48667f9eb6effef3d3e0af3eddc9a6cd5e46607c49db52fe97caa0a4f0
+size 2074522
diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e6.vtu b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e6.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..2b26b2977367107d382721e269e13b8d811d8bc3
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e6.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:be8ddb9571c1b57b0935f61ae9593d678d11d0b46cbc4a65c78ed393af0f49ba
+size 12419971
diff --git a/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e6_expected.vtu b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e6_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..24eef5373a752e6556ba67b872aa1ab15b084d29
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/circle_radius_1/line_1_lines_1e6_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:0548dc4336277d035d3049442faf2972dd98b1cbfb5b364602848ce68eb32930
+size 20336416
diff --git a/Tests/lfs-data/Elliptic/square_1x1_GroundWaterFlow/square_1e6_with_nodal_sources.prj b/Tests/lfs-data/Elliptic/square_1x1_GroundWaterFlow/square_1e6_with_nodal_sources.prj
new file mode 100644
index 0000000000000000000000000000000000000000..fc58c74366310608744fbb6f68bdf6c4977a6052
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/square_1x1_GroundWaterFlow/square_1e6_with_nodal_sources.prj
@@ -0,0 +1,123 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+    <mesh>square_1x1_quad_1e6.vtu</mesh>
+    <geometry>square_1x1.gml</geometry>
+    <processes>
+        <process>
+            <name>GW23</name>
+            <type>GROUNDWATER_FLOW</type>
+            <integration_order>2</integration_order>
+            <hydraulic_conductivity>K</hydraulic_conductivity>
+            <process_variables>
+                <process_variable>pressure</process_variable>
+            </process_variables>
+            <secondary_variables>
+                <secondary_variable type="static" internal_name="darcy_velocity" output_name="v"/>
+            </secondary_variables>
+        </process>
+    </processes>
+    <time_loop>
+        <processes>
+            <process ref="GW23">
+                <nonlinear_solver>basic_picard</nonlinear_solver>
+                <convergence_criterion>
+                    <type>DeltaX</type>
+                    <norm_type>NORM2</norm_type>
+                    <abstol>1.e-6</abstol>
+                </convergence_criterion>
+                <time_discretization>
+                    <type>BackwardEuler</type>
+                </time_discretization>
+                <output>
+                    <variables>
+                        <variable> pressure </variable>
+                        <variable> v      </variable>
+                    </variables>
+                </output>
+                <time_stepping>
+                    <type>SingleStep</type>
+                </time_stepping>
+            </process>
+        </processes>
+        <output>
+            <type>VTK</type>
+            <prefix>square_1e6</prefix>
+        </output>
+    </time_loop>
+    <parameters>
+        <parameter>
+            <name>K</name>
+            <type>Constant</type>
+            <value>1</value>
+        </parameter>
+        <parameter>
+            <name>p0</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+        <parameter>
+            <name>p_Dirichlet_left</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+        <parameter>
+            <name>p_Dirichlet_right</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+    </parameters>
+    <process_variables>
+        <process_variable>
+            <name>pressure</name>
+            <components>1</components>
+            <order>1</order>
+            <initial_condition>p0</initial_condition>
+            <boundary_conditions>
+                <boundary_condition>
+                    <geometrical_set>square_1x1_geometry</geometrical_set>
+                    <geometry>left</geometry>
+                    <type>Dirichlet</type>
+                    <parameter>p_Dirichlet_left</parameter>
+                </boundary_condition>
+                <boundary_condition>
+                    <geometrical_set>square_1x1_geometry</geometrical_set>
+                    <geometry>right</geometry>
+                    <type>Dirichlet</type>
+                    <parameter>p_Dirichlet_right</parameter>
+                </boundary_condition>
+            </boundary_conditions>
+            <source_terms>
+                <source_term>
+                    <geometrical_set>square_1x1_geometry</geometrical_set>
+                    <geometry>middle_point</geometry>
+                    <type>Nodal</type>
+                    <value>1</value>
+                </source_term>
+            </source_terms>
+        </process_variable>
+    </process_variables>
+    <nonlinear_solvers>
+        <nonlinear_solver>
+            <name>basic_picard</name>
+            <type>Picard</type>
+            <max_iter>10</max_iter>
+            <linear_solver>general_linear_solver</linear_solver>
+        </nonlinear_solver>
+    </nonlinear_solvers>
+    <linear_solvers>
+        <linear_solver>
+            <name>general_linear_solver</name>
+            <lis>-i cg -p jacobi -tol 1e-16 -maxiter 10000</lis>
+            <eigen>
+                <solver_type>CG</solver_type>
+                <precon_type>DIAGONAL</precon_type>
+                <max_iteration_step>10000</max_iteration_step>
+                <error_tolerance>1e-16</error_tolerance>
+            </eigen>
+            <petsc>
+                <prefix>gw</prefix>
+                <parameters>-gw_ksp_type cg -gw_pc_type bjacobi -gw_ksp_rtol 1e-16 -gw_ksp_max_it 10000</parameters>
+            </petsc>
+        </linear_solver>
+    </linear_solvers>
+</OpenGeoSysProject>
diff --git a/Tests/lfs-data/Elliptic/square_1x1_GroundWaterFlow/square_1x1.gml b/Tests/lfs-data/Elliptic/square_1x1_GroundWaterFlow/square_1x1.gml
index b64525fffa2590f930ad45d289e6f631596ae279..3e5bb2a80a8f6bbdc4b7776b696d7a6af0d43b51 100644
--- a/Tests/lfs-data/Elliptic/square_1x1_GroundWaterFlow/square_1x1.gml
+++ b/Tests/lfs-data/Elliptic/square_1x1_GroundWaterFlow/square_1x1.gml
@@ -1,3 +1,3 @@
 version https://git-lfs.github.com/spec/v1
-oid sha256:9a6e7b692b55d4f0b5e29b055a4945bc3e57a822ba3df69ee94a23e83ace6172
-size 1017
+oid sha256:a3f52cd5a8058d967f8ead873e1afa323af39266e538bcd2c35af08683a81dad
+size 1083
diff --git a/Tests/lfs-data/Elliptic/square_1x1_GroundWaterFlow/square_1x1_quad_1e6_nodal_sources_expected.vtu b/Tests/lfs-data/Elliptic/square_1x1_GroundWaterFlow/square_1x1_quad_1e6_nodal_sources_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..7f8270e97169edabb3a7d2db3be28e07ee7e2a1a
--- /dev/null
+++ b/Tests/lfs-data/Elliptic/square_1x1_GroundWaterFlow/square_1x1_quad_1e6_nodal_sources_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:e17abe1e21fc25ab2a760a9db709335fe1e5c8b9eb51b005d306151d7b97762a
+size 20360466