From fee0a6f8e7c77a0d22c6f5e798d52d9a2bd85871 Mon Sep 17 00:00:00 2001
From: Guillaume Bittoun <guillaume.bittoun@gmail.com>
Date: Tue, 18 Nov 2014 21:30:03 +0000
Subject: [PATCH] simplex_interpolator: bug fix

---
 src/lol/math/simplex_interpolator.h | 9 ++-------
 1 file changed, 2 insertions(+), 7 deletions(-)

diff --git a/src/lol/math/simplex_interpolator.h b/src/lol/math/simplex_interpolator.h
index 02648229..2b280791 100644
--- a/src/lol/math/simplex_interpolator.h
+++ b/src/lol/math/simplex_interpolator.h
@@ -128,16 +128,12 @@ protected:
 
     inline void InitBase()
     {
-        this->samples.SetSize(vec_t<int, 2>(N, N));
-        this->samples.SetSize(vec_t<int, 2>(N, N));
-
         for (int i = 0 ; i < N ; ++i)
         {
             for (int j = i ; j < N ; ++j)
             {
-                this->base[i][j] = sqrt((i+2)/((float)(2*i+2))) / (j > i ? i+2 : 1);
-                this->base_inverse[i][j] = sqrt((2*j+2) / ((float)(j+2))) * (j > i ? (1 / (float)(j+1)) : 1);
-                this->diagonal[i] += (this->base[i][j]);
+                this->base[j][i] = sqrt((i+2)/((float)(2*i+2))) / (j > i ? i+2 : 1);
+                this->base_inverse[j][i] = sqrt((2*j+2) / ((float)(j+2))) * (j > i ? (-1 / (float)(j+1)) : 1);
             }
         }
     }
@@ -147,7 +143,6 @@ protected:
 
     arraynd<N, T> samples;
 
-    vec_t<T, N> diagonal;
 };
 
 }