| 109 | } |
| 110 | |
| 111 | Bayes_base::Float |
| 112 | Information_scheme::predict (Linrz_predict_model& f) |
| 113 | /* Extended linrz information prediction |
| 114 | * Computation is through state to accommodate linearised model |
| 115 | */ |
| 116 | { |
| 117 | update (); // x,X required |
| 118 | x = f.f(x); // Extended Kalman state predict is f(x) directly |
| 119 | // Predict information matrix, and state covariance |
| 120 | noalias(X) = prod_SPD(f.Fx,X, tempX); |
| 121 | noalias(X) += prod_SPD(f.G, f.q, tempX); |
| 122 | |
| 123 | // Information |
| 124 | Float rcond = UdUinversePD (Y, X); |
| 125 | rclimit.check_PD(rcond, "X not PD in predict"); |
| 126 | // Predict information state |
| 127 | noalias(y) = prod(Y,x); |
| 128 | return rcond; |
| 129 | } |
| 130 | |
| 131 | Float Information_scheme::predict (Linear_invertable_predict_model& f, Linear_predict_byproducts& b) |
| 132 | /* Linear information predict |
nothing calls this directly
no test coverage detected