f(x,y,z) = 1.0 / sqrt(x*x + y*y + z*z);
(charges)
| 21 | return parser.parse_args(); |
| 22 | |
| 23 | def test_function(charges): |
| 24 | """ |
| 25 | f(x,y,z) = 1.0 / sqrt(x*x + y*y + z*z); |
| 26 | """ |
| 27 | return lambda pts: np.sum([1.0 / numpy.linalg.norm(pts - c, axis=1) |
| 28 | for c in charges], axis=0); |
| 29 | |
| 30 | def test_function_grad(charges): |
| 31 | """ |