From d869c6bd24da3240efc34400b02c3e701332253c Mon Sep 17 00:00:00 2001
From: Valentin Helml <valentin.helml@tuwien.ac.at>
Date: Tue, 14 Mar 2023 14:06:32 +0100
Subject: [PATCH] Update p_laplace.m
---
experiments/p_laplace.m | 8 +++-----
1 file changed, 3 insertions(+), 5 deletions(-)
diff --git a/experiments/p_laplace.m b/experiments/p_laplace.m
index 9e6c2fc..e4168f2 100644
--- a/experiments/p_laplace.m
+++ b/experiments/p_laplace.m
@@ -5,7 +5,7 @@ lambda = 0.1;
deltaZ = 0.1;
linearization = 'kacanov';
geometry = 'Lshape';
-p = 2;
+p = 3;
%% initialisation
[nElem, eta, E] = deal(zeros(1,1000));
@@ -21,9 +21,7 @@ Du = Gradient(u);
g = Constant(mesh, 1);
[blf, lf] = setProblemData(mesh, fes, Du, g, linearization, p);
eDensity = CompositeFunction(@(t) muIntegral(vectorProduct(t, t),p), Du);
-u.setData(0);
-%u.setFreeData(1);
-u.setData(rand(1,length(u.data)));
+u.setFreeData(rand(1,length(getFreeDofs(u.fes))));
%% nested iteration
P = LoFeProlongation(fes);
@@ -128,7 +126,7 @@ function updateDataU(u, deltaZ, A, F, linearization)
u.setData(u.data + deltaZ*v.data);
case "kacanov"
%u.setData(A \ F);
- u.setFreeData(A(freeDofs,freeDofs) \ F(freeDofs))
+ u.setFreeData(A(freeDofs,freeDofs) \ F(freeDofs));
case "newton"
v.setFreeData(A(freeDofs,freeDofs) \ F(freeDofs));
u.setData(u.data + v.data);
--
GitLab