From 67acac80d920767585b5581a1d0a09572787a351 Mon Sep 17 00:00:00 2001 From: Gertjan van den Burg Date: Mon, 7 Nov 2016 14:34:23 +0100 Subject: compute ZAZ in blocks to increase numerical precision --- include/gensvm_base.h | 2 ++ 1 file changed, 2 insertions(+) (limited to 'include') diff --git a/include/gensvm_base.h b/include/gensvm_base.h index e4c2be1..912c477 100644 --- a/include/gensvm_base.h +++ b/include/gensvm_base.h @@ -140,6 +140,8 @@ struct GenWork { ///< (K-1) x (m+1) working matrix for the Z'*B calculation double *ZAZ; ///< (m+1) x (m+1) working matrix for the Z'*A*Z calculation + double *tmpZAZ; + ///< (m+1) x (m+1) temporary working matrix for the Z'*A*Z calculation double *ZV; ///< n x (K-1) working matrix for the Z * V calculation double *beta; -- cgit v1.2.3