Mercurial > octave
view libinterp/corefcn/xdiv.cc @ 16892:68fc671a9339
maint: Collapse interpfcn and interpfcn-core directories into corefcn directory.
* libgui/src/module.mk: Remove -I references to interp-core, interpfcn, add
reference to corefcn.
* libinterp/Makefile.am: Remove -I references to interp-core, interpfcn, add
reference to corefcn.
* libinterp/corefcn/module.mk: Add files from interp-core, interpfcn to
build system. Copy over special rules from module.mk files in interp-core
andd interpfcn.
* src/Makefile.am: Replace references to interp-core, interpfcn with those
to corefcn.
* libinterp/corefcn/Cell.cc, libinterp/corefcn/Cell.h,
libinterp/corefcn/action-container.h, libinterp/corefcn/c-file-ptr-stream.cc,
libinterp/corefcn/c-file-ptr-stream.h, libinterp/corefcn/comment-list.cc,
libinterp/corefcn/comment-list.h, libinterp/corefcn/cutils.c,
libinterp/corefcn/cutils.h, libinterp/corefcn/data.cc,
libinterp/corefcn/data.h, libinterp/corefcn/debug.cc,
libinterp/corefcn/debug.h, libinterp/corefcn/defaults.cc,
libinterp/corefcn/defaults.in.h, libinterp/corefcn/defun-dld.h,
libinterp/corefcn/defun-int.h, libinterp/corefcn/defun.cc,
libinterp/corefcn/defun.h, libinterp/corefcn/dirfns.cc,
libinterp/corefcn/dirfns.h, libinterp/corefcn/display.cc,
libinterp/corefcn/display.h, libinterp/corefcn/dynamic-ld.cc,
libinterp/corefcn/dynamic-ld.h, libinterp/corefcn/error.cc,
libinterp/corefcn/error.h, libinterp/corefcn/event-queue.h,
libinterp/corefcn/file-io.cc, libinterp/corefcn/file-io.h,
libinterp/corefcn/gl-render.cc, libinterp/corefcn/gl-render.h,
libinterp/corefcn/gl2ps-renderer.cc, libinterp/corefcn/gl2ps-renderer.h,
libinterp/corefcn/gl2ps.c, libinterp/corefcn/gl2ps.h,
libinterp/corefcn/graphics.cc, libinterp/corefcn/graphics.in.h,
libinterp/corefcn/gripes.cc, libinterp/corefcn/gripes.h,
libinterp/corefcn/help.cc, libinterp/corefcn/help.h,
libinterp/corefcn/hook-fcn.cc, libinterp/corefcn/hook-fcn.h,
libinterp/corefcn/input.cc, libinterp/corefcn/input.h,
libinterp/corefcn/jit-ir.cc, libinterp/corefcn/jit-ir.h,
libinterp/corefcn/jit-typeinfo.cc, libinterp/corefcn/jit-typeinfo.h,
libinterp/corefcn/jit-util.cc, libinterp/corefcn/jit-util.h,
libinterp/corefcn/load-path.cc, libinterp/corefcn/load-path.h,
libinterp/corefcn/load-save.cc, libinterp/corefcn/load-save.h,
libinterp/corefcn/ls-ascii-helper.cc, libinterp/corefcn/ls-ascii-helper.h,
libinterp/corefcn/ls-hdf5.cc, libinterp/corefcn/ls-hdf5.h,
libinterp/corefcn/ls-mat-ascii.cc, libinterp/corefcn/ls-mat-ascii.h,
libinterp/corefcn/ls-mat4.cc, libinterp/corefcn/ls-mat4.h,
libinterp/corefcn/ls-mat5.cc, libinterp/corefcn/ls-mat5.h,
libinterp/corefcn/ls-oct-ascii.cc, libinterp/corefcn/ls-oct-ascii.h,
libinterp/corefcn/ls-oct-binary.cc, libinterp/corefcn/ls-oct-binary.h,
libinterp/corefcn/ls-utils.cc, libinterp/corefcn/ls-utils.h,
libinterp/corefcn/matherr.c, libinterp/corefcn/mex.cc, libinterp/corefcn/mex.h,
libinterp/corefcn/mexproto.h, libinterp/corefcn/mxarray.in.h,
libinterp/corefcn/oct-errno.h, libinterp/corefcn/oct-errno.in.cc,
libinterp/corefcn/oct-fstrm.cc, libinterp/corefcn/oct-fstrm.h,
libinterp/corefcn/oct-hdf5.h, libinterp/corefcn/oct-hist.cc,
libinterp/corefcn/oct-hist.h, libinterp/corefcn/oct-iostrm.cc,
libinterp/corefcn/oct-iostrm.h, libinterp/corefcn/oct-lvalue.cc,
libinterp/corefcn/oct-lvalue.h, libinterp/corefcn/oct-map.cc,
libinterp/corefcn/oct-map.h, libinterp/corefcn/oct-obj.cc,
libinterp/corefcn/oct-obj.h, libinterp/corefcn/oct-prcstrm.cc,
libinterp/corefcn/oct-prcstrm.h, libinterp/corefcn/oct-procbuf.cc,
libinterp/corefcn/oct-procbuf.h, libinterp/corefcn/oct-stdstrm.h,
libinterp/corefcn/oct-stream.cc, libinterp/corefcn/oct-stream.h,
libinterp/corefcn/oct-strstrm.cc, libinterp/corefcn/oct-strstrm.h,
libinterp/corefcn/oct.h, libinterp/corefcn/octave-link.cc,
libinterp/corefcn/octave-link.h, libinterp/corefcn/pager.cc,
libinterp/corefcn/pager.h, libinterp/corefcn/pr-output.cc,
libinterp/corefcn/pr-output.h, libinterp/corefcn/procstream.cc,
libinterp/corefcn/procstream.h, libinterp/corefcn/profiler.cc,
libinterp/corefcn/profiler.h, libinterp/corefcn/pt-jit.cc,
libinterp/corefcn/pt-jit.h, libinterp/corefcn/sighandlers.cc,
libinterp/corefcn/sighandlers.h, libinterp/corefcn/siglist.c,
libinterp/corefcn/siglist.h, libinterp/corefcn/sparse-xdiv.cc,
libinterp/corefcn/sparse-xdiv.h, libinterp/corefcn/sparse-xpow.cc,
libinterp/corefcn/sparse-xpow.h, libinterp/corefcn/symtab.cc,
libinterp/corefcn/symtab.h, libinterp/corefcn/sysdep.cc,
libinterp/corefcn/sysdep.h, libinterp/corefcn/toplev.cc,
libinterp/corefcn/toplev.h, libinterp/corefcn/txt-eng-ft.cc,
libinterp/corefcn/txt-eng-ft.h, libinterp/corefcn/txt-eng.h,
libinterp/corefcn/unwind-prot.cc, libinterp/corefcn/unwind-prot.h,
libinterp/corefcn/utils.cc, libinterp/corefcn/utils.h,
libinterp/corefcn/variables.cc, libinterp/corefcn/variables.h,
libinterp/corefcn/workspace-element.h, libinterp/corefcn/xdiv.cc,
libinterp/corefcn/xdiv.h, libinterp/corefcn/xgl2ps.c,
libinterp/corefcn/xnorm.cc, libinterp/corefcn/xnorm.h,
libinterp/corefcn/xpow.cc, libinterp/corefcn/xpow.h,
libinterp/corefcn/zfstream.cc, libinterp/corefcn/zfstream.h:
Files moved from interp-core and interpfcn directories.
* libinterp/interp-core/Cell.cc, libinterp/interp-core/Cell.h,
libinterp/interp-core/action-container.h,
libinterp/interp-core/c-file-ptr-stream.cc,
libinterp/interp-core/c-file-ptr-stream.h,
libinterp/interp-core/comment-list.cc, libinterp/interp-core/comment-list.h,
libinterp/interp-core/cutils.c, libinterp/interp-core/cutils.h,
libinterp/interp-core/defun-dld.h, libinterp/interp-core/defun-int.h,
libinterp/interp-core/display.cc, libinterp/interp-core/display.h,
libinterp/interp-core/dynamic-ld.cc, libinterp/interp-core/dynamic-ld.h,
libinterp/interp-core/event-queue.h, libinterp/interp-core/gl-render.cc,
libinterp/interp-core/gl-render.h, libinterp/interp-core/gl2ps-renderer.cc,
libinterp/interp-core/gl2ps-renderer.h, libinterp/interp-core/gl2ps.c,
libinterp/interp-core/gl2ps.h, libinterp/interp-core/gripes.cc,
libinterp/interp-core/gripes.h, libinterp/interp-core/jit-ir.cc,
libinterp/interp-core/jit-ir.h, libinterp/interp-core/jit-typeinfo.cc,
libinterp/interp-core/jit-typeinfo.h, libinterp/interp-core/jit-util.cc,
libinterp/interp-core/jit-util.h, libinterp/interp-core/ls-ascii-helper.cc,
libinterp/interp-core/ls-ascii-helper.h, libinterp/interp-core/ls-hdf5.cc,
libinterp/interp-core/ls-hdf5.h, libinterp/interp-core/ls-mat-ascii.cc,
libinterp/interp-core/ls-mat-ascii.h, libinterp/interp-core/ls-mat4.cc,
libinterp/interp-core/ls-mat4.h, libinterp/interp-core/ls-mat5.cc,
libinterp/interp-core/ls-mat5.h, libinterp/interp-core/ls-oct-binary.cc,
libinterp/interp-core/ls-oct-binary.h, libinterp/interp-core/ls-utils.cc,
libinterp/interp-core/ls-utils.h, libinterp/interp-core/matherr.c,
libinterp/interp-core/mex.cc, libinterp/interp-core/mex.h,
libinterp/interp-core/mexproto.h, libinterp/interp-core/module.mk,
libinterp/interp-core/mxarray.in.h, libinterp/interp-core/oct-errno.h,
libinterp/interp-core/oct-errno.in.cc, libinterp/interp-core/oct-fstrm.cc,
libinterp/interp-core/oct-fstrm.h, libinterp/interp-core/oct-hdf5.h,
libinterp/interp-core/oct-iostrm.cc, libinterp/interp-core/oct-iostrm.h,
libinterp/interp-core/oct-lvalue.cc, libinterp/interp-core/oct-lvalue.h,
libinterp/interp-core/oct-map.cc, libinterp/interp-core/oct-map.h,
libinterp/interp-core/oct-obj.cc, libinterp/interp-core/oct-obj.h,
libinterp/interp-core/oct-prcstrm.cc, libinterp/interp-core/oct-prcstrm.h,
libinterp/interp-core/oct-procbuf.cc, libinterp/interp-core/oct-procbuf.h,
libinterp/interp-core/oct-stdstrm.h, libinterp/interp-core/oct-stream.cc,
libinterp/interp-core/oct-stream.h, libinterp/interp-core/oct-strstrm.cc,
libinterp/interp-core/oct-strstrm.h, libinterp/interp-core/oct.h,
libinterp/interp-core/procstream.cc, libinterp/interp-core/procstream.h,
libinterp/interp-core/pt-jit.cc, libinterp/interp-core/pt-jit.h,
libinterp/interp-core/siglist.c, libinterp/interp-core/siglist.h,
libinterp/interp-core/sparse-xdiv.cc, libinterp/interp-core/sparse-xdiv.h,
libinterp/interp-core/sparse-xpow.cc, libinterp/interp-core/sparse-xpow.h,
libinterp/interp-core/txt-eng-ft.cc, libinterp/interp-core/txt-eng-ft.h,
libinterp/interp-core/txt-eng.h, libinterp/interp-core/unwind-prot.cc,
libinterp/interp-core/unwind-prot.h, libinterp/interp-core/xdiv.cc,
libinterp/interp-core/xdiv.h, libinterp/interp-core/xgl2ps.c,
libinterp/interp-core/xnorm.cc, libinterp/interp-core/xnorm.h,
libinterp/interp-core/xpow.cc, libinterp/interp-core/xpow.h,
libinterp/interp-core/zfstream.cc, libinterp/interp-core/zfstream.h,
libinterp/interpfcn/data.cc, libinterp/interpfcn/data.h,
libinterp/interpfcn/debug.cc, libinterp/interpfcn/debug.h,
libinterp/interpfcn/defaults.cc, libinterp/interpfcn/defaults.in.h,
libinterp/interpfcn/defun.cc, libinterp/interpfcn/defun.h,
libinterp/interpfcn/dirfns.cc, libinterp/interpfcn/dirfns.h,
libinterp/interpfcn/error.cc, libinterp/interpfcn/error.h,
libinterp/interpfcn/file-io.cc, libinterp/interpfcn/file-io.h,
libinterp/interpfcn/graphics.cc, libinterp/interpfcn/graphics.in.h,
libinterp/interpfcn/help.cc, libinterp/interpfcn/help.h,
libinterp/interpfcn/hook-fcn.cc, libinterp/interpfcn/hook-fcn.h,
libinterp/interpfcn/input.cc, libinterp/interpfcn/input.h,
libinterp/interpfcn/load-path.cc, libinterp/interpfcn/load-path.h,
libinterp/interpfcn/load-save.cc, libinterp/interpfcn/load-save.h,
libinterp/interpfcn/ls-oct-ascii.cc, libinterp/interpfcn/ls-oct-ascii.h,
libinterp/interpfcn/module.mk, libinterp/interpfcn/oct-hist.cc,
libinterp/interpfcn/oct-hist.h, libinterp/interpfcn/octave-link.cc,
libinterp/interpfcn/octave-link.h, libinterp/interpfcn/pager.cc,
libinterp/interpfcn/pager.h, libinterp/interpfcn/pr-output.cc,
libinterp/interpfcn/pr-output.h, libinterp/interpfcn/profiler.cc,
libinterp/interpfcn/profiler.h, libinterp/interpfcn/sighandlers.cc,
libinterp/interpfcn/sighandlers.h, libinterp/interpfcn/symtab.cc,
libinterp/interpfcn/symtab.h, libinterp/interpfcn/sysdep.cc,
libinterp/interpfcn/sysdep.h, libinterp/interpfcn/toplev.cc,
libinterp/interpfcn/toplev.h, libinterp/interpfcn/utils.cc,
libinterp/interpfcn/utils.h, libinterp/interpfcn/variables.cc,
libinterp/interpfcn/variables.h, libinterp/interpfcn/workspace-element.h:
deleted files.
author | Rik <rik@octave.org> |
---|---|
date | Wed, 03 Jul 2013 17:43:48 -0700 |
parents | libinterp/interp-core/xdiv.cc@2fc554ffbc28 |
children | d63878346099 |
line wrap: on
line source
/* Copyright (C) 1993-2012 John W. Eaton Copyright (C) 2008 Jaroslav Hajek Copyright (C) 2009-2010 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 <cassert> #include "Array-util.h" #include "CMatrix.h" #include "dMatrix.h" #include "CNDArray.h" #include "dNDArray.h" #include "fCMatrix.h" #include "fMatrix.h" #include "fCNDArray.h" #include "fNDArray.h" #include "oct-cmplx.h" #include "dDiagMatrix.h" #include "fDiagMatrix.h" #include "CDiagMatrix.h" #include "fCDiagMatrix.h" #include "quit.h" #include "error.h" #include "xdiv.h" static inline bool result_ok (octave_idx_type info) { assert (info != -1); return (info != -2); } static void solve_singularity_warning (double rcond) { warning_with_id ("Octave:singular-matrix-div", "matrix singular to machine precision, rcond = %g", rcond); } template <class T1, class T2> bool mx_leftdiv_conform (const T1& a, const T2& b, blas_trans_type blas_trans) { octave_idx_type a_nr = blas_trans == blas_no_trans ? a.rows () : a.cols (); octave_idx_type b_nr = b.rows (); if (a_nr != b_nr) { octave_idx_type a_nc = blas_trans == blas_no_trans ? a.cols () : a.rows (); octave_idx_type b_nc = b.cols (); gripe_nonconformant ("operator \\", a_nr, a_nc, b_nr, b_nc); return false; } return true; } #define INSTANTIATE_MX_LEFTDIV_CONFORM(T1, T2) \ template bool mx_leftdiv_conform (const T1&, const T2&, blas_trans_type) INSTANTIATE_MX_LEFTDIV_CONFORM (Matrix, Matrix); INSTANTIATE_MX_LEFTDIV_CONFORM (Matrix, ComplexMatrix); INSTANTIATE_MX_LEFTDIV_CONFORM (ComplexMatrix, Matrix); INSTANTIATE_MX_LEFTDIV_CONFORM (ComplexMatrix, ComplexMatrix); template <class T1, class T2> bool mx_div_conform (const T1& a, const T2& b) { octave_idx_type a_nc = a.cols (); octave_idx_type b_nc = b.cols (); if (a_nc != b_nc) { octave_idx_type a_nr = a.rows (); octave_idx_type b_nr = b.rows (); gripe_nonconformant ("operator /", a_nr, a_nc, b_nr, b_nc); return false; } return true; } #define INSTANTIATE_MX_DIV_CONFORM(T1, T2) \ template bool mx_div_conform (const T1&, const T2&) INSTANTIATE_MX_DIV_CONFORM (Matrix, Matrix); INSTANTIATE_MX_DIV_CONFORM (Matrix, ComplexMatrix); INSTANTIATE_MX_DIV_CONFORM (ComplexMatrix, Matrix); INSTANTIATE_MX_DIV_CONFORM (ComplexMatrix, ComplexMatrix); // Right division functions. // // op2 / op1: m cm // +-- +---+----+ // matrix | 1 | 3 | // +---+----+ // complex_matrix | 2 | 4 | // +---+----+ // -*- 1 -*- Matrix xdiv (const Matrix& a, const Matrix& b, MatrixType &typ) { if (! mx_div_conform (a, b)) return Matrix (); octave_idx_type info; double rcond = 0.0; Matrix result = b.solve (typ, a.transpose (), info, rcond, solve_singularity_warning, true, blas_trans); return result.transpose (); } // -*- 2 -*- ComplexMatrix xdiv (const Matrix& a, const ComplexMatrix& b, MatrixType &typ) { if (! mx_div_conform (a, b)) return ComplexMatrix (); octave_idx_type info; double rcond = 0.0; ComplexMatrix result = b.solve (typ, a.transpose (), info, rcond, solve_singularity_warning, true, blas_trans); return result.transpose (); } // -*- 3 -*- ComplexMatrix xdiv (const ComplexMatrix& a, const Matrix& b, MatrixType &typ) { if (! mx_div_conform (a, b)) return ComplexMatrix (); octave_idx_type info; double rcond = 0.0; ComplexMatrix result = b.solve (typ, a.transpose (), info, rcond, solve_singularity_warning, true, blas_trans); return result.transpose (); } // -*- 4 -*- ComplexMatrix xdiv (const ComplexMatrix& a, const ComplexMatrix& b, MatrixType &typ) { if (! mx_div_conform (a, b)) return ComplexMatrix (); octave_idx_type info; double rcond = 0.0; ComplexMatrix result = b.solve (typ, a.transpose (), info, rcond, solve_singularity_warning, true, blas_trans); return result.transpose (); } // Funny element by element division operations. // // op2 \ op1: s cs // +-- +---+----+ // matrix | 1 | 3 | // +---+----+ // complex_matrix | 2 | 4 | // +---+----+ Matrix x_el_div (double a, const Matrix& b) { octave_idx_type nr = b.rows (); octave_idx_type nc = b.columns (); Matrix result (nr, nc); for (octave_idx_type j = 0; j < nc; j++) for (octave_idx_type i = 0; i < nr; i++) { octave_quit (); result (i, j) = a / b (i, j); } return result; } ComplexMatrix x_el_div (double a, const ComplexMatrix& b) { octave_idx_type nr = b.rows (); octave_idx_type nc = b.columns (); ComplexMatrix result (nr, nc); for (octave_idx_type j = 0; j < nc; j++) for (octave_idx_type i = 0; i < nr; i++) { octave_quit (); result (i, j) = a / b (i, j); } return result; } ComplexMatrix x_el_div (const Complex a, const Matrix& b) { octave_idx_type nr = b.rows (); octave_idx_type nc = b.columns (); ComplexMatrix result (nr, nc); for (octave_idx_type j = 0; j < nc; j++) for (octave_idx_type i = 0; i < nr; i++) { octave_quit (); result (i, j) = a / b (i, j); } return result; } ComplexMatrix x_el_div (const Complex a, const ComplexMatrix& b) { octave_idx_type nr = b.rows (); octave_idx_type nc = b.columns (); ComplexMatrix result (nr, nc); for (octave_idx_type j = 0; j < nc; j++) for (octave_idx_type i = 0; i < nr; i++) { octave_quit (); result (i, j) = a / b (i, j); } return result; } // Funny element by element division operations. // // op2 \ op1: s cs // +-- +---+----+ // N-d array | 1 | 3 | // +---+----+ // complex N-d array | 2 | 4 | // +---+----+ NDArray x_el_div (double a, const NDArray& b) { NDArray result (b.dims ()); for (octave_idx_type i = 0; i < b.length (); i++) { octave_quit (); result (i) = a / b (i); } return result; } ComplexNDArray x_el_div (double a, const ComplexNDArray& b) { ComplexNDArray result (b.dims ()); for (octave_idx_type i = 0; i < b.length (); i++) { octave_quit (); result (i) = a / b (i); } return result; } ComplexNDArray x_el_div (const Complex a, const NDArray& b) { ComplexNDArray result (b.dims ()); for (octave_idx_type i = 0; i < b.length (); i++) { octave_quit (); result (i) = a / b (i); } return result; } ComplexNDArray x_el_div (const Complex a, const ComplexNDArray& b) { ComplexNDArray result (b.dims ()); for (octave_idx_type i = 0; i < b.length (); i++) { octave_quit (); result (i) = a / b (i); } return result; } // Left division functions. // // op2 \ op1: m cm // +-- +---+----+ // matrix | 1 | 3 | // +---+----+ // complex_matrix | 2 | 4 | // +---+----+ // -*- 1 -*- Matrix xleftdiv (const Matrix& a, const Matrix& b, MatrixType &typ, blas_trans_type transt) { if (! mx_leftdiv_conform (a, b, transt)) return Matrix (); octave_idx_type info; double rcond = 0.0; return a.solve (typ, b, info, rcond, solve_singularity_warning, true, transt); } // -*- 2 -*- ComplexMatrix xleftdiv (const Matrix& a, const ComplexMatrix& b, MatrixType &typ, blas_trans_type transt) { if (! mx_leftdiv_conform (a, b, transt)) return ComplexMatrix (); octave_idx_type info; double rcond = 0.0; return a.solve (typ, b, info, rcond, solve_singularity_warning, true, transt); } // -*- 3 -*- ComplexMatrix xleftdiv (const ComplexMatrix& a, const Matrix& b, MatrixType &typ, blas_trans_type transt) { if (! mx_leftdiv_conform (a, b, transt)) return ComplexMatrix (); octave_idx_type info; double rcond = 0.0; return a.solve (typ, b, info, rcond, solve_singularity_warning, true, transt); } // -*- 4 -*- ComplexMatrix xleftdiv (const ComplexMatrix& a, const ComplexMatrix& b, MatrixType &typ, blas_trans_type transt) { if (! mx_leftdiv_conform (a, b, transt)) return ComplexMatrix (); octave_idx_type info; double rcond = 0.0; return a.solve (typ, b, info, rcond, solve_singularity_warning, true, transt); } static void solve_singularity_warning (float rcond) { warning ("matrix singular to machine precision, rcond = %g", rcond); warning ("attempting to find minimum norm solution"); } INSTANTIATE_MX_LEFTDIV_CONFORM (FloatMatrix, FloatMatrix); INSTANTIATE_MX_LEFTDIV_CONFORM (FloatMatrix, FloatComplexMatrix); INSTANTIATE_MX_LEFTDIV_CONFORM (FloatComplexMatrix, FloatMatrix); INSTANTIATE_MX_LEFTDIV_CONFORM (FloatComplexMatrix, FloatComplexMatrix); INSTANTIATE_MX_DIV_CONFORM (FloatMatrix, FloatMatrix); INSTANTIATE_MX_DIV_CONFORM (FloatMatrix, FloatComplexMatrix); INSTANTIATE_MX_DIV_CONFORM (FloatComplexMatrix, FloatMatrix); INSTANTIATE_MX_DIV_CONFORM (FloatComplexMatrix, FloatComplexMatrix); // Right division functions. // // op2 / op1: m cm // +-- +---+----+ // matrix | 1 | 3 | // +---+----+ // complex_matrix | 2 | 4 | // +---+----+ // -*- 1 -*- FloatMatrix xdiv (const FloatMatrix& a, const FloatMatrix& b, MatrixType &typ) { if (! mx_div_conform (a, b)) return FloatMatrix (); octave_idx_type info; float rcond = 0.0; FloatMatrix result = b.solve (typ, a.transpose (), info, rcond, solve_singularity_warning, true, blas_trans); return result.transpose (); } // -*- 2 -*- FloatComplexMatrix xdiv (const FloatMatrix& a, const FloatComplexMatrix& b, MatrixType &typ) { if (! mx_div_conform (a, b)) return FloatComplexMatrix (); octave_idx_type info; float rcond = 0.0; FloatComplexMatrix result = b.solve (typ, a.transpose (), info, rcond, solve_singularity_warning, true, blas_trans); return result.transpose (); } // -*- 3 -*- FloatComplexMatrix xdiv (const FloatComplexMatrix& a, const FloatMatrix& b, MatrixType &typ) { if (! mx_div_conform (a, b)) return FloatComplexMatrix (); octave_idx_type info; float rcond = 0.0; FloatComplexMatrix result = b.solve (typ, a.transpose (), info, rcond, solve_singularity_warning, true, blas_trans); return result.transpose (); } // -*- 4 -*- FloatComplexMatrix xdiv (const FloatComplexMatrix& a, const FloatComplexMatrix& b, MatrixType &typ) { if (! mx_div_conform (a, b)) return FloatComplexMatrix (); octave_idx_type info; float rcond = 0.0; FloatComplexMatrix result = b.solve (typ, a.transpose (), info, rcond, solve_singularity_warning, true, blas_trans); return result.transpose (); } // Funny element by element division operations. // // op2 \ op1: s cs // +-- +---+----+ // matrix | 1 | 3 | // +---+----+ // complex_matrix | 2 | 4 | // +---+----+ FloatMatrix x_el_div (float a, const FloatMatrix& b) { octave_idx_type nr = b.rows (); octave_idx_type nc = b.columns (); FloatMatrix result (nr, nc); for (octave_idx_type j = 0; j < nc; j++) for (octave_idx_type i = 0; i < nr; i++) { octave_quit (); result (i, j) = a / b (i, j); } return result; } FloatComplexMatrix x_el_div (float a, const FloatComplexMatrix& b) { octave_idx_type nr = b.rows (); octave_idx_type nc = b.columns (); FloatComplexMatrix result (nr, nc); for (octave_idx_type j = 0; j < nc; j++) for (octave_idx_type i = 0; i < nr; i++) { octave_quit (); result (i, j) = a / b (i, j); } return result; } FloatComplexMatrix x_el_div (const FloatComplex a, const FloatMatrix& b) { octave_idx_type nr = b.rows (); octave_idx_type nc = b.columns (); FloatComplexMatrix result (nr, nc); for (octave_idx_type j = 0; j < nc; j++) for (octave_idx_type i = 0; i < nr; i++) { octave_quit (); result (i, j) = a / b (i, j); } return result; } FloatComplexMatrix x_el_div (const FloatComplex a, const FloatComplexMatrix& b) { octave_idx_type nr = b.rows (); octave_idx_type nc = b.columns (); FloatComplexMatrix result (nr, nc); for (octave_idx_type j = 0; j < nc; j++) for (octave_idx_type i = 0; i < nr; i++) { octave_quit (); result (i, j) = a / b (i, j); } return result; } // Funny element by element division operations. // // op2 \ op1: s cs // +-- +---+----+ // N-d array | 1 | 3 | // +---+----+ // complex N-d array | 2 | 4 | // +---+----+ FloatNDArray x_el_div (float a, const FloatNDArray& b) { FloatNDArray result (b.dims ()); for (octave_idx_type i = 0; i < b.length (); i++) { octave_quit (); result (i) = a / b (i); } return result; } FloatComplexNDArray x_el_div (float a, const FloatComplexNDArray& b) { FloatComplexNDArray result (b.dims ()); for (octave_idx_type i = 0; i < b.length (); i++) { octave_quit (); result (i) = a / b (i); } return result; } FloatComplexNDArray x_el_div (const FloatComplex a, const FloatNDArray& b) { FloatComplexNDArray result (b.dims ()); for (octave_idx_type i = 0; i < b.length (); i++) { octave_quit (); result (i) = a / b (i); } return result; } FloatComplexNDArray x_el_div (const FloatComplex a, const FloatComplexNDArray& b) { FloatComplexNDArray result (b.dims ()); for (octave_idx_type i = 0; i < b.length (); i++) { octave_quit (); result (i) = a / b (i); } return result; } // Left division functions. // // op2 \ op1: m cm // +-- +---+----+ // matrix | 1 | 3 | // +---+----+ // complex_matrix | 2 | 4 | // +---+----+ // -*- 1 -*- FloatMatrix xleftdiv (const FloatMatrix& a, const FloatMatrix& b, MatrixType &typ, blas_trans_type transt) { if (! mx_leftdiv_conform (a, b, transt)) return FloatMatrix (); octave_idx_type info; float rcond = 0.0; return a.solve (typ, b, info, rcond, solve_singularity_warning, true, transt); } // -*- 2 -*- FloatComplexMatrix xleftdiv (const FloatMatrix& a, const FloatComplexMatrix& b, MatrixType &typ, blas_trans_type transt) { if (! mx_leftdiv_conform (a, b, transt)) return FloatComplexMatrix (); octave_idx_type info; float rcond = 0.0; return a.solve (typ, b, info, rcond, solve_singularity_warning, true, transt); } // -*- 3 -*- FloatComplexMatrix xleftdiv (const FloatComplexMatrix& a, const FloatMatrix& b, MatrixType &typ, blas_trans_type transt) { if (! mx_leftdiv_conform (a, b, transt)) return FloatComplexMatrix (); octave_idx_type info; float rcond = 0.0; return a.solve (typ, b, info, rcond, solve_singularity_warning, true, transt); } // -*- 4 -*- FloatComplexMatrix xleftdiv (const FloatComplexMatrix& a, const FloatComplexMatrix& b, MatrixType &typ, blas_trans_type transt) { if (! mx_leftdiv_conform (a, b, transt)) return FloatComplexMatrix (); octave_idx_type info; float rcond = 0.0; return a.solve (typ, b, info, rcond, solve_singularity_warning, true, transt); } // Diagonal matrix division. template <class MT, class DMT> MT mdm_div_impl (const MT& a, const DMT& d) { if (! mx_div_conform (a, d)) return MT (); octave_idx_type m = a.rows (), n = d.rows (), l = d.length (); MT x (m, n); typedef typename DMT::element_type S; typedef typename MT::element_type T; const T *aa = a.data (); const S *dd = d.data (); T *xx = x.fortran_vec (); for (octave_idx_type j = 0; j < l; j++) { const S del = dd[j]; if (del != S ()) for (octave_idx_type i = 0; i < m; i++) xx[i] = aa[i] / del; else for (octave_idx_type i = 0; i < m; i++) xx[i] = T (); aa += m; xx += m; } for (octave_idx_type i = l*m; i < n*m; i++) xx[i] = T (); return x; } // Right division functions. // // op2 / op1: dm cdm // +-- +---+----+ // matrix | 1 | | // +---+----+ // complex_matrix | 2 | 3 | // +---+----+ // -*- 1 -*- Matrix xdiv (const Matrix& a, const DiagMatrix& b) { return mdm_div_impl (a, b); } // -*- 2 -*- ComplexMatrix xdiv (const ComplexMatrix& a, const DiagMatrix& b) { return mdm_div_impl (a, b); } // -*- 3 -*- ComplexMatrix xdiv (const ComplexMatrix& a, const ComplexDiagMatrix& b) { return mdm_div_impl (a, b); } // Right division functions, float type. // // op2 / op1: dm cdm // +-- +---+----+ // matrix | 1 | | // +---+----+ // complex_matrix | 2 | 3 | // +---+----+ // -*- 1 -*- FloatMatrix xdiv (const FloatMatrix& a, const FloatDiagMatrix& b) { return mdm_div_impl (a, b); } // -*- 2 -*- FloatComplexMatrix xdiv (const FloatComplexMatrix& a, const FloatDiagMatrix& b) { return mdm_div_impl (a, b); } // -*- 3 -*- FloatComplexMatrix xdiv (const FloatComplexMatrix& a, const FloatComplexDiagMatrix& b) { return mdm_div_impl (a, b); } template <class MT, class DMT> MT dmm_leftdiv_impl (const DMT& d, const MT& a) { if (! mx_leftdiv_conform (d, a, blas_no_trans)) return MT (); octave_idx_type m = d.cols (), n = a.cols (), k = a.rows (), l = d.length (); MT x (m, n); typedef typename DMT::element_type S; typedef typename MT::element_type T; const T *aa = a.data (); const S *dd = d.data (); T *xx = x.fortran_vec (); for (octave_idx_type j = 0; j < n; j++) { for (octave_idx_type i = 0; i < l; i++) xx[i] = dd[i] != S () ? aa[i] / dd[i] : T (); for (octave_idx_type i = l; i < m; i++) xx[i] = T (); aa += k; xx += m; } return x; } // Left division functions. // // op2 \ op1: m cm // +---+----+ // diag_matrix | 1 | 2 | // +---+----+ // complex_diag_matrix | | 3 | // +---+----+ // -*- 1 -*- Matrix xleftdiv (const DiagMatrix& a, const Matrix& b) { return dmm_leftdiv_impl (a, b); } // -*- 2 -*- ComplexMatrix xleftdiv (const DiagMatrix& a, const ComplexMatrix& b) { return dmm_leftdiv_impl (a, b); } // -*- 3 -*- ComplexMatrix xleftdiv (const ComplexDiagMatrix& a, const ComplexMatrix& b) { return dmm_leftdiv_impl (a, b); } // Left division functions, float type. // // op2 \ op1: m cm // +---+----+ // diag_matrix | 1 | 2 | // +---+----+ // complex_diag_matrix | | 3 | // +---+----+ // -*- 1 -*- FloatMatrix xleftdiv (const FloatDiagMatrix& a, const FloatMatrix& b) { return dmm_leftdiv_impl (a, b); } // -*- 2 -*- FloatComplexMatrix xleftdiv (const FloatDiagMatrix& a, const FloatComplexMatrix& b) { return dmm_leftdiv_impl (a, b); } // -*- 3 -*- FloatComplexMatrix xleftdiv (const FloatComplexDiagMatrix& a, const FloatComplexMatrix& b) { return dmm_leftdiv_impl (a, b); } // Diagonal by diagonal matrix division. template <class MT, class DMT> MT dmdm_div_impl (const MT& a, const DMT& d) { if (! mx_div_conform (a, d)) return MT (); octave_idx_type m = a.rows (), n = d.rows (), k = d.cols (); octave_idx_type l = std::min (m, n), lk = std::min (l, k); MT x (m, n); typedef typename DMT::element_type S; typedef typename MT::element_type T; const T *aa = a.data (); const S *dd = d.data (); T *xx = x.fortran_vec (); for (octave_idx_type i = 0; i < lk; i++) xx[i] = dd[i] != S () ? aa[i] / dd[i] : T (); for (octave_idx_type i = lk; i < l; i++) xx[i] = T (); return x; } // Right division functions. // // op2 / op1: dm cdm // +-- +---+----+ // diag_matrix | 1 | | // +---+----+ // complex_diag_matrix | 2 | 3 | // +---+----+ // -*- 1 -*- DiagMatrix xdiv (const DiagMatrix& a, const DiagMatrix& b) { return dmdm_div_impl (a, b); } // -*- 2 -*- ComplexDiagMatrix xdiv (const ComplexDiagMatrix& a, const DiagMatrix& b) { return dmdm_div_impl (a, b); } // -*- 3 -*- ComplexDiagMatrix xdiv (const ComplexDiagMatrix& a, const ComplexDiagMatrix& b) { return dmdm_div_impl (a, b); } // Right division functions, float type. // // op2 / op1: dm cdm // +-- +---+----+ // diag_matrix | 1 | | // +---+----+ // complex_diag_matrix | 2 | 3 | // +---+----+ // -*- 1 -*- FloatDiagMatrix xdiv (const FloatDiagMatrix& a, const FloatDiagMatrix& b) { return dmdm_div_impl (a, b); } // -*- 2 -*- FloatComplexDiagMatrix xdiv (const FloatComplexDiagMatrix& a, const FloatDiagMatrix& b) { return dmdm_div_impl (a, b); } // -*- 3 -*- FloatComplexDiagMatrix xdiv (const FloatComplexDiagMatrix& a, const FloatComplexDiagMatrix& b) { return dmdm_div_impl (a, b); } template <class MT, class DMT> MT dmdm_leftdiv_impl (const DMT& d, const MT& a) { if (! mx_leftdiv_conform (d, a, blas_no_trans)) return MT (); octave_idx_type m = d.cols (), n = a.cols (), k = d.rows (); octave_idx_type l = std::min (m, n), lk = std::min (l, k); MT x (m, n); typedef typename DMT::element_type S; typedef typename MT::element_type T; const T *aa = a.data (); const S *dd = d.data (); T *xx = x.fortran_vec (); for (octave_idx_type i = 0; i < lk; i++) xx[i] = dd[i] != S () ? aa[i] / dd[i] : T (); for (octave_idx_type i = lk; i < l; i++) xx[i] = T (); return x; } // Left division functions. // // op2 \ op1: dm cdm // +---+----+ // diag_matrix | 1 | 2 | // +---+----+ // complex_diag_matrix | | 3 | // +---+----+ // -*- 1 -*- DiagMatrix xleftdiv (const DiagMatrix& a, const DiagMatrix& b) { return dmdm_leftdiv_impl (a, b); } // -*- 2 -*- ComplexDiagMatrix xleftdiv (const DiagMatrix& a, const ComplexDiagMatrix& b) { return dmdm_leftdiv_impl (a, b); } // -*- 3 -*- ComplexDiagMatrix xleftdiv (const ComplexDiagMatrix& a, const ComplexDiagMatrix& b) { return dmdm_leftdiv_impl (a, b); } // Left division functions, float type. // // op2 \ op1: dm cdm // +---+----+ // diag_matrix | 1 | 2 | // +---+----+ // complex_diag_matrix | | 3 | // +---+----+ // -*- 1 -*- FloatDiagMatrix xleftdiv (const FloatDiagMatrix& a, const FloatDiagMatrix& b) { return dmdm_leftdiv_impl (a, b); } // -*- 2 -*- FloatComplexDiagMatrix xleftdiv (const FloatDiagMatrix& a, const FloatComplexDiagMatrix& b) { return dmdm_leftdiv_impl (a, b); } // -*- 3 -*- FloatComplexDiagMatrix xleftdiv (const FloatComplexDiagMatrix& a, const FloatComplexDiagMatrix& b) { return dmdm_leftdiv_impl (a, b); }