octave: libinterp/corefcn/__pchip_deriv_

annotate libinterp/corefcn/__pchip_deriv__.cc @ 22322:93b3cdd36854

move most f77 function decls to separate header files * liboctave/numeric/lo-amos-proto.h, liboctave/numeric/lo-arpack-proto.h, liboctave/numeric/lo-blas-proto.h, liboctave/numeric/lo-fftpack-proto.h, liboctave/numeric/lo-lapack-proto.h, liboctave/numeric/lo-qrupdate-proto.h, liboctave/numeric/lo-ranlib-proto.h, liboctave/numeric/lo-slatec-proto.h: New files. * liboctave/numeric/module.mk: Update. * __pchip_deriv__.cc, dot.cc, interpreter.cc, ordschur.cc, qz.cc, CColVector.cc, CMatrix.cc, CNDArray.cc, CRowVector.cc, CSparse.cc, dColVector.cc, dMatrix.cc, dNDArray.cc, dRowVector.cc, dSparse.cc, fCColVector.cc, fCMatrix.cc, fCNDArray.cc, fCRowVector.cc, fColVector.cc, fMatrix.cc, fNDArray.cc, fRowVector.cc, EIG.cc, aepbalance.cc, chol.cc, eigs-base.cc, fEIG.cc, gepbalance.cc, gsvd.cc, hess.cc, lo-specfun.cc, lu.cc, oct-rand.cc, qr.cc, qrp.cc, randpoisson.cc, schur.cc, sparse-qr.cc, svd.cc: Use new header files.

author	John W. Eaton <jwe@octave.org>
date	Wed, 17 Aug 2016 00:18:08 -0400
parents	59cadee1c74b
children	bac0d6f07a3e

rev	line source
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11553 diff changeset	1 /*
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	2
19697 4197fc428c7d maint: Update copyright notices for 2015. John W. Eaton <jwe@octave.org> parents: 18605 diff changeset	3 Copyright (C) 2002-2015 Kai Habel
11523 fd0a3ac60b0e update copyright notices John W. Eaton <jwe@octave.org> parents: 11518 diff changeset	4 Copyright (C) 2008-2009 Jaroslav Hajek
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	5
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	6 This file is part of Octave.
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	7
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	8 Octave is free software; you can redistribute it and/or modify it
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	9 under the terms of the GNU General Public License as published by the
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6945 diff changeset	10 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6945 diff changeset	11 option) any later version.
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	12
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	13 Octave is distributed in the hope that it will be useful, but WITHOUT
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	16 for more details.
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	17
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	18 You should have received a copy of the GNU General Public License
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6945 diff changeset	19 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6945 diff changeset	20 <http://www.gnu.org/licenses/>.
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	21
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	22 */
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	23
21724 aba2e6293dd8 use "#if ..." consistently instead of "#ifdef" and "#ifndef" John W. Eaton <jwe@octave.org> parents: 21301 diff changeset	24 #if defined (HAVE_CONFIG_H)
21301 40de9f8f23a6 Use '#include "config.h"' rather than <config.h>. Rik <rik@octave.org> parents: 21200 diff changeset	25 # include "config.h"
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	26 #endif
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	27
22322 93b3cdd36854 move most f77 function decls to separate header files John W. Eaton <jwe@octave.org> parents: 22133 diff changeset	28 #include "lo-slatec-proto.h"
93b3cdd36854 move most f77 function decls to separate header files John W. Eaton <jwe@octave.org> parents: 22133 diff changeset	29
15039 e753177cde93 maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory Rik <rik@octave.org> parents: 14846 diff changeset	30 #include "defun.h"
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	31 #include "error.h"
21100 e39e05d90788 Switch gripe_XXX to either err_XXX or warn_XXX naming scheme. Rik <rik@octave.org> parents: 21078 diff changeset	32 #include "errwarn.h"
20940 48b2ad5ee801 maint: Rename oct-obj.[cc\|h] to ovl.[cc\|h] for clarity. Rik <rik@octave.org> parents: 20853 diff changeset	33 #include "ovl.h"
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	34 #include "utils.h"
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	35 #include "f77-fcn.h"
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	36
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6549 diff changeset	37 // Wrapper for SLATEC/PCHIP function DPCHIM to calculate the derivates
6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6549 diff changeset	38 // for piecewise polynomials.
6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6549 diff changeset	39
15039 e753177cde93 maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory Rik <rik@octave.org> parents: 14846 diff changeset	40 DEFUN (__pchip_deriv__, args, ,
21966 112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	41 doc: /* -- texinfo --
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	42 @deftypefn {} {} __pchip_deriv__ (@var{x}, @var{y}, @var{dim})
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	43 Undocumented internal function.
112b20240c87 move docstrings in C++ files out of C strings and into comments John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	44 @end deftypefn */)
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	45 {
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	46 octave_value retval;
21078 49852ff04747 maint: Remove unnecessary declarations of retval. Rik <rik@octave.org> parents: 20940 diff changeset	47 int nargin = args.length ();
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	48
18112 b560bac0fca2 maint: Don't use space between 'args' and '(' when doing indexing. Rik <rik@octave.org> parents: 17787 diff changeset	49 bool rows = (nargin == 3 && args(2).uint_value () == 2);
8712 010e15c7be9a support pchip method in interp2 Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	50
010e15c7be9a support pchip method in interp2 Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	51 if (nargin >= 2)
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	52 {
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	53 if (args(0).is_single_type () \|\| args(1).is_single_type ())
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	54 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	55 FloatColumnVector xvec (args(0).float_vector_value ());
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	56 FloatMatrix ymat (args(1).float_matrix_value ());
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	57
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19873 diff changeset	58 octave_idx_type nx = xvec.numel ();
14462 af552038cc52 Fix segfault in pchip when input array is too small (bug #35835). Jordi Gutiérrez Hermoso <jordigh@octave.org> parents: 14138 diff changeset	59
af552038cc52 Fix segfault in pchip when input array is too small (bug #35835). Jordi Gutiérrez Hermoso <jordigh@octave.org> parents: 14138 diff changeset	60 if (nx < 2)
20831 35241c4b696c eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20232 diff changeset	61 error ("__pchip_deriv__: X must be at least of length 2");
14462 af552038cc52 Fix segfault in pchip when input array is too small (bug #35835). Jordi Gutiérrez Hermoso <jordigh@octave.org> parents: 14138 diff changeset	62
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	63 octave_idx_type nyr = ymat.rows ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	64 octave_idx_type nyc = ymat.columns ();
5838 376e02b2ce70 [project @ 2006-06-01 20:23:53 by jwe] jwe parents: 5837 diff changeset	65
8712 010e15c7be9a support pchip method in interp2 Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	66 if (nx != (rows ? nyc : nyr))
20831 35241c4b696c eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20232 diff changeset	67 error ("__pchip_deriv__: X and Y dimension mismatch");
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	68
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	69 FloatMatrix dmat (nyr, nyc);
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	70
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	71 octave_idx_type ierr;
8712 010e15c7be9a support pchip method in interp2 Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	72 const octave_idx_type incfd = rows ? nyr : 1;
19787 4ae2987c5f66 avoid variable possibly clobbered by longjmp warning from GCC John W. Eaton <jwe@octave.org> parents: 19780 diff changeset	73 volatile const octave_idx_type inc = rows ? 1 : nyr;
4ae2987c5f66 avoid variable possibly clobbered by longjmp warning from GCC John W. Eaton <jwe@octave.org> parents: 19780 diff changeset	74 volatile octave_idx_type k = 0;
8712 010e15c7be9a support pchip method in interp2 Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	75
19873 5261186f1b00 __pchip_deriv__.cc: Stop variable might be clobbered by 'longjmp' warning. Rik <rik@octave.org> parents: 19787 diff changeset	76 for (volatile octave_idx_type i = (rows ? nyr : nyc); i > 0; i--)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	77 {
18603 e7369c43ebe0 Fix segmentation fault for interp2 and pchip method (bug #41838). Stefan Mahr <dac922@gmx.de> parents: 17787 diff changeset	78 F77_XFCN (pchim, PCHIM, (nx, xvec.data (),
19787 4ae2987c5f66 avoid variable possibly clobbered by longjmp warning from GCC John W. Eaton <jwe@octave.org> parents: 19780 diff changeset	79 ymat.data () + k * inc,
4ae2987c5f66 avoid variable possibly clobbered by longjmp warning from GCC John W. Eaton <jwe@octave.org> parents: 19780 diff changeset	80 dmat.fortran_vec () + k * inc,
4ae2987c5f66 avoid variable possibly clobbered by longjmp warning from GCC John W. Eaton <jwe@octave.org> parents: 19780 diff changeset	81 incfd, &ierr));
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	82
19787 4ae2987c5f66 avoid variable possibly clobbered by longjmp warning from GCC John W. Eaton <jwe@octave.org> parents: 19780 diff changeset	83 k++;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	84
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	85 if (ierr < 0)
20831 35241c4b696c eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20232 diff changeset	86 error ("__pchip_deriv__: PCHIM failed with ierr = %i", ierr);
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	87 }
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	88
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	89 retval = dmat;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	90 }
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	91 else
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	92 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	93 ColumnVector xvec (args(0).vector_value ());
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	94 Matrix ymat (args(1).matrix_value ());
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	95
20232 a9574e3c6e9e Deprecate Array::length() and Sparse::length() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 19873 diff changeset	96 octave_idx_type nx = xvec.numel ();
14462 af552038cc52 Fix segfault in pchip when input array is too small (bug #35835). Jordi Gutiérrez Hermoso <jordigh@octave.org> parents: 14138 diff changeset	97
af552038cc52 Fix segfault in pchip when input array is too small (bug #35835). Jordi Gutiérrez Hermoso <jordigh@octave.org> parents: 14138 diff changeset	98 if (nx < 2)
20831 35241c4b696c eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20232 diff changeset	99 error ("__pchip_deriv__: X must be at least of length 2");
14462 af552038cc52 Fix segfault in pchip when input array is too small (bug #35835). Jordi Gutiérrez Hermoso <jordigh@octave.org> parents: 14138 diff changeset	100
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	101 octave_idx_type nyr = ymat.rows ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	102 octave_idx_type nyc = ymat.columns ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	103
8712 010e15c7be9a support pchip method in interp2 Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	104 if (nx != (rows ? nyc : nyr))
20831 35241c4b696c eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20232 diff changeset	105 error ("__pchip_deriv__: X and Y dimension mismatch");
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	106
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	107 Matrix dmat (nyr, nyc);
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	108
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	109 octave_idx_type ierr;
8712 010e15c7be9a support pchip method in interp2 Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	110 const octave_idx_type incfd = rows ? nyr : 1;
19787 4ae2987c5f66 avoid variable possibly clobbered by longjmp warning from GCC John W. Eaton <jwe@octave.org> parents: 19780 diff changeset	111 volatile const octave_idx_type inc = rows ? 1 : nyr;
4ae2987c5f66 avoid variable possibly clobbered by longjmp warning from GCC John W. Eaton <jwe@octave.org> parents: 19780 diff changeset	112 volatile octave_idx_type k = 0;
8712 010e15c7be9a support pchip method in interp2 Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	113
19873 5261186f1b00 __pchip_deriv__.cc: Stop variable might be clobbered by 'longjmp' warning. Rik <rik@octave.org> parents: 19787 diff changeset	114 for (volatile octave_idx_type i = (rows ? nyr : nyc); i > 0; i--)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	115 {
18603 e7369c43ebe0 Fix segmentation fault for interp2 and pchip method (bug #41838). Stefan Mahr <dac922@gmx.de> parents: 17787 diff changeset	116 F77_XFCN (dpchim, DPCHIM, (nx, xvec.data (),
19787 4ae2987c5f66 avoid variable possibly clobbered by longjmp warning from GCC John W. Eaton <jwe@octave.org> parents: 19780 diff changeset	117 ymat.data () + k * inc,
4ae2987c5f66 avoid variable possibly clobbered by longjmp warning from GCC John W. Eaton <jwe@octave.org> parents: 19780 diff changeset	118 dmat.fortran_vec () + k * inc,
4ae2987c5f66 avoid variable possibly clobbered by longjmp warning from GCC John W. Eaton <jwe@octave.org> parents: 19780 diff changeset	119 incfd, &ierr));
4ae2987c5f66 avoid variable possibly clobbered by longjmp warning from GCC John W. Eaton <jwe@octave.org> parents: 19780 diff changeset	120 k++;
5838 376e02b2ce70 [project @ 2006-06-01 20:23:53 by jwe] jwe parents: 5837 diff changeset	121
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	122 if (ierr < 0)
20831 35241c4b696c eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20232 diff changeset	123 error ("__pchip_deriv__: DPCHIM failed with ierr = %i", ierr);
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	124 }
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	125
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	126 retval = dmat;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	127 }
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	128 }
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	129
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	130 return retval;
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	131 }
12805 3641167e5b75 codesprint: .cc helper functions do not need tests Rik <octave@nomad.inbox5.com>* parents: 11586 diff changeset	132
3641167e5b75 codesprint: .cc helper functions do not need tests Rik <octave@nomad.inbox5.com>* parents: 11586 diff changeset	133 /*
3641167e5b75 codesprint: .cc helper functions do not need tests Rik <octave@nomad.inbox5.com>* parents: 11586 diff changeset	134 ## No test needed for internal helper function.
3641167e5b75 codesprint: .cc helper functions do not need tests Rik <octave@nomad.inbox5.com>* parents: 11586 diff changeset	135 %!assert (1)
3641167e5b75 codesprint: .cc helper functions do not need tests Rik <octave@nomad.inbox5.com>* parents: 11586 diff changeset	136 */

Mercurial > octave

annotate libinterp/corefcn/__pchip_deriv__.cc @ 22322:93b3cdd36854