2d solver A matrix.

.gitignore

@@ -197,4 +197,11 @@ Thumbs.db
 !examples/read_and_plot_mesh/mesh/*
 !dg_maxwell/tests/mesh/*
 !dg_maxwell/tests/msh_parser/mesh/*
+examples/results
 !docs/quazar_logo.png
+!dg_maxwell/tests/wave_equation_2d/files/*
+
+*.kate-swp
+
+*.ipynb_checkpoints*
+

.travis.yml

@@ -54,4 +54,6 @@ install:
 
 script:
     - python -c "import dg_maxwell"
-    - pytest --verbose -r P --color=yes
+    - pytest --verbose -r P --color=yes --cov dg_maxwell
+after_success:
+  - coveralls

README.md

@@ -3,7 +3,7 @@
 [![Build Status](https://travis-ci.org/QuazarTech/DG_Maxwell.svg?branch=master)](https://travis-ci.org/QuazarTech/DG_Maxwell)
 [![Codacy Badge](https://api.codacy.com/project/badge/Grade/e8733cdbf1454af0ac35ae5b2d017d9f)](https://www.codacy.com/app/aman2official/DG_Maxwell_2?utm_source=github.com&utm_medium=referral&utm_content=QuazarTech/DG_Maxwell&utm_campaign=Badge_Grade)
 [![Documentation Status](http://readthedocs.org/projects/dg-maxwell/badge/?version=latest)](http://dg-maxwell.readthedocs.io/en/latest/?badge=latest)
-
+[![Coverage Status](https://coveralls.io/repos/github/QuazarTech/DG_Maxwell/badge.svg?branch=master)](https://coveralls.io/github/QuazarTech/DG_Maxwell?branch=master)
 ```
 _____   _____   __  __                          _ _ 
 |  __ \ / ____| |  \/  |                        | | |
@@ -67,6 +67,8 @@ solving the wave equation in a rectangular domain in the
 - [arrayfire](http://arrayfire.org)
 - [gmshtranslate](https://github.com/amanabt/gmshtranslator)
 - [texlive](https://www.tug.org/texlive/)
+- [coveralls](https://pypi.python.org/pypi/coveralls)
+- [python-coveralls](https://pypi.python.org/pypi/python-coveralls/)
 
 ## Documentation
 The documenation for the project can be generated using `sphinx`
@@ -79,12 +81,12 @@ The unit tests are located in the
 
 ### Dependencies
 - [pytest](https://docs.pytest.org/en/latest/#)
-
+- [pytest-cov](https://pypi.python.org/pypi/pytest-cov)
 ### Running Unit Tests
 To run the unit tests, enter the following commands
 ```
 $ cd path/to/the/DG_Maxwell/repo
-$ pytest
+$ pytest --verbose -r P --color=yes --cov dg_maxwell
 ```
 * The parameters of the simulation are stored in params.py in
   the app folder, These can be changed accordingly.

dg_maxwell/__init__.py

@@ -2,6 +2,7 @@
 # -*- coding: utf-8 -*-
 
 import arrayfire as af
-af.set_backend('cpu')
 
 from dg_maxwell import params
+
+af.set_backend(params.backend)

dg_maxwell/isoparam.py

@@ -2,8 +2,11 @@
 # -*- coding: utf-8 -*-
 
 import arrayfire as af
-af.set_backend('cpu')
+
 from dg_maxwell import utils
+from dg_maxwell import params
+
+af.set_backend(params.backend)
 
 def isoparam_1D(x_nodes, xi):
     '''
@@ -12,27 +15,27 @@ def isoparam_1D(x_nodes, xi):
 
     Parameters
     ----------
-    
+
     x_nodes : arrayfire.Array [2 1 1 1]
               Element nodes.
-    
+
     xi      : arrayfire.Array [N 1 1 1]
               Value of :math:`\\xi` coordinate for which the corresponding
               :math:`x` coordinate is to be found.
-    
+
     Returns
     -------
     x : arrayfire.Array
         :math:`x` value in the element corresponding to :math:`\\xi`.
     '''
     N_0 = (1 - xi) / 2
     N_1 = (1 + xi) / 2
-    
+
     N0_x0 = af.broadcast(utils.multiply, N_0, x_nodes[0])
     N1_x1 = af.broadcast(utils.multiply, N_1, x_nodes[1])
-    
+
     x = N0_x0 + N1_x1
-    
+
     return x
 
 def isoparam_x_2D(x_nodes, xi, eta):
@@ -43,7 +46,7 @@ def isoparam_x_2D(x_nodes, xi, eta):
 
     Here :math:`P_i` corresponds to :math:`(\\xi_i, \\eta_i)` coordinates,
     :math:`i \in \\{0, 1, ..., 7\\}` respectively, where,
-    
+
     .. math:: (\\xi_0, \\eta_0) &\equiv (-1,  1) \\\\
               (\\xi_1, \\eta_1) &\equiv (-1,  0) \\\\
               (\\xi_2, \\eta_2) &\equiv (-1, -1) \\\\
@@ -52,12 +55,12 @@ def isoparam_x_2D(x_nodes, xi, eta):
               (\\xi_5, \\eta_5) &\equiv ( 1,  0) \\\\
               (\\xi_6, \\eta_6) &\equiv ( 1,  1) \\\\
               (\\xi_7, \\eta_7) &\equiv ( 0,  1)
-              
+
     Parameters
     ----------
     x_nodes : np.ndarray [8]
               :math:`x` nodes.
-              
+
     xi      : float
             :math:`\\xi` coordinate for which :math:`x` has to be found.
 
@@ -68,7 +71,6 @@ def isoparam_x_2D(x_nodes, xi, eta):
     -------
     x : float
         :math:`x` coordinate corresponding to :math:`(\\xi, \\eta)` coordinate.
-
     '''
     N_0 = (-1.0 / 4.0) * (1 - xi)  * (1 + eta) * (1 + xi - eta)
     N_1 = (1.0 / 2.0)  * (1 - xi)  * (1 - eta**2)
@@ -95,12 +97,12 @@ def isoparam_y_2D(y_nodes, xi, eta):
     '''
     This function allows isoparametric mapping of a :math:`2^{nd}` order
     element with :math:`8` nodes
-    
+
     .. math:: (P_0, P_1, P_2, P_3, P_4, P_5, P_6, P_7)
-    
+
     Here :math:`P_i` corresponds to :math:`(\\xi_i, \\eta_i)` coordinates,
     :math:`i \in \\{0, 1, ..., 7\\}` respectively, where,
-    
+
     .. math:: (\\xi_0, \\eta_0) &\equiv (-1,  1) \\\\
               (\\xi_1, \\eta_1) &\equiv (-1,  0) \\\\
               (\\xi_2, \\eta_2) &\equiv (-1, -1) \\\\
@@ -109,12 +111,12 @@ def isoparam_y_2D(y_nodes, xi, eta):
               (\\xi_5, \\eta_5) &\equiv ( 1,  0) \\\\
               (\\xi_6, \\eta_6) &\equiv ( 1,  1) \\\\
               (\\xi_7, \\eta_7) &\equiv ( 0,  1)
-              
+
     Parameters
     ----------
     y_nodes : np.ndarray [8]
               :math:`y` nodes.
-              
+
     xi      : float
             :math:`\\xi` coordinate for which :math:`y` has to be found.
 
@@ -125,6 +127,5 @@ def isoparam_y_2D(y_nodes, xi, eta):
     -------
     float
         :math:`y` coordinate corresponding to :math:`(\\xi, \\eta)` coordinate.
-
     '''
     return isoparam_x_2D(y_nodes, xi, eta)