WIP: Add test computation for bending isometries using reducedcubichermitetrianglebasis

problems/L-clamped-Plate.py

@@ -44,37 +44,20 @@ parameterSet.interpolationMethod = 'embedded'
 
 
 
-parameterSet.conforming_DiscreteJacobian = 1
+parameterSet.conforming_DiscreteJacobian = 0
 
 #############################################
 #  Solver parameters
 #############################################
 
 # Tolerance of the multigrid solver
-parameterSet.tolerance = 1e-15
+parameterSet.tolerance = 1e-12
 # Maximum number of multigrid iterations
-parameterSet.maxTrustRegionSteps = 1000
-# Initial trust-region radius
-parameterSet.initialTrustRegionRadius = 0.5
-# Number of multigrid iterations per trust-region step
-parameterSet.numIt = 200
-# Number of presmoothing steps
-parameterSet.nu1 = 3
-# Number of postsmoothing steps
-parameterSet.nu2 = 3
-# Number of coarse grid corrections
-parameterSet.mu = 1
-# Number of base solver iterations
-parameterSet.baseIt = 100
-# Tolerance of the multigrid solver
-parameterSet.mgTolerance = 1e-10
-# Tolerance of the base grid solver
-parameterSet.baseTolerance = 1e-8
+parameterSet.maxProximalNewtonSteps = 1000
+# Initial regularization
+parameterSet.initialRegularization = 1.0
 # Measure convergence
-parameterSet.instrumented = 1
-# Norm used for stopping criterion
-# parameterSet.norm = 'H1-Semi'
-parameterSet.norm = 'infinity'
+parameterSet.instrumented = 0
 
 ############################
 #   Problem specifications
@@ -118,6 +101,16 @@ def dirichlet_indicator(x) :
 #         return False
 
 
+#boundary-values/derivative function
+def boundaryValues(x):
+    return [x[0], x[1], 0]
+
+
+# def boundaryValuesDerivative(x):
+#     return ((1,0,0),
+#             (0,1,0),
+#             (0,0,1))
+
 
 
 #############################################