From de449709e34c619be468356ea56d110fc94ef1e3 Mon Sep 17 00:00:00 2001
From: Jukka Aho <jukka.aho@kapsi.fi>
Date: Mon, 14 Sep 2015 18:50:17 +0300
Subject: [PATCH] hajoitus

---
 src/equations.jl    | 9 +++++++++
 src/hierarchical.jl | 2 +-
 2 files changed, 10 insertions(+), 1 deletion(-)

diff --git a/src/equations.jl b/src/equations.jl
index ba07347..22232f3 100644
--- a/src/equations.jl
+++ b/src/equations.jl
@@ -50,6 +50,15 @@ function get_detJ(el::Element, ip::IntegrationPoint)
         return det(J)
     end
 end
+function get_detJ(el::Element, xi::Vector)
+    J = get_Jacobian(el, xi)
+    n, m = size(J)
+    if n != m # for manifolds
+        return norm(J)
+    else
+        return det(J)
+    end
+end
 
 """
 Integrate f over element
diff --git a/src/hierarchical.jl b/src/hierarchical.jl
index 98bd4a0..f0775f3 100644
--- a/src/hierarchical.jl
+++ b/src/hierarchical.jl
@@ -29,7 +29,7 @@ end
 Return derivative of Legendgre polynomial of order n to inverval ξ ∈  [-1, 1]
 """
 function get_legendre_polynomial_derivative(n::Int)
-    dP(xi) = 2^n*sum([k*xi.^(k-1)*bin(n, k)*bin(1/2*(n+k-1), n) for k=0:n])
+    dP(xi) = 2^n*sum([k*xi.^(k-1)*bin(n, k)*bin(1/2*(n+k-1), n) for k=1:n])
     dP
 end