From c843f5a8b5d0d07fed7c3a29d69fffe7f656d380 Mon Sep 17 00:00:00 2001
From: Andreas Kloeckner <inform@tiker.net>
Date: Sun, 1 Nov 2015 22:10:38 -0600
Subject: [PATCH] Move warped-rect generation to meshmode
---
examples/geometry.py | 20 ++------------------
1 file changed, 2 insertions(+), 18 deletions(-)
diff --git a/examples/geometry.py b/examples/geometry.py
index 51831b4e..e10254bf 100644
--- a/examples/geometry.py
+++ b/examples/geometry.py
@@ -33,24 +33,8 @@ def main(write_output=True):
cl_ctx = cl.create_some_context()
queue = cl.CommandQueue(cl_ctx)
- from meshmode.mesh.generation import generate_regular_rect_mesh
- mesh = generate_regular_rect_mesh(a=(-0.5, -0.5), b=(0.5, 0.5),
- n=(6, 6), order=4)
-
- def m(x):
- result = np.empty_like(x)
- result[0] = (
- 1.5*x[0] + np.cos(x[0])
- + 0.1*np.sin(10*x[1]))
- result[1] = (
- 0.05*np.cos(10*x[0])
- + 1.3*x[1] + np.sin(x[1]))
- if len(x) == 3:
- result[2] = x[2]
- return result
-
- from meshmode.mesh.processing import map_mesh
- mesh = map_mesh(mesh, m)
+ from meshmode.mesh.generation import generate_warped_rect_mesh
+ mesh = generate_warped_rect_mesh(dim=2, order=4, n=6)
discr = Discretization(cl_ctx, mesh, order=4)
--
GitLab