Mercurial > octave
view liboctave/util/sun-utils.h @ 21176:791dcb32b657
revamp sparse QR factorizatino classes
* sparse-qr-inst.cc, sparse-qr.cc, sparse-qr.h: New files, adapted
from SparseCmplxQR.cc, SparseCmplxQR.h, SparseQR.cc, and SparseQR.h.
(sparse_qr): New template class.
* SparseCmplxQR.cc, SparseCmplxQR.h, SparseQR.cc, SparseQR.h: Delete.
* liboctave/numeric/module.mk: Update.
* dmperm.cc, qr.cc, CSparse.cc, dSparse.cc, sparse-dmsolve.cc:
Use new sparse_qr class.
* oct-sparse.h (CXSPARSE_DNAME, CXSPARSE_ZNAME):
Move macros here from SparseQR.h and SparseCmplxQR.h.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Tue, 02 Feb 2016 12:34:06 -0500 |
parents | f7084eae3318 |
children | 1473547f50f5 |
line wrap: on
line source
/* Copyright (C) 1993-2015 John W. Eaton 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/>. */ #if ! defined (octave_sun_utils_h) #define octave_sun_utils_h 1 // This is only needed to dereference pointers to doubles if mixing // GCC and Sun SPARC f77/cc compiled code. See the GCC manual (where the // function access_double() is described) and the Sun f77 manual, // which explains that doubles are not always aligned on 8 byte // boundaries. #if defined (__sparc) && defined (__GNUC__) inline double access_double (double *unaligned_ptr) { union d2i { double d; int i[2]; }; union d2i *p = (union d2i *) unaligned_ptr; union d2i u; u.i[0] = p->i[0]; u.i[1] = p->i[1]; return u.d; } inline void assign_double (double *unaligned_ptr, double value) { union d2i { double d; int i[2]; }; double *ptr = &value; union d2i *v = (union d2i *) ptr; union d2i *p = (union d2i *) unaligned_ptr; p->i[0] = v->i[0]; p->i[1] = v->i[1]; } #endif #endif