| 51 | } |
| 52 | |
| 53 | int actual_df(const VectorXd &x, MatrixXd &fjac) const |
| 54 | { |
| 55 | double tmp1, tmp2, tmp3, tmp4; |
| 56 | for (int i = 0; i < values(); i++) |
| 57 | { |
| 58 | tmp1 = i+1; |
| 59 | tmp2 = 16 - i - 1; |
| 60 | tmp3 = (i>=8)? tmp2 : tmp1; |
| 61 | tmp4 = (x[1]*tmp2 + x[2]*tmp3); tmp4 = tmp4*tmp4; |
| 62 | fjac(i,0) = -1; |
| 63 | fjac(i,1) = tmp1*tmp2/tmp4; |
| 64 | fjac(i,2) = tmp1*tmp3/tmp4; |
| 65 | } |
| 66 | return 0; |
| 67 | } |
| 68 | }; |
| 69 | |
| 70 | void test_forward() |
no outgoing calls
no test coverage detected