Commit 1d749739 authored by Tilman Aleman's avatar Tilman Aleman
Browse files

Merge branch 'main' of git.rwth-aachen.de:tilmandiego.aleman/reproduction

parents e76d934f ea4519ef
for ORDER in 1 2 3
for ORDER in 1 2 3 4 5
do
for C in 1 1.1 1.5
for C in 1 1.0005 1.001
do
./killing.py --refinements 5 -c $C --order $ORDER --plot --saveplot
......@@ -8,7 +8,7 @@ do
done
done
for ORDER in 1 2 3
for ORDER in 1 2 3 4 5
do
for C in "1" "0.5" "0.1"
......
......@@ -153,8 +153,9 @@ def computeEV(phi, mesh, lset_approx, deformation, order=2):
vertexdofs = BitArray(VhG.ndof)
vertexdofs[:] = False
print("doing stuff for precond")
"""
print("doing stuff for precond")
for v in mesh.vertices:
......@@ -170,10 +171,14 @@ def computeEV(phi, mesh, lset_approx, deformation, order=2):
# blocks.append(vdofs)
<<<<<<< HEAD
blocks = [set(d for el in mesh[v].elements for d in VhG.GetDofNrs(el) if freedofs[d]) for v in mesh.vertices]
=======
blocks = [set(VhG.GetDofNrs(el)) for el in VhG.Elements(VOL)]
>>>>>>> ea4519effe74102893ea15789e6c4633f4f61bbf
# blocks = [reduce(lambda x,y: x or y,[set(d for d in VhG.GetDofNrs(el) if freedofs[d])for el in mesh[v].elements]) for v in mesh.vertices]
"""
def blockcr(FEspace):
for v in FEspace.mesh.vectices:
vdofs = set()
......@@ -184,6 +189,7 @@ def computeEV(phi, mesh, lset_approx, deformation, order=2):
# blocks = pool.map(partial(calcstuff, mesh=mesh), range(len(mesh.vertices)))
vertexdofs &= VhG.FreeDofs()
"""
"""
print("done with stuff for precond")
n = Normalize(grad(lset_approx))
ntilde = GridFunction(VhG)
......
......@@ -275,9 +275,11 @@ if __name__ == "__main__":
lset_approx, deformation = refineMesh(mesh, phi, order)
vlam, ev = computeEV(phi,mesh,lset_approx, deformation, order)
print("eigenvalues:", vlam)
# vlam, ev = computeEV2(Vh, VhG, a, b, pre)
for k in range(3):
vlams[k].append(vlam[k])
print("vlamk: ", vlam[k])
if i < 2 and i != 0:
symmetries = 0
for k in range(3):
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment