From: Francis Russell <francis@unchartedbackwaters.co.uk>
Date: Sat, 19 May 2012 16:44:15 +0000 (+0100)
Subject: Fix reciprocal vector calculation in Laplacian.
X-Git-Url: https://git.unchartedbackwaters.co.uk/w/?a=commitdiff_plain;h=46aba98d4c0f625623b7047395db2332126a4948;p=francis%2Fofc.git

Fix reciprocal vector calculation in Laplacian.
---

diff --git a/src/ofc/generators/onetep/Laplacian.scala b/src/ofc/generators/onetep/Laplacian.scala
index 445337e..ff3ef49 100644
--- a/src/ofc/generators/onetep/Laplacian.scala
+++ b/src/ofc/generators/onetep/Laplacian.scala
@@ -35,7 +35,7 @@ class Laplacian(op: Field)  extends Field {
         var component : Expression[FloatType] = new FloatLiteral(0.0)
         for(vec <- 0 to 2) {
           val vector = OnetepTypes.CellInfo.public % OnetepTypes.CellInfo.latticeReciprocal(vec)
-          component = component + vector % OnetepTypes.Point.coord(dim) * new Conversion[IntType, FloatType](indices(dim))
+          component = component + vector % OnetepTypes.Point.coord(dim) * new Conversion[IntType, FloatType](indices(vec) - 1)
         }
         loops(0) += new AssignStatement(reciprocalVector(dim), component)
       }