From dda165d0288ca4cd585355917c84c29ec1731bf0 Mon Sep 17 00:00:00 2001
From: Mark Jenkinson <mark@fmrib.ox.ac.uk>
Date: Mon, 13 Nov 2000 15:38:00 +0000
Subject: [PATCH] changed getdiag() to diag()
---
glmrand.cc | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/glmrand.cc b/glmrand.cc
index 9c2fb91..6ab42f2 100644
--- a/glmrand.cc
+++ b/glmrand.cc
@@ -185,13 +185,15 @@ namespace SIGPROC {
// e is the estimate of the variance of the timeseries, sigma^2
Matrix varmat = (r.Columns(batch_pos, numrand+1).t()*r.Columns(batch_pos, numrand+1))*var_on_e/sizeTS;
ColumnVector vartemp;
- getdiag(vartemp, varmat);
+ //getdiag(vartemp, varmat); // obsolete fn
+ vartemp = diag(varmat); // MJ NOTE: new fn
var.Rows(batch_pos, numrand+1) = vartemp;
}
else
{
varmatfull = (r.Columns(batch_pos, batch_pos+batch_size-1).t()*r.Columns(batch_pos, batch_pos+batch_size-1))*var_on_e/sizeTS;
- getdiag(vartempfull, varmatfull);
+ //getdiag(vartempfull, varmatfull); // obsolete fn
+ vartempfull = diag(varmatfull); // MJ NOTE: new fn
var.Rows(batch_pos, batch_pos+batch_size-1) = vartempfull;
}
batch_pos += batch_size;
--
GitLab