I don't know the code well, so I am not 100% sure this is correct, but glancing around I see += / -= used for these *temp objects, so I think it is (the alternative is ntemp = -wt[p1];, which is what the current code is doing, but with clearer spacing).
As flagged by our compiler.
I don't know the code well, so I am not 100% sure this is correct, but glancing around I see
+=
/-=
used for these*temp
objects, so I think it is (the alternative isntemp = -wt[p1];
, which is what the current code is doing, but with clearer spacing).