Remove z1.

memmett · memmett · commit 3f315b382ed1 · 2015-04-08T07:57:56.000-06:00
diff --git a/src/pfasst/quadrature/polynomial_impl.hpp b/src/pfasst/quadrature/polynomial_impl.hpp
@@ -73,10 +73,9 @@ namespace pfasst
       size_t n = c.size() - 1;
 
       // initial guess
-      vector<complex<CoeffT>> z0(n), z1(n);
+      vector<complex<CoeffT>> z0(n);
       for (size_t j = 0; j < n; j++) {
         z0[j] = pow(complex<double>(0.4, 0.9), j);
-        z1[j] = z0[j];
       }
 
       // durand-kerner-weierstrass iterations
@@ -92,7 +91,6 @@ namespace pfasst
           }
           z0[i] = z0[i] - num / den;
         }
-        z1 = z0;
       }
 
       vector<CoeffT> roots(n);

Original file line number	Diff line number	Diff line change
`@@ -73,10 +73,9 @@ namespace pfasst`
`73`	`73`	`size_t n = c.size() - 1;`
`74`	`74`
`75`	`75`	`// initial guess`
`76`		`- vector<complex<CoeffT>> z0(n), z1(n);`
	`76`	`+ vector<complex<CoeffT>> z0(n);`
`77`	`77`	`for (size_t j = 0; j < n; j++) {`
`78`	`78`	`z0[j] = pow(complex<double>(0.4, 0.9), j);`
`79`		`- z1[j] = z0[j];`
`80`	`79`	`}`
`81`	`80`
`82`	`81`	`// durand-kerner-weierstrass iterations`
`@@ -92,7 +91,6 @@ namespace pfasst`
`92`	`91`	`}`
`93`	`92`	`z0[i] = z0[i] - num / den;`
`94`	`93`	`}`
`95`		`- z1 = z0;`
`96`	`94`	`}`
`97`	`95`
`98`	`96`	`vector<CoeffT> roots(n);`