| 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 ); |