Opened 14 years ago
#1401 new defect
gsl messages
| Reported by: | Jorge | Owned by: | Víctor de Buen Remiro |
|---|---|---|---|
| Priority: | normal | Milestone: | Mantainance |
| Component: | Kernel | Version: | head |
| Severity: | critical | Keywords: | |
| Cc: |
Description
The messages emited by gsl are not shown in the log, try the following code:
Real _.m = 10;
Matrix Xb = Rand(_.m,1,-120,120);
Matrix z_ = Xb + Gaussian(_.m,1,0,1);
Matrix p_ = F01(z_);
Matrix _.Y = GE(Rand(_.m,1,0,1),0.5); //GE(p_,0.5);
Matrix _inf = Constant(_.m,1,1/0);
Matrix _zero = Constant(_.m,1,0);
Matrix _one = Constant(1,1,1);
Matrix a = IfMat(_.Y, - Xb, -_inf);
//WriteLn("TRACE [@Filter.Probit::eval.Filter.Probit] 3");
Matrix b = IfMat(_.Y, +_inf, - Xb);
Matrix e_ = RandTruncatedMultNormal(_zero, _one, a, b);
the output of gsl is:
WARNING in bmath/gsl_ext/rtnorm.c, gsl_rtbnorm_combo, 157: Mean extremely far from truncation point. Returning truncation point WARNING in bmath/gsl_ext/rtnorm.c, gsl_rtanorm_combo, 201: Mean extremely far from truncation point. Returning truncation point
Note: See
TracTickets for help on using
tickets.
