| 126 | } |
| 127 | |
| 128 | VectorXf computeFunction( EnergyFunction & efun, const VectorXf & x, const VectorXf & dx, int n_samples ) { |
| 129 | VectorXf r( n_samples ); |
| 130 | VectorXf tmp = x; |
| 131 | for( int i=0; i<n_samples; i++ ) |
| 132 | r[i] = efun.gradient( x+i*dx, tmp ); |
| 133 | return r; |
| 134 | } |
| 135 |