Commit: 2e2bb6e1c5d6d899c15ec5ebd68312bdcb1acdd1 - MolDS (git)

Commit MetaInfo

• Révision: 2e2bb6e1c5d6d899c15ec5ebd68312bdcb1acdd1 (tree)
• l'heure: 2014-01-17 17:17:18
• Auteur: Katsuhiko Nishimra <ktns.87@gmai...>
• Commiter: Katsuhiko Nishimra

Message de Log

Add wrapper functions of BLAS dnrm2 routine. #30379

git-svn-id: https://svn.sourceforge.jp/svnroot/molds/trunk@1653 1136aad2-a195-0410-b898-f5ea1d11b9d8

Change Summary

• modified: src/wrappers/Blas.cpp
• modified: src/wrappers/Blas.h

Modification

--- a/src/wrappers/Blas.cpp
+++ b/src/wrappers/Blas.cpp
@@ -132,6 +132,25 @@ double Blas::Ddot(molds_blas_int n,
    return cblas_ddot(n, x, incrementX, y, incrementY);
 }
 
+// returns sqrt(sum(Xi^2))
+// vectorX: n-vector
+double Blas::Dnrm2(molds_blas_int n,
+                   double const* vectorX) const{
+   molds_blas_int incrementX = 1;
+   return this->Dnrm2(n,vectorX,incrementX);
+}
+// returns sqrt(sum(Xi^2))
+// vectorX: n-vector
+double Blas::Dnrm2(molds_blas_int n,
+                   double const* vectorX,
+                   molds_blas_int incrementX) const{
+   if(n<=0 || vectorX == NULL || incrementX <= 0){
+      return 0.0;
+   }
+   double* x=const_cast<double*>(vectorX);
+   return cblas_dnrm2(n,x,incrementX);
+}
+
 // returns max(abs(vectorX[i]))
 //    vectorX: n-vector
 double Blas::Damax(molds_blas_int n,

--- a/src/wrappers/Blas.h
+++ b/src/wrappers/Blas.h
@@ -45,6 +45,11 @@ public:
    double Ddot(molds_blas_int n,
                double const* vectorX, molds_blas_int incrementX,
                double const* vectorY, molds_blas_int incrementY)const;
+   double Dnrm2(molds_blas_int n,
+                double const* vectorX) const;
+   double Dnrm2(molds_blas_int n,
+                double const* vectorX,
+                molds_blas_int incrementX) const;
    double Damax(molds_blas_int n,
                       double const* vectorX) const;
    double Damax(molds_blas_int n,