26 use zero::ZeroFn; |
26 use zero::ZeroFn; |
27 use alg_tools::mapping::Sum; |
27 use alg_tools::mapping::Sum; |
28 use alg_tools::iterate::{AlgIteratorOptions, Verbose}; |
28 use alg_tools::iterate::{AlgIteratorOptions, Verbose}; |
29 |
29 |
30 let points = [ |
30 let points = [ |
|
31 //OnCube::new(F1, Loc([0.5, 0.5])), |
|
32 //OnCube::new(F2, Loc([0.5, 0.5])), |
|
33 //OnCube::new(F4, Loc([0.1, 0.1])), |
31 OnCube::new(F1, Loc([0.5, 0.5])), |
34 OnCube::new(F1, Loc([0.5, 0.5])), |
|
35 OnCube::new(F3, Loc([0.5, 0.5])), |
32 OnCube::new(F2, Loc([0.5, 0.5])), |
36 OnCube::new(F2, Loc([0.5, 0.5])), |
33 OnCube::new(F4, Loc([0.1, 0.1])), |
|
34 ]; |
37 ]; |
35 |
38 |
36 //let x = points[0].clone(); |
39 //let x = points[0].clone(); |
37 let x = OnCube::new(F6, Loc([0.5, 0.5])); |
40 // OnCube::new(F3, Loc([0.5, 0.5])); goes to opposite side |
|
41 let x = OnCube::new(F3, Loc([0.5, 0.4])); |
38 let f = Sum::new(points.into_iter().map(DistToSquaredDiv2)); |
42 let f = Sum::new(points.into_iter().map(DistToSquaredDiv2)); |
39 let g = ZeroFn::new(); |
43 let g = ZeroFn::new(); |
40 let τ = 0.1; |
44 let τ = 0.1; |
41 let iter = AlgIteratorOptions{ |
45 let iter = AlgIteratorOptions{ |
42 max_iter : 100, |
46 max_iter : 100, |