Changeset 2b448e9 in tol

Timestamp:

Mar 9, 2009, 8:57:02 AM (18 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

1941-docker-image-to-test, devel, master

Children:

516e91c

Parents:

d11b5b0

Message:

New method Drop

git-svn-id: ​https://www.tol-project.org/svn/tolp/trunk@704 be4846ae-f051-483a-831f-b64e7d08527f

File:

• tol/btol/vmatrix_type/vmat_arith.cpp (modified) (2 diffs)

tol/btol/vmatrix_type/vmat_arith.cpp

@@ -67 +67 @@
 DefineMonary (ACosH,      fRR_acosh, false, "ACosH");
 DefineMonary (ATanH,      fRR_atanh, true,  "ATanH");
-DefineBinaryR(Drop,       fR2R_drop, true,  "Drop");
 DefineBinaryR(Sum,        fR2R_sum,  true,  "+ Real");
 DefineBinaryR(Rest,       fR2R_rest, true,  "- Real");
@@ -137 +136 @@
 };
 
+
+////////////////////////////////////////////////////////////////////////////////
+  int BVMat::Drop(double drop)
+//Importing method
+////////////////////////////////////////////////////////////////////////////////
+{
+  int r, c, k, result=0;
+  double* x;
+  switch(code_) {
+  case(ESC_blasRdense  ) : 
+  {
+    r = s_.blasRdense_->nrow;
+    c = s_.blasRdense_->ncol;
+    x = (double*)(s_.blasRdense_->x);
+    for (k=0; k<r*c; k++) 
+    {
+      if(fabs(*x)<=drop) { *x = 0.0; }
+    }    
+    break;
+  }
+  case(ESC_chlmRsparse ) : 
+  {
+    cholmod_drop(drop,s_.chlmRsparse_, common_);
+    cholmod_sort(s_.chlmRsparse_,   common_);
+    break;
+  }
+  case(ESC_chlmRtriplet) : 
+  {
+    BVMat sprs;
+    cRt2cRs(*this,sprs);
+    cholmod_drop(drop,sprs.s_.chlmRsparse_, common_);
+    cholmod_sort(sprs.s_.chlmRsparse_,   common_);
+    cRs2cRt(sprs,*this);
+    break; 
+  }
+  default: 
+    err_invalid_subtype("Drop",*this); 
+    result = -1; }
+  return(result);
+};
+
+////////////////////////////////////////////////////////////////////////////////
+  BVMat BVMat::Drop(const BVMat& A, double drop)
+////////////////////////////////////////////////////////////////////////////////
+{
+  BVMat v;
+  v.Copy(A);
+  v.Drop(drop);
+  return(v);
+}
 
 ////////////////////////////////////////////////////////////////////////////////