Andreas Dedner (0994e799) at 29 Feb 09:50
remove .dev versioning requirements for dune-grid
Andreas Dedner (cc3c4d47) at 29 Feb 09:47
remove .dev versioning requirements for dune-grid
Andreas Dedner (a1010a33) at 24 Feb 12:02
increase required cmake version
@samuel.burbulla are you still occasional looking into this
Hi,
I am currently trying to do something with interface identification. I had been hinted towards MMesh which has been helpful.
While I found out how to do what I wanted to do, I also encountered a bug whereby adding two new vertices within a single element via refineEdge
only adds the first.
This can be seen from the relatively minimal code:
from dune.mmesh import mmesh
from dune.fem.view import adaptiveLeafGridView as adaptive
from dune.fem import parameter, adapt
import numpy as np
from ufl import SpatialCoordinate, conditional
parameter.append( { "fem.adaptation.method": "callback" } )
dim = 2
vertices = np.array( [(-1,-1), (1,-1), (1,1), (-1,1),(0,0)] )
triangles = np.array( [(0,1,4), (1,2,4), (2,4,3), (3,4,0)] )
gridView = mmesh({"vertices": vertices, "simplices": triangles}, dim)
gridView = adaptive( gridView )
hgrid = gridView.hierarchicalGrid
from dune.fem.space import lagrange
spaceScalar = lagrange(gridView, order = 1); x = SpatialCoordinate(spaceScalar)
phi = spaceScalar.interpolate(1- x[0]**2 - x[1]**2 ,name = "phi")
for e in gridView.elements:
for k, edge in enumerate(e.subEntities(1)):
hgrid.refineEdge(e, k)
adapt([phi])
We see that not every edge has been bisected here.
For my project, I do not wish to refine every edge, I have some criteria based on phi. My solution is to loop over the edges and adapt twice - perhaps something like that might be of use.
Best,
Philip
Samuel Burbulla (9b19a653) at 13 Apr 07:01
Merge pull request #32 from samuelburbulla/tkoch-master-patch-10131
... and 3 more commits
Closed in favour of https://github.com/samuelburbulla/dune-mmesh/pull/32
Samuel Burbulla (d0a3ba4f) at 03 Jan 14:51
Feature/coding style (#31)
Samuel Burbulla (b52ede6f) at 13 Dec 08:18
Add dune-grid 2.10.0.dev0 to python requires.
It seems like the configure step runs through but then the build step wants to build dune-fem even though it hasn't been found.
dunecontrol --opts=$DUNE_OPTS_FILE configure
dunecontrol --opts=$DUNE_OPTS_FILE make -j8
This is when building a Dockerfile (https://git.iws.uni-stuttgart.de/dumux-repositories/dumux-docker-ci/-/blob/master/full/ubuntu/Dockerfile) with Dumux and dune-mmesh (1.4). See failing pipeline here: (https://git.iws.uni-stuttgart.de/dumux-repositories/dumux-docker-ci/-/jobs/142789).
Is dune-fem a requirement for dune-mmesh now?
the issue seems to be resolved thanks!
Samuel Burbulla (3e0a4908) at 12 Dec 14:09
Add dune-grid with version to python requires.
Samuel Burbulla (b67c7f97) at 09 Dec 20:01
Merge pull request #30 from samuelburbulla/feature/fix-docker
... and 2 more commits
Samuel Burbulla (00418a84) at 09 Dec 13:00
Fixup duplicate tests.