--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/liboctave/sparse-base-chol.h	Sun Oct 23 19:09:33 2005 +0000
@@ -0,0 +1,159 @@
+/*
+
+Copyright (C) 2005 David Bateman
+Copyright (C) 1998-2005 Andy Adler
+
+Octave is free software; you can redistribute it and/or modify it
+under the terms of the GNU General Public License as published by the
+Free Software Foundation; either version 2, or (at your option) any
+later version.
+
+Octave is distributed in the hope that it will be useful, but WITHOUT
+ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+for more details.
+
+You should have received a copy of the GNU General Public License
+along with this program; see the file COPYING.  If not, write to the
+Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor,
+Boston, MA 02110-1301, USA.
+
+*/
+
+#if !defined (octave_sparse_base_chol_h)
+#define octave_sparse_base_chol_h 1
+
+#ifdef HAVE_CONFIG_H
+#include <config.h>
+#endif
+
+#include "oct-sparse.h"
+#include "dColVector.h"
+
+template <class chol_type, class chol_elt, class p_type>
+class
+sparse_base_chol
+{
+protected:
+  class sparse_base_chol_rep
+  {
+  public:
+    sparse_base_chol_rep (void) : count (1), Lsparse (NULL), 
+				  is_pd (false), minor_p (0) { }
+
+    sparse_base_chol_rep (const chol_type& a, 
+			  const bool natural) : count (1)
+      { init (a, natural); }
+
+    sparse_base_chol_rep (const chol_type& a, octave_idx_type& info, 
+			  const bool natural) : count (1)
+      { info = init (a, natural); }
+
+#ifndef HAVE_CHOLMOD
+    ~sparse_base_chol_rep (void) { }
+#else
+    ~sparse_base_chol_rep (void) 
+      { CHOLMOD_NAME(free_sparse) (&Lsparse, &Common); }
+#endif
+
+    cholmod_sparse * L (void) const { return Lsparse; }
+
+    octave_idx_type P (void) const 
+      { return (minor_p == static_cast<octave_idx_type>(Lsparse->ncol) ? 
+		0 : minor_p + 1); }
+
+    ColumnVector perm (void) const { return perms + 1; }
+
+    p_type Q (void) const;
+
+    bool is_positive_definite (void) const { return is_pd; }
+
+    double rcond (void) const { return cond; }
+
+    int count;
+
+  private:
+    cholmod_sparse *Lsparse;
+
+    cholmod_common Common;
+
+    bool is_pd;
+
+    octave_idx_type minor_p;
+
+    ColumnVector perms;
+
+    double cond;
+
+    octave_idx_type init (const chol_type& a, bool natural = true);
+
+    void drop_zeros (const cholmod_sparse* S);
+
+    // No assignment
+    sparse_base_chol_rep& operator = (const sparse_base_chol_rep& a);
+  };
+
+ private:
+  sparse_base_chol_rep *rep;
+  
+public:
+
+  sparse_base_chol (void) : rep (new typename 
+    sparse_base_chol<chol_type, chol_elt, p_type>::sparse_base_chol_rep ()) { }
+
+  sparse_base_chol (const chol_type& a, const bool n) : rep (new typename 
+    sparse_base_chol<chol_type, chol_elt, p_type>::
+	sparse_base_chol_rep (a, n)) { }
+
+  sparse_base_chol (const chol_type& a, octave_idx_type& info, const bool n) :
+    rep (new typename sparse_base_chol<chol_type, chol_elt, p_type>::
+	sparse_base_chol_rep (a, info, n)) { }
+
+  sparse_base_chol (const sparse_base_chol<chol_type, chol_elt, p_type>& a) : 
+    rep (a.rep) { rep->count++; }
+
+  ~sparse_base_chol (void) 
+    {
+      if (--rep->count <= 0)
+	delete rep;
+    }
+
+  sparse_base_chol& operator = (const sparse_base_chol& a)
+    {
+      if (this != &a)
+	{
+	  if (--rep->count <= 0)
+	    delete rep;
+
+	  rep = a.rep;
+	  rep->count++;
+	}
+
+      return *this;
+    }
+
+  chol_type L (void) const;
+
+  chol_type R (void) const { return L().hermitian (); }
+
+  octave_idx_type P (void) const { return rep->P(); }
+
+  ColumnVector perm (void) const { return rep->perm(); }
+
+  p_type Q (void) const { return rep->Q(); }
+
+  bool is_positive_definite (void) const 
+    { return rep->is_positive_definite(); }
+
+  double rcond (void) const { return rep->rcond(); }
+
+  chol_type inverse (void) const;
+};
+
+#endif
+
+/*
+;;; Local Variables: ***
+;;; mode: C++ ***
+;;; End: ***
+*/