Transposed a1_phys, a2_phys and the returns in pull_2d_hcurl

Author: max-models

psydac/feec/pull_push.py

@@ -121,8 +121,8 @@ def f1_logical(eta1, eta2):
         a2_phys = f2(x, y)
 
         J_T_value = F.jacobian(eta1, eta2).T
-        value_1   = J_T_value[..., 0, 0] * a1_phys + J_T_value[..., 0, 1] * a2_phys
-        return value_1
+        value_1   = J_T_value[..., 0, 0] * a1_phys.T + J_T_value[..., 0, 1] * a2_phys.T
+        return value_1.T
 
     def f2_logical(eta1, eta2):
         x, y = F(eta1, eta2)
@@ -131,8 +131,8 @@ def f2_logical(eta1, eta2):
         a2_phys = f2(x, y)
 
         J_T_value = F.jacobian(eta1, eta2).T
-        value_2   = J_T_value[..., 1, 0] * a1_phys + J_T_value[..., 1, 1] * a2_phys
-        return value_2
+        value_2   = J_T_value[..., 1, 0] * a1_phys.T + J_T_value[..., 1, 1] * a2_phys.T
+        return value_2.T
 
     return f1_logical, f2_logical