CSGCGALMeshGenerator3D() parameters

Issue #87 new
adam.connolly@kcl.ac.uk created an issue

Hello,

I am unable to feed parameters to the CSGCGALMeshGenerator3D() mesher. The following code gives identical meses for mesh1 and mesh2...

from dolfin import *
from mshr import *
from fenics import *

domain=Box(Point(-1,1,-1), Point(1,-1,1))

gen = CSGCGALMeshGenerator3D()
gen.parameters["facet_angle"] = 30.0

# Parameter set 1:
gen.parameters["facet_size"] = 0.5
gen.parameters["edge_size"] = 0.5

mesh1 = gen.generate(CSGCGALDomain3D(domain))
f1 =  File("mesh1.pvd")
f1 << mesh1

# Parameter set 2:
gen.parameters["facet_size"] = 1.0
gen.parameters["edge_size"] = 1.0

mesh2 = gen.generate(CSGCGALDomain3D(domain))
f2 =  File("mesh2.pvd")
f2 << mesh2

Is this a bug or am I doing something wrong? Thanks, Adam

Comments (3)

  1. adam.connolly@kcl.ac.uk reporter

    I reverted my FENICS docker image from "quay.io/fenicsproject/stable:current" to "quay.io/fenicsproject/stable:2017.1.0.r1" and this works as it should.

    Check the difference in the source to narrow down the bug.

  2. Log in to comment