summaryrefslogtreecommitdiff
path: root/lib
diff options
context:
space:
mode:
Diffstat (limited to 'lib')
-rw-r--r--lib/orthogonal.h2
-rw-r--r--lib/vector.h17
2 files changed, 2 insertions, 17 deletions
diff --git a/lib/orthogonal.h b/lib/orthogonal.h
index 34cd44e..85d0e11 100644
--- a/lib/orthogonal.h
+++ b/lib/orthogonal.h
@@ -24,6 +24,7 @@ class orthogonal_t : public std::array<std::array<T, q>, q> {
vector_t<q, T> act(const vector_t <q, T>& v) const {
vector_t <q, T> v_rot;
+ v_rot.fill(0);
if (is_reflection) {
double prod = 0;
@@ -80,6 +81,7 @@ class orthogonal_t : public std::array<std::array<T, q>, q> {
return this->act(v); // reflections are their own inverse
} else {
vector_t <q, T> v_rot;
+ v_rot.fill(0);
for (q_t i = 0; i < q; i++) {
for (q_t j = 0; j < q; j++) {
diff --git a/lib/vector.h b/lib/vector.h
index 2f4077a..72633a8 100644
--- a/lib/vector.h
+++ b/lib/vector.h
@@ -7,23 +7,6 @@
#include "types.h"
-/* The following is the minimum definition of a spin class.
- *
- * The class must contain an M_t and an F_t for holding the sum of an
- * integer number of spins and a double-weighted number of spins,
- * respectively.
- *
- * void init(X_t *p);
- * void free_spin(X_t p);
- * X_t copy(X_t x);
- * void add(M_t *x1, int factor, X_t x2);
- * void add(F_t *x1, double factor, X_t x2);
- * M_t scalar_multiple(int factor, X_t x);
- * double norm_squared(F_t x);
- * void write_magnetization(M_t M, FILE *outfile);
- *
- */
-
template <q_t q, class T>
class vector_t : public std::array<T, q> {
public: