# Path: sucemfem/Consts.py
# Z0 = mu0*c0                              


import sys
import numpy as N
import os
import dolfin as dol
import scipy.sparse
import point_source
from dolfin import dot, cross, curl, inner, dx, ds
from sucemfem.Consts import eps0, mu0, c0
from scipy.sparse import csr_matrix
from pyamg import smoothed_aggregation_solver
from numpy import intc
## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
## GNU General Public License for more details.
##
## You should have received a copy of the GNU General Public License
## along with SUCEM. If not, see <http://www.gnu.org/licenses/>. 
##
## Contact: cemagga@gmail.com 
# Authors:
# Neilen Marais <nmarais@gmail.com>
from __future__ import division

sys.path.append('../../')

reload(point_source)

# Define material and frequency
eps_r = 1;
mu_r = 1;
freq = 1e9;
source_coord = [0.5,0.5,0.5]
source_value = N.array([0,0,1.])

# Define mesh
mesh = dol.UnitCube(11,11,11)

# Define function space
order = 1
V = dol.FunctionSpace(mesh, "Nedelec 1st kind H(curl)", order)

# Define basis and bilinear form

# Path: sucemfem/Consts.py
# Z0 = mu0*c0                              


# Path: sandbox/ABC/abc_test.py
import sys
import numpy as N
import os
import dolfin as dol
import scipy.sparse
import point_source
from dolfin import dot, cross, curl, inner, dx, ds
from sucemfem.Consts import eps0, mu0, c0
from scipy.sparse import csr_matrix
from pyamg import smoothed_aggregation_solver
from numpy import intc

## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
## GNU General Public License for more details.
##
## You should have received a copy of the GNU General Public License
## along with SUCEM. If not, see <http://www.gnu.org/licenses/>. 
##
## Contact: cemagga@gmail.com 
# Authors:
# Neilen Marais <nmarais@gmail.com>
from __future__ import division

sys.path.append('../../')

reload(point_source)

# Define material and frequency
eps_r = 1;
mu_r = 1;
freq = 1e9;
source_coord = [0.5,0.5,0.5]
source_value = N.array([0,0,1.])

# Define mesh
mesh = dol.UnitCube(11,11,11)

# Define function space
order = 1
V = dol.FunctionSpace(mesh, "Nedelec 1st kind H(curl)", order)

# Define basis and bilinear form

k_0 = 2*N.pi*freq/c0