bugs - Same region, different meshes

I'm playing around with ToElementMesh and I noticed the following behaviour that puzzles me. Consider the two regions

<< NDSolve`FEM`

Ω1 = 
  Polygon[{{-1/2, -1}, {-1/2, 1}, {1/2, 1}, {1/2, -1}}];
Ω2 = Rectangle[{-1/2, -1}, {1/2, 1}];

RegionEqual[Ω1, Ω2] returns True, but

mesh1 = ToElementMesh[Ω1, 
   "BoundaryMeshGenerator" -> "Continuation", MaxCellMeasure -> 0.05, 
   "MaxBoundaryCellMeasure" -> 0.05];
mesh2 = ToElementMesh[Ω2, 
   "BoundaryMeshGenerator" -> "Continuation", MaxCellMeasure -> 0.05, 
   "MaxBoundaryCellMeasure" -> 0.05];
Show[RegionPlot[Ω1], mesh1["Wireframe"]]
Show[RegionPlot[Ω2], mesh2["Wireframe"]]

yields two different meshes:

Questions:

Why this behaviour? From the documentation I couldn't find an explanation for this.

How can I control the mesh generation to get a uniform behaviour?

Answer

This is a shortcoming in ToElementMesh. I filed this as a bug, but it will take a bit of time to fix.

As a workaround you could use

mesh1 = 
  ToElementMesh[
    BoundaryDiscretizeRegion[Ω1, 
      MaxCellMeasure -> {"Edges" -> 0.01}], 
    "BoundaryMeshGenerator" -> "Continuation", MaxCellMeasure -> 0.05, 
    "MaxBoundaryCellMeasure" -> 0.01];
Show[RegionPlot[Ω1], mesh1["Wireframe"]]

Comments

plotting - How to draw lines between specified dots on ListPlot?

I would like to create a plot where I have unconnected dots and some connected. So far, I have figured out how to draw the dots. My code is the following: ListPlot[{{1, 1}, {2, 2}, {3, 3}, {4, 4}, {1, 4}, {2, 5}, {3, 6}, {4, 7}, {1, 7}, {2, 8}, {3, 9}, {4, 10}, {1, 10}, {2, 11}, {3, 12}, {4,13}, {2.5, 7}}, Ticks -> {{1, 2, 3, 4}, None}, AxesStyle -> Thin, TicksStyle -> Directive[Black, Bold, 12], Mesh -> Full] I have thought using ListLinePlot command, but I don't know how to specify to the command to draw only selected lines between the dots. Do have any suggestions/hints on how to do that? Thank you. Answer One possibility would be to use Epilog with Line : ListPlot[ {{1, 1}, {2, 2}, {3, 3}, {4, 4}, {1, 4}, {2, 5}, {3, 6}, {4, 7}, {1, 7}, {2, 8}, {3, 9}, {4, 10}, {1, 10}, {2, 11}, {3, 12}, {4, 13}, {2.5, 7}}, Ticks -> {{1, 2, 3, 4}, None}, AxesStyle -> Thin, TicksStyle -> Directive[Black, Bold, 12], Mesh -> Full, Epilog -> { Line[ ...

dynamic - How can I make a clickable ArrayPlot that returns input?

I would like to create a dynamic ArrayPlot so that the rectangles, when clicked, provide the input. Can I use ArrayPlot for this? Or is there something else I should have to use? Answer ArrayPlot is much more than just a simple array like Grid : it represents a ranged 2D dataset, and its visualization can be finetuned by options like DataReversed and DataRange . These features make it quite complicated to reproduce the same layout and order with Grid . Here I offer AnnotatedArrayPlot which comes in handy when your dataset is more than just a flat 2D array. The dynamic interface allows highlighting individual cells and possibly interacting with them. AnnotatedArrayPlot works the same way as ArrayPlot and accepts the same options plus Enabled , HighlightCoordinates , HighlightStyle and HighlightElementFunction . data = {{Missing["HasSomeMoreData"], GrayLevel[ 1], {RGBColor[0, 1, 1], RGBColor[0, 0, 1], GrayLevel[1]}, RGBColor[0, 1, 0]}, {GrayLevel[0], GrayLevel...

functions - What best practices or performance considerations are there for choosing between Cases, Position, Pick and Select?

Cases , Select , Pick and Position each have different syntaxes and purposes, but there are times when you can express the same calculation equivalently using either of them. So with this input: test = RandomInteger[{-25, 25}, {20, 2}] {{-15, 13}, {-8, 16}, {-8, -19}, {7, 6}, {-21, 9}, {-3, -25}, {21, -18}, {4, 4}, {2, -2}, {-24, 8}, {-17, -8}, {4, -18}, {22, -24}, {-4, -3}, {21, 0}, {19, 18}, {-23, -8}, {23, -25}, {14, -2}, {-1, -13}} You can get the following equivalent results: Cases[test, {_, _?Positive}] {{-15, 13}, {-8, 16}, {7, 6}, {-21, 9}, {4, 4}, {-24, 8}, {19, 18}} Select[test, #[[2]] > 0 &] {{-15, 13}, {-8, 16}, {7, 6}, {-21, 9}, {4, 4}, {-24, 8}, {19, 18}} Pick[test, Sign[test[[All, 2]] ], 1] {{-15, 13}, {-8, 16}, {7, 6}, {-21, 9}, {4, 4}, {-24, 8}, {19, 18}} test[[Flatten@Position[test[[All, 2]], _?Positive] ]] {{-15, 13}, {-8, 16}, {7, 6}, {-21, 9}, {4, 4}, {-24, 8}, {19, 18}} Are there performance or other considerations that should guide which you shou...

Blog

Search This Blog