Add grid parameter support to discretize_domain function

- Add grid parameter to discretize_domain() in discretization.py
- Add grid parameter to Geometry.from_topological_domain() in geometry.py

Author: ArasuCandassamy

psydac/api/discretization.py

@@ -553,8 +553,7 @@ def discretize_space(V, domain_h, *, degree=None, multiplicity=None, knots=None,
 
 
 #==============================================================================
-def discretize_domain(domain, *, filename=None, ncells=None, periodic=None, comm=None, mpi_dims_mask=None):
-
+def discretize_domain(domain, *, filename=None, ncells=None, periodic=None, comm=None, mpi_dims_mask=None, grid=None):
     if comm is not None:
         # Create a copy of the communicator
         comm = comm.Dup()
@@ -569,7 +568,7 @@ def discretize_domain(domain, *, filename=None, ncells=None, periodic=None, comm
         return Geometry(filename=filename, comm=comm)
 
     elif ncells:
-        return Geometry.from_topological_domain(domain, ncells, periodic=periodic, comm=comm, mpi_dims_mask=mpi_dims_mask)
+        return Geometry.from_topological_domain(domain, ncells, periodic=periodic, comm=comm, mpi_dims_mask=mpi_dims_mask, grid=grid)
 
 #==============================================================================
 def discretize(a, *args, **kwargs):