Changeset 8f6db64 in tol

Timestamp:

Feb 3, 2009, 9:11:48 AM (18 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

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

Children:

1d83266

Parents:

7c166a2

Message:

New static method force_natural_order
Setting configuration options for

Stop when non positive definite is found
Store allways as L*L', never as L*D*L'
Allow fast factor to sparse convertion

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

File:

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

tol/btol/vmatrix_type/vmat.cpp

@@ -72 +72 @@
   return(done);
 }  
+
+////////////////////////////////////////////////////////////////////////////////
+bool BVMat::force_natural_order(bool force)
+////////////////////////////////////////////////////////////////////////////////
+{
+  static bool force_ = false;
+  bool old = force_;
+  force_ = force;
+  if(force)
+  {
+    common_->nmethods            = 1 ;
+    common_->method[0].ordering  = CHOLMOD_NATURAL ;
+    common_->postorder           = false;
+  }
+  else
+  {
+    common_->nmethods           = 3 ;
+    common_->method[0].ordering = CHOLMOD_GIVEN ;
+    common_->method[1].ordering = CHOLMOD_AMD ;
+    common_->method[2].ordering = CHOLMOD_METIS ;
+    common_->postorder          = true;
+  }
+  return(old);
+}
+
 ////////////////////////////////////////////////////////////////////////////////
 void BVMat::restore_cholmod_common()
@@ -118 +143 @@
   common_->error_handler  = cholmod_error_handler;
   common_->try_catch      = NULL;
+
+  //Stop when non positive definite is found
+  common_->quick_return_if_not_posdef = true;
+
+  //Store allways as L*L', never as L*D*L'
+  common_->final_ll = true ;
+
+  //Needed to allow fast factor to sparse convertion
+  common_->final_pack      = true ;
+  common_->final_super     = false ;
+  common_->final_monotonic = true;
+
+  //We need specific final storement
+  common_->final_asis = false ;
+
+  //Natural order is not forced by default
+  force_natural_order(false);
+
 /* * /
 //MatLab lchol options