From d763bb33360e7dc6ecd3f8d0f2edb6bac1d91f42 Mon Sep 17 00:00:00 2001
From: Olivier Ali <olivier.ali@inria.fr>
Date: Thu, 12 Jan 2023 21:22:09 +0100
Subject: [PATCH] Minor update (comments) in the CustomDomain class.
---
src/bvpy/domains/custom_domain.py | 8 +++-----
1 file changed, 3 insertions(+), 5 deletions(-)
diff --git a/src/bvpy/domains/custom_domain.py b/src/bvpy/domains/custom_domain.py
index 553a49f..a3a0ae0 100644
--- a/src/bvpy/domains/custom_domain.py
+++ b/src/bvpy/domains/custom_domain.py
@@ -83,14 +83,12 @@ class CustomDomain(AbstractDomain, BuiltInModel):
points = np.array(points)
cells = np.array(cells)
- # self.mesh = fe.Mesh()
- self.mesh = None
self.dim = points.shape[1]
- print(f'self.dim = {self.dim}')
-
# editor = fe.MeshEditor()
+
+ # OA: What is the purpose of the following block ?
if self._cell_type == 'line':
topo_dim = 1 # OA: changed variable name top -> topo_dim
# OA: what is the use of this variable ?
@@ -98,6 +96,7 @@ class CustomDomain(AbstractDomain, BuiltInModel):
topo_dim = 2
elif self._cell_type == 'tetra':
topo_dim = 3
+
# editor.open(self.mesh, self._cell_type, topo_dim, self.dim)
# editor.init_vertices(len(points))
# editor.init_cells(len(cells))
@@ -107,7 +106,6 @@ class CustomDomain(AbstractDomain, BuiltInModel):
# editor.close()
- print(f'pts used to create mesh = {points[:, :self.dim]}')
_cell = ufl.Cell(self._cell_type, geometric_dimension=self.dim)
_domain = ufl.Mesh(ufl.VectorElement("CG", _cell, 1))
--
GitLab