Merge pull request #49 from gridap/release-0.6

Release 0.6.1

NEWS.md

@@ -4,6 +4,12 @@ All notable changes to this project will be documented in this file.
 The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/),
 and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
 
+## [0.6.1] - 2021-06-24
+
+### Added
+- Method `AnalyticalGeometry(::Function)` to simplify the setup of geometries from user-defined functions.
+- `export AnalyticalGeometry` in `GridapEmbedded`.
+
 ## [0.6.0] - 2021-01-12
 
 ### Changed

Project.toml

@@ -1,7 +1,7 @@
 name = "GridapEmbedded"
 uuid = "8838a6a3-0006-4405-b874-385995508d5d"
 authors = ["Francesc Verdugo <[email protected]>"]
-version = "0.6.0"
+version = "0.6.1"
 
 [deps]
 AbstractTrees = "1520ce14-60c1-5f80-bbc7-55ef81b5835c"

src/Exports.jl

@@ -20,6 +20,7 @@ end
 @publish Interfaces OUT
 
 @publish LevelSetCutters LevelSetCutter
+@publish LevelSetCutters AnalyticalGeometry
 @publish LevelSetCutters popcorn
 @publish LevelSetCutters doughnut
 @publish LevelSetCutters tube

src/LevelSetCutters/AnalyticalGeometries.jl

@@ -14,6 +14,11 @@ struct BoundingBox{D,T}
   pmax::Point{D,T}
 end
 
+function AnalyticalGeometry(f::Function)
+  tree = Leaf((f,string(nameof(f)),nothing))
+  AnalyticalGeometry(tree)
+end
+
 # Factories
 
 function doughnut(R,r;x0=zero(Point{3,typeof(R)}),name="doughnut")

src/LevelSetCutters/LevelSetCutters.jl

@@ -31,6 +31,7 @@ using Gridap.Polynomials
 using Gridap.Visualization
 
 export LevelSetCutter
+export AnalyticalGeometry
 export doughnut
 export popcorn
 export tube

test/LevelSetCuttersTests/AnalyticalGeometriesTests.jl

@@ -18,6 +18,12 @@ end
 
 tmpdir() do d
 
+geo0 = AnalyticalGeometry(x->.9*x[1]^2+1.1*x[2]^2-0.6)
+model = CartesianDiscreteModel((-1,1,-1,1),(10,10))
+cutgeo = cut(model,geo0)
+Γ = EmbeddedBoundary(cutgeo)
+writevtk(Γ,joinpath(d,"user"))
+
 geo0 = popcorn()
 test_geometry(geo0)