fixed missing argument in set_up_bulk_ghost_penalty_lhs

gridap · Nov 28, 2024 · d82cd53 · d82cd53
1 parent 6a60233
commit d82cd53
Show file tree

Hide file tree

Showing 3 changed files with 9 additions and 6 deletions.
diff --git a/bulk_ghost_penalty_canvas.jl b/bulk_ghost_penalty_canvas.jl
@@ -234,7 +234,7 @@ for (i,cells) in enumerate(aggregate_to_local_cells)
     end
 end
 
-function set_up_bulk_ghost_penalty_lhs(aggregates_bounding_box_model,
+function set_up_bulk_ghost_penalty_lhs(aggregate_to_local_cells,
                                        agg_cells_to_aggregate,
                                        ref_agg_cell_to_ref_bb_map,
                                        dΩagg_cells,
@@ -262,14 +262,14 @@ function set_up_bulk_ghost_penalty_lhs(aggregates_bounding_box_model,
     lazy_map(ass_lhs_map,aggregate_to_local_cells)
 end 
 
-p_lhs=set_up_bulk_ghost_penalty_lhs(aggregates_bounding_box_model,
+p_lhs=set_up_bulk_ghost_penalty_lhs(aggregate_to_local_cells,
                                   agg_cells_to_aggregate,
                                   ref_agg_cell_to_ref_bb_map,
                                   dΩagg_cells,
                                   pbb,
                                   qbb)
 
-u_lhs=set_up_bulk_ghost_penalty_lhs(aggregates_bounding_box_model,
+u_lhs=set_up_bulk_ghost_penalty_lhs(aggregate_to_local_cells,
                                   agg_cells_to_aggregate,
                                   ref_agg_cell_to_ref_bb_map,
                                   dΩagg_cells,

diff --git a/bulk_ghost_penalty_canvas_div.jl b/bulk_ghost_penalty_canvas_div.jl
@@ -13,6 +13,7 @@ order = 2
 uex(x) = -VectorValue(2*x[1],2*x[2])
 pex(x) = (x[1]^2 + x[2]^2)
 divuex(x) = -4.0
+
 # Select geometry
 R = 0.2
 geom = disk(R, x0=Point(0.5,0.5))
@@ -84,14 +85,14 @@ vbb=get_fe_basis(Vbb)
 
 aggregate_to_local_cells=setup_aggregate_to_local_cells(aggregate_to_cells)
 
-p_lhs=set_up_bulk_ghost_penalty_lhs(aggregates_bounding_box_model,
+p_lhs=set_up_bulk_ghost_penalty_lhs(aggregate_to_local_cells,
                                   agg_cells_to_aggregate,
                                   ref_agg_cell_to_ref_bb_map,
                                   dΩagg_cells,
                                   pbb,
                                   qbb)
 
-u_lhs=set_up_bulk_ghost_penalty_lhs(aggregates_bounding_box_model,
+u_lhs=set_up_bulk_ghost_penalty_lhs(aggregate_to_local_cells,
                                   agg_cells_to_aggregate,
                                   ref_agg_cell_to_ref_bb_map,
                                   dΩagg_cells,

diff --git a/bulk_ghost_penalty_stab_tools.jl b/bulk_ghost_penalty_stab_tools.jl
@@ -112,6 +112,8 @@ function div_penalty_stabilization_collect_cell_matrix(agg_cells_to_aggregate,
                                                        U_agg_cells_to_aggregate_dof_ids,
                                                        h_U,
                                                        γ)
+    Ωagg_cells=dΩagg_cells.quad.trian
+
     ## +DIVU STABILIZATION 
     ## Compute ∫( (div_dv)*(div_du-div_du_l2_proj_agg_cells))*dΩ_agg_cells (simplified, equivalent version)
     # Manually set up the arrays that collect_cell_matrix would return automatically
@@ -212,7 +214,7 @@ function div_penalty_stabilization_collect_cell_matrix(agg_cells_to_aggregate,
     w,r,c
 end
 
-function set_up_bulk_ghost_penalty_lhs(aggregates_bounding_box_model,
+function set_up_bulk_ghost_penalty_lhs(aggregate_to_local_cells,
     agg_cells_to_aggregate,
     ref_agg_cell_to_ref_bb_map,
     dΩagg_cells,