Plotting surfaces doesn't work

[Garth Wells]

It's not possible to view the mesh or contour field properly in 3D case. To reproduce, run elasticity demo

  demo/undocumented/elasticity

and press 'S' on plot screen. It does nothing.

[Anders Logg (logg)]

Yes it does. It toggles the wireframe/solid mode for the arrows.

--
Anders

On Tue, Aug 07, 2012 at 10:39:47PM -0000, Garth Wells wrote:
> Public bug reported:
>
> It's not possible to view the mesh or contour field properly in 3D case.
> To reproduce, run elasticity demo
>
> demo/undocumented/elasticity
>
> and press 'S' on plot screen. It does nothing.
>
> ** Affects: dolfin
> Importance: Undecided
> Status: New
>

[Joachim Haga (jobh)]

I have added a shortcut key 'm' that views the mesh (wireframe) also for glyps (arrows) mode. Is that sufficient, or do we also need a shortcut key to toggle between warped (displaced) mesh and glyphs (arrows)?

[Joachim Haga (jobh)]

Not sure if the "fix" of showing mesh along with glyphs (present in trunk) is sufficient to close this bug. Please advise.

[Anders Logg (logg)]

I think the plotting for the elasticity demo looks perfect now. Open new bug if problem persists.

[Joachim Haga (jobh)]

Just a couple of final notes:

1) This plot probably looks better when plotted as displacement (mode="displacement").

2) I see a problem with glyph scaling, which seems to be caused by mesh->hmax() returning a very high value on the gear mesh. Maybe a numerical precision thing. This:

In [1]: from dolfin import *
In [2]: m=Mesh("gear.xml.gz")
In [3]: c=Cell(m,18847)
In [4]: c.diameter()
Out[4]: 36.0014830310018
In [5]: c.volume()
Out[5]: 7.373890632542822e-11

[Anders Logg (logg)]

Yes, very suspicious.

Can you open a new bug report for this? (Fixing the gear grid).

--
Anders

On Mon, Nov 19, 2012 at 02:31:36PM -0000, Joachim Haga wrote:
> Just a couple of final notes:
>
> 1) This plot probably looks better when plotted as displacement
> (mode="displacement").
>
> 2) I see a problem with glyph scaling, which seems to be caused by
> mesh->hmax() returning a very high value on the gear mesh. Maybe a
> numerical precision thing. This:
>
> In [1]: from dolfin import *
> In [2]: m=Mesh("gear.xml.gz")
> In [3]: c=Cell(m,18847)
> In [4]: c.diameter()
> Out[4]: 36.0014830310018
> In [5]: c.volume()
> Out[5]: 7.373890632542822e-11
>