From 46aba98d4c0f625623b7047395db2332126a4948 Mon Sep 17 00:00:00 2001 From: Francis Russell Date: Sat, 19 May 2012 17:44:15 +0100 Subject: [PATCH] Fix reciprocal vector calculation in Laplacian. --- src/ofc/generators/onetep/Laplacian.scala | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) 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) } -- 2.47.3