Commit 5f9cbdad authored by Oleg Dzhimiev's avatar Oleg Dzhimiev

testing reiteration

parent 68f03f5f
...@@ -215,6 +215,9 @@ numbers.calculus.GaussNewton_forHeading = function(v,n,r,dr,eps,w){ ...@@ -215,6 +215,9 @@ numbers.calculus.GaussNewton_forHeading = function(v,n,r,dr,eps,w){
}; };
} }
console.log("V0:")
console.log(v0)
while(!stop){ while(!stop){
counter++ counter++
...@@ -233,6 +236,9 @@ numbers.calculus.GaussNewton_forHeading = function(v,n,r,dr,eps,w){ ...@@ -233,6 +236,9 @@ numbers.calculus.GaussNewton_forHeading = function(v,n,r,dr,eps,w){
diff.push(r(i,v1)-r(i,v0)) diff.push(r(i,v1)-r(i,v0))
} }
console.log("V1:")
console.log(v1)
console.log("residuals old:") console.log("residuals old:")
console.log(rs0) console.log(rs0)
console.log("residuals new:") console.log("residuals new:")
...@@ -255,9 +261,16 @@ numbers.calculus.GaussNewton_forHeading = function(v,n,r,dr,eps,w){ ...@@ -255,9 +261,16 @@ numbers.calculus.GaussNewton_forHeading = function(v,n,r,dr,eps,w){
} }
if (sscale!=0){ if (sscale!=0){
console.log("SCALE = "+sscale)
for(var i=0;i<n;i++){ for(var i=0;i<n;i++){
v1[i] = v0[i] + diff[i]*sscale v1[i] = v0[i] + diff[i]*sscale
} }
console.log("Corrected V1:")
console.log(v1)
} }
var s0 = sigma(v0,n,r) var s0 = sigma(v0,n,r)
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment