diff --git a/Geo/SVector3.h b/Geo/SVector3.h index 96e796a3298e645cca237ba4aecf16f7edab2b3a..97683d7b094c6921f8eda20df7ffaa7cfb609001 100644 --- a/Geo/SVector3.h +++ b/Geo/SVector3.h @@ -91,6 +91,11 @@ class SVector3 { } const double* data() const {return P.data();} double* data() {return P.data();} + double axpy(const double a, const SVector3& y){ + for (int i=0; i<3; i++){ + P[i] += (a*y[i]); + } + } }; inline double dot(const SVector3 &a, const SVector3 &b)