Skip to content

Commit

Permalink
Mixed poisson in work
Browse files Browse the repository at this point in the history
  • Loading branch information
niravshah241 committed Jul 14, 2024
1 parent e7fa2ca commit 7a43074
Show file tree
Hide file tree
Showing 2 changed files with 5 additions and 5 deletions.
4 changes: 4 additions & 0 deletions demo/mixed_poisson_dlrbnicsx/dlrbnicsx_mixed_poisson.py
Original file line number Diff line number Diff line change
Expand Up @@ -261,13 +261,15 @@ def norm_error_u(self, u_true, u_rb):
print(sigma_sol.x.array, np.linalg.norm(sigma_sol.x.array))
print(u_sol.x.array, np.linalg.norm(u_sol.x.array))

'''
with dolfinx.io.XDMFFile(mesh.comm, "parametric_mixed_poisson/sigma.xdmf", "w") as sol_file:
sol_file.write_mesh(mesh)
sol_file.write_function(sigma_sol)
with dolfinx.io.XDMFFile(mesh.comm, "parametric_mixed_poisson/u.xdmf", "w") as sol_file:
sol_file.write_mesh(mesh)
sol_file.write_function(u_sol)
'''

# TODO more benchmarking with dolfinx implementation for correctness

Expand Down Expand Up @@ -559,6 +561,7 @@ def generate_ann_output_set(problem, reduced_problem, input_set, mode=None):
solution_sigma_error = dolfinx.fem.Function(problem_parametric._Q)
solution_sigma_error.x.array[:] = abs(solution_sigma.x.array - rb_solution_sigma.x.array)

'''
with dolfinx.io.XDMFFile(mesh.comm, "dlrbnicsx_solution/fem_sigma_online_mu.xdmf",
"w") as solution_file:
solution_file.write_mesh(mesh)
Expand All @@ -568,3 +571,4 @@ def generate_ann_output_set(problem, reduced_problem, input_set, mode=None):
"w") as solution_file:
solution_file.write_mesh(mesh)
solution_file.write_function(rb_solution_sigma)
'''
6 changes: 1 addition & 5 deletions demo/mixed_poisson_dlrbnicsx/mesh_data/mesh.py
Original file line number Diff line number Diff line change
@@ -1,10 +1,6 @@
import dolfinx
import gmsh
import sys
import math
import numpy as np

from mpi4py import MPI
import dolfinx

gmsh.initialize()
gmsh.model.add("3d_mesh")
Expand Down

0 comments on commit 7a43074

Please sign in to comment.