# Inconsistent results when using boundary operator for maxwell's equations

**URL:** <https://bempp.discourse.group/t/inconsistent-results-when-using-boundary-operator-for-maxwells-equations/152>\
**Category:** Uncategorized\
**Created:** [15 December 2022 15:02 UTC](https://bempp.discourse.group/t/inconsistent-results-when-using-boundary-operator-for-maxwells-equations/152 "2022-12-15T15:02:41Z")\
**Posts on this page:** 1\
**Page:** 1

<div class="post-metadata">

**Author:** ![chris](https://avatars.discourse-cdn.com/v4/letter/c/958977/32.png) [@chris](https://bempp.discourse.group/u/chris)\
**Post date:** [15 December 2022 15:02 UTC](https://bempp.discourse.group/t/inconsistent-results-when-using-boundary-operator-for-maxwells-equations/152/1 "2022-12-15T15:02:41Z")

</div>

Hello,

The results that I get when using the solver are inconsistent. I have the following two examples. I use the following code for the two examples:

```python
import bempp.api
import numpy as np

# example 1: plate with four RWGs / triangles
vertices = np.array([[0, 0.004, 0.004, 0, 0.002],[0, 0, 0.004, 0.004, 0.002],[0,0,0,0,0]],dtype=np.float64)
elements = np.array([[0,1,2,3],[1,2,3,0],[4,4,4,4]],dtype=np.uint32)
grid = bempp.api.Grid(vertices, elements)

# example 2: one RWG / two triangles of the plate of example 1
#vertices = np.array(
# [[0.0,0.004,0.004,0.002],[0.0,0.0,0.004,0.002],[0.0,0.0,0.0,0.0]],
# dtype=np.float64)
#elements = np.array(
# [[0,1],[1,2],[3,3]], dtype=np.uint32)
#grid = bempp.api.Grid(vertices, elements)

E0 = np.array([0.0, 1.0, 0.0]) # polarization of electric field
kv = np.array([0.0, 0, -314.1593]) # wave vector
k = np.linalg.norm(kv) # wave number

div_space = bempp.api.function_space(grid, "RWG", 0)
curl_space = bempp.api.function_space(grid, "SNC", 0)

def incident_field(x):
    return E0*np.exp(1j*np.dot(kv,x))

@bempp.api.complex_callable
def tangential_trace(x, n, domain_index, result):
    incident_field = E0*np.exp(1j*np.dot(kv,x))
    result[:] = np.cross(incident_field, n)

trace_fun = bempp.api.GridFunction(div_space, fun=tangential_trace, dual_space=curl_space)

from bempp.api.operators.boundary import maxwell
elec = maxwell.electric_field(div_space, div_space, curl_space, k)

print(trace_fun._projections) # values of incident field term
z = elec.weak_form()
print(z._impl) # impedance matrix

```

Example 1 contains a plate of 4 triangles / RWGs. This gives a 4x4 BEM matrix. Example 2 contains one RWG of the plate of example 1. The diagonal elements of the 4x4 BEM matrix in example 1 have the value:  
1.44531462e-07+2.71073075e-06j  
The single element BEM matrix in example 2 has the value:  
1.60467455e-07+1.40209657e-05j  
These values are not the same. I think the values of the diagonal elements of the BEM matrix in example 1 should be the same as the value of the single element of the BEM matrix in example 2.

I hope I formatted the code correctly.

Best regards,  
Christiaan