![Saree90](/assets/img/avatar_default.png)
Mat_<double> coord_sum(1,ndim,0.0),buf(1,ndim,0.0),y(1,ndim,0.0); nfunk = 0; for(i=0;i<ndim+1;++i) { y(i) = f->calc(p[i]); } y has only ndim elements, while the loop goes over ndim+1
Mat_<double> coord_sum(1,ndim,0.0),buf(1,ndim,0.0),y(1,ndim,0.0); nfunk = 0; for(i=0;i<ndim+1;++i) { y(i) = f->calc(p[i]); } y has only ndim elements, while the loop goes over ndim+1