89 | | CHECK_CLOSE_EX ( hepdf->evallog ( zeron ), mEp.evallogcond ( zeron, zero ), tolerance ); |
| 89 | double lpz = hepdf->evallog ( zeron ); |
| 90 | CHECK_CLOSE_EX ( lpz, mEp.evallogcond ( zeron, zero ), tolerance ); |
| 91 | |
| 92 | vec lpzv(1); |
| 93 | lpzv(0) = lpz; |
| 94 | |
| 95 | mat zero1n ( hepdf->dimension(), 1 ); |
| 96 | for ( int i = 0; i < zero1n.rows(); ++i ) { |
| 97 | zero1n ( i, 0 ) = 0; |
| 98 | } |
| 99 | |
| 100 | vec lpzv_act = hepdf->evallog_m ( zero1n ); |
| 101 | CHECK_CLOSE_EX ( lpzv, lpzv_act, tolerance ); |
| 102 | |
| 103 | Array<vec> zeroa(3); |
| 104 | lpzv = vec( zeroa.size() ); |
| 105 | for ( int i = 0; i < zeroa.size(); ++i ) { |
| 106 | zeroa(i) = zeron; |
| 107 | lpzv(i) = lpz; |
| 108 | } |
| 109 | |
| 110 | lpzv_act = hepdf->evallog_m ( zeroa ); |
| 111 | CHECK_CLOSE_EX ( lpzv, lpzv_act, tolerance ); |