CSG mesh generation cell size broken

Ok. Do you have a simple example that reproduces this?

On 14 March 2013 09:12, Benjamin Kehlet <email address hidden> wrote:
> Ok. Do you have a simple example that reproduces this?
>

from dolfin import *
sphere1 = Sphere(Point(0.0, 0.0, 0.0), 1.0)
sphere0 = Sphere(Point(0.0, 0.0, 0.0), 0.2)
g3d = sphere1 - sphere0;
mesh3d = Mesh(g3d, 32);

and put

cout << "*** Cell size: " << r << ", " << cell_size << ", " <<
(double) parameters["edge_size"] << endl;

on line 150 of CSGCGALMeshGenerator3D.cpp.

Sometimes 'r' is roughly right, sometimes it's zero and sometimes it's huge.

Garth

> ** Changed in: dolfin
> Assignee: (unassigned) => Benjamin Kehlet (benjamik)
>
> --
> You received this bug notification because you are a member of DOLFIN
> Core Team, which is subscribed to DOLFIN.
> https://bugs.launchpad.net/bugs/1154991
>
> Title:
> CSG mesh generation cell size broken
>
> To manage notifications about this bug go to:
> https://bugs.launchpad.net/dolfin/+bug/1154991/+subscriptions

2013/3/14 Garth Wells <email address hidden>

> On 14 March 2013 09:12, Benjamin Kehlet <email address hidden>
> wrote:
> > Ok. Do you have a simple example that reproduces this?
> >
>
>
> from dolfin import *
> sphere1 = Sphere(Point(0.0, 0.0, 0.0), 1.0)
> sphere0 = Sphere(Point(0.0, 0.0, 0.0), 0.2)
> g3d = sphere1 - sphere0;
> mesh3d = Mesh(g3d, 32);
>
>
> and put
>
> cout << "*** Cell size: " << r << ", " << cell_size << ", " <<
> (double) parameters["edge_size"] << endl;
>
> on line 150 of CSGCGALMeshGenerator3D.cpp.
>
> Sometimes 'r' is roughly right, sometimes it's zero and sometimes it's
> huge.
>
> Garth
>

I'm not able to reproduce this on my (32 bit) laptop. I'll try tomorrow on
another computer.

Benjamin

>
>
> > ** Changed in: dolfin
> > Assignee: (unassigned) => Benjamin Kehlet (benjamik)
> >
> > --
> > You received this bug notification because you are a member of DOLFIN
> > Core Team, which is subscribed to DOLFIN.
> > https://bugs.launchpad.net/bugs/1154991
> >
> > Title:
> > CSG mesh generation cell size broken
> >
> > To manage notifications about this bug go to:
> > https://bugs.launchpad.net/dolfin/+bug/1154991/+subscriptions
>
> --
> You received this bug notification because you are a bug assignee.
> https://bugs.launchpad.net/bugs/1154991
>
> Title:
> CSG mesh generation cell size broken
>
> Status in DOLFIN:
> New
>
> Bug description:
> The 3D CSG mesh generation gives erratic results when the cell size if
> given. The problem appears to be with the line
>
> const double r = PolyhedronUtils::getBoundingSphereRadius(p);
>
> in CSGCGALMeshGenerator3D::generate(Mesh& mesh). The return value r is
> not deterministic, and often returns extremely large values, e.g.
> 10^54
>
> To manage notifications about this bug go to:
> https://bugs.launchpad.net/dolfin/+bug/1154991/+subscriptions
>