Changeset e85cfd for src/solver/givens.hpp

Timestamp:

Apr 9, 2013, 9:45:13 AM (13 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

0bd47e

Parents:

4a709e

Message:

Work.

File:

• src/solver/givens.hpp (modified) (2 diffs)

src/solver/givens.hpp

@@ -42 +42 @@
 public:
   Givens(bool register_ = true) :
-    T(register_)
-  {}
+  T(register_)
+{}
 
   Givens(int size, bool register_ = true) :
-    T(size, register_)
+  T(size, register_)
   {}
 
@@ -63 +63 @@
       if (fabs(this->Mat(j,i)) > DBL_EPSILON) {
 
-        t = 1.0 / sqrt(this->Mat(i,i)*this->Mat(i,i) + this->Mat(j,i)*this->Mat(j,i));
-        s = t * this->Mat(j,i);
-        c = t * this->Mat(i,i);
+        t = 1.0 / sqrt(this->Mat(i,i)*this->Mat(i,i) + this->Mat(j,i)*this->Mat(j,i));
+        s = t * this->Mat(j,i);
+        c = t * this->Mat(i,i);
 
-        for (int k=i; k<n; k++) {
+        for (int k=i; k<n; k++) {
 
-          t = c * this->Mat(i,k) + s * this->Mat(j,k);
+          t = c * this->Mat(i,k) + s * this->Mat(j,k);
 
-          if (k != i)
-            this->Mat(j,k) = c * this->Mat(j,k) - s * this->Mat(i,k);
+          if (k != i)
+            this->Mat(j,k) = c * this->Mat(j,k) - s * this->Mat(i,k);
 
-          this->Mat(i,k) = t;
+          this->Mat(i,k) = t;
 
-        }
+        }
 
-        t = c * this->Rhs(i) + s * this->Rhs(j);
+        t = c * this->Rhs(i) + s * this->Rhs(j);
 
-        this->Rhs(j) = c * this->Rhs(j) - s * this->Rhs(i);
+        this->Rhs(j) = c * this->Rhs(j) - s * this->Rhs(i);
 
-        this->Rhs(i) = t;
+        this->Rhs(i) = t;
 
-        this->Mat(j,i) = 0.0;
+        this->Mat(j,i) = 0.0;
 
       }