[Toon-members] TooN Lapack_Cholesky.h

[Edward Rosten]

CVSROOT:        /cvsroot/toon
Module name:    TooN
Changes by:     Edward Rosten <edrosten>        09/10/19 14:15:57

Modified files:
        .              : Lapack_Cholesky.h 

Log message:
        Single precision LAPACK Cholesky.

CVSWeb URLs:
http://cvs.savannah.gnu.org/viewcvs/TooN/Lapack_Cholesky.h?cvsroot=toon&r1=1.3&r2=1.4

Patches:
Index: Lapack_Cholesky.h
===================================================================
RCS file: /cvsroot/toon/TooN/Lapack_Cholesky.h,v
retrieving revision 1.3
retrieving revision 1.4
diff -u -b -r1.3 -r1.4
--- Lapack_Cholesky.h   19 Oct 2009 14:08:52 -0000      1.3
+++ Lapack_Cholesky.h   19 Oct 2009 14:15:56 -0000      1.4
@@ -100,7 +100,7 @@
        void do_compute(){
                int N = my_cholesky.num_rows();
                int info;
-               dpotrf_("L", &N, (double*) my_cholesky_lapack.my_data, &N, 
&info);
+               potrf_("L", &N, my_cholesky_lapack.my_data, &N, &info);
                for (int i=0;i<N;i++) {
                  int j;
                  for (j=0;j<=i;j++) {
@@ -128,7 +128,7 @@
                int N=my_cholesky.num_rows();
                int NRHS=1;
                int info;
-               dpotrs_("L", &N, &NRHS, my_cholesky_lapack.my_data, &N, 
result.my_data, &N, &info);     
+               potrs_("L", &N, &NRHS, my_cholesky_lapack.my_data, &N, 
result.my_data, &N, &info);     
                assert(info==0);
                return result;
        }
@@ -141,7 +141,7 @@
                int N=my_cholesky.num_rows();
                int NRHS=m.num_cols();
                int info;
-               dpotrs_("L", &N, &NRHS, my_cholesky_lapack.my_data, &N, 
result.my_data, &N, &info);     
+               potrs_("L", &N, &NRHS, my_cholesky_lapack.my_data, &N, 
result.my_data, &N, &info);     
                assert(info==0);
                return result;
        }
@@ -167,7 +167,7 @@
                M=my_cholesky_lapack;
                int N = my_cholesky.num_rows();
                int info;
-               dpotri_("L", &N, M.my_data, &N, &info);
+               potri_("L", &N, M.my_data, &N, &info);
                assert(info == 0);
                for (int i=1;i<N;i++) {
                  for (int j=0;j<i;j++) {
@@ -178,7 +178,7 @@
        }
 
 private:
-       Matrix<Size,Size,double> my_cholesky_lapack;     
+       Matrix<Size,Size,Precision> my_cholesky_lapack;     
        Matrix<Size,Size,Precision> my_cholesky;     
        int my_rank;
 };