fix bugs in Gradient and test

This commit is contained in:
soypat
2022-08-04 11:23:44 -03:00
committed by Dan Kortschak
parent 1fbf8fb724
commit 177c87adcd
2 changed files with 6 additions and 6 deletions

View File

@@ -104,9 +104,9 @@ func Gradient(p, step Vec, field func(Vec) float64) Vec {
dy := Vec{Y: step.Y}
dz := Vec{Z: step.Z}
return Vec{
X: field(Add(p, dx)) - field(Sub(p, dx)),
Y: field(Add(p, dy)) - field(Sub(p, dy)),
Z: field(Add(p, dz)) - field(Sub(p, dz)),
X: (field(Add(p, dx)) - field(Sub(p, dx))) / (2 * step.X),
Y: (field(Add(p, dy)) - field(Sub(p, dy))) / (2 * step.Y),
Z: (field(Add(p, dz)) - field(Sub(p, dz))) / (2 * step.Z),
}
}