Add parent to para scaling to Galerkin IGA example

markriegler · j042 · commit fc2c83deebf3 · 2025-04-05T10:08:24.000+02:00
diff --git a/examples/iga/galerkin_laplace_problem.py b/examples/iga/galerkin_laplace_problem.py
@@ -95,6 +95,8 @@ def map_positions(positions, x_min, x_max, y_min, y_max):
             ukv[1][j + 1],
         )
         det_jacs = np.linalg.det(geometry.jacobian(mapped_positions))
+        # Scale jacobian by element size
+        det_jacs *= (ukv[0][i + 1] - ukv[0][i]) * (ukv[1][j + 1] - ukv[1][j])
 
         # RHS
         # q : quadrature point | d: dim

Original file line number	Diff line number	Diff line change
`@@ -95,6 +95,8 @@ def map_positions(positions, x_min, x_max, y_min, y_max):`
`95`	`95`	`ukv[1][j + 1],`
`96`	`96`	`)`
`97`	`97`	`det_jacs = np.linalg.det(geometry.jacobian(mapped_positions))`
	`98`	`+ # Scale jacobian by element size`
	`99`	`+ det_jacs = (ukv[0][i + 1] - ukv[0][i]) (ukv[1][j + 1] - ukv[1][j])`
`98`	`100`
`99`	`101`	`# RHS`
`100`	`102`	`# q : quadrature point \| d: dim`