Revision 32, 0.8 kB
(checked in by smidl, 17 years ago)
|
test KF : estimation of R in KF is not possible! Likelihood of y_t is growing when R -> 0
|
Rev | Line | |
---|
[32] | 1 | # This file defines handy gdb macros for printing out ITPP types |
---|
| 2 | # To use it, add this line to your ~/.gdbinit : |
---|
| 3 | # source gdb_bdm |
---|
| 4 | |
---|
| 5 | define dm |
---|
| 6 | set $i=0 |
---|
| 7 | |
---|
| 8 | set $M = $arg0 |
---|
| 9 | set $rs = $M.no_rows |
---|
| 10 | set $col = $M.no_cols |
---|
| 11 | |
---|
| 12 | while $i < $rs |
---|
| 13 | set $j=0 |
---|
| 14 | printf "[" |
---|
| 15 | while $j< $col |
---|
| 16 | printf "%f,", $M.data[$j*$col+$i] |
---|
| 17 | set $j++ |
---|
| 18 | end |
---|
| 19 | printf "]" |
---|
| 20 | printf "\n" |
---|
| 21 | set $i++ |
---|
| 22 | end |
---|
| 23 | end |
---|
| 24 | |
---|
| 25 | define dme |
---|
| 26 | set $i=0 |
---|
| 27 | |
---|
| 28 | set $M = $arg0 |
---|
| 29 | set $rs = $M.no_rows |
---|
| 30 | set $col = $M.no_cols |
---|
| 31 | |
---|
| 32 | while $i < $rs |
---|
| 33 | set $j=0 |
---|
| 34 | printf "[" |
---|
| 35 | while $j< $col |
---|
| 36 | printf "%e,", $M.data[$j*$col+$i] |
---|
| 37 | set $j++ |
---|
| 38 | end |
---|
| 39 | printf "]" |
---|
| 40 | printf "\n" |
---|
| 41 | set $i++ |
---|
| 42 | end |
---|
| 43 | end |
---|
| 44 | |
---|
| 45 | define dv |
---|
| 46 | print $arg0.data[0]@($arg0.datasize) |
---|
| 47 | end |
---|
| 48 | |
---|
| 49 | define dl |
---|
| 50 | dm $M.L |
---|
| 51 | dv $M.D |
---|
| 52 | end |
---|