Add evaluation kernels for single FemField

psydac/core/field_evaluation_kernels.py

@@ -4,6 +4,37 @@
 # =============================================================================
 # Field evaluation functions
 # =============================================================================
+# -----------------------------------------------------------------------------
+# 0: Evaluation of single 3d field at single point
+# -----------------------------------------------------------------------------
+@template(name='T', types=[float, complex])
+def eval_field_3d_once(local_coeffs : 'T[:,:,:]', 
+                       local_bases_0: 'float[:]',
+                       local_bases_1: 'float[:]',
+                       local_bases_2: 'float[:]'):
+    """
+    Parameters
+    ----------
+    local_coeffs: ndarray of floats
+        Active (local) coefficients of the fields in all directions
+
+    local_bases: list of ndarrays
+        Active (local) 1D-basis functions values at the point of evaluation. 
+    """
+    n1, n2, n3 = local_coeffs.shape
+
+    res = local_coeffs[0, 0, 0] - local_coeffs[0, 0, 0]
+    for i1 in range(n1):
+        for i2 in range(n2):
+            for i3 in range(n3):
+                res += (local_coeffs[i1, i2, i3] *
+                        local_bases_0[i1] *
+                        local_bases_1[i2] *
+                        local_bases_2[i3])
+
+    return res
+
+
 # -----------------------------------------------------------------------------
 # 1: Regular tensor grid without weight
 # -----------------------------------------------------------------------------

psydac/core/tests/test_field_evaluation_kernel.py

@@ -642,3 +642,7 @@ def test_pushforwards_hcurl(ldim):
         pushforward_3d_hcurl(field_to_push, inv_jacobians, out)
 
     assert np.allclose(expected, out, atol=ATOL, rtol=RTOL)
+
+if __name__ == '__main__':
+    test_regular_jacobians('bent_pipe.h5', 2)
+    test_regular_evaluations([np.sort(np.concatenate((np.zeros(3), np.random.random(9), np.ones(3)))) for i in range(2)], 2, [2] * 2,2)