casa-mp-base

diff --git a/src/base-lu.cc b/src/base-lu.cc

new file mode 100644

index 0000000..c27d868

--- /dev/null

+++ b/src/base-lu.cc

@@ -0,0 +1,191 @@

+/*

+

+Copyright (C) 1996-2015 John W. Eaton

+Copyright (C) 2009 VZLU Prague

+

+This file is part of Octave.

+

+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 3 of the License, 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 Octave; see the file COPYING.  If not, see

+<http://www.gnu.org/licenses/>.

+

+*/

+

+#ifdef HAVE_CONFIG_H

+#include <config.h>

+#endif

+

+#include "base-lu.h"

+

+template <class lu_type>

+base_lu<lu_type>::base_lu (const lu_type& l, const lu_type& u,

+                           const PermMatrix& p)

+  : a_fact (u), l_fact (l), ipvt (p.transpose ().col_perm_vec ())

+{

+  if (l.columns () != u.rows ())

+    (*current_liboctave_error_handler) ("lu: dimension mismatch");

+}

+

+template <class lu_type>

+bool

+base_lu <lu_type> :: packed (void) const

+{

+  return l_fact.dims () == dim_vector ();

+}

+

+template <class lu_type>

+void

+base_lu <lu_type> :: unpack (void)

+{

+  if (packed ())

+    {

+      l_fact = L ();

+      a_fact = U (); // FIXME: sub-optimal

+      ipvt = getp ();

+    }

+}

+