octave: libinterp/corefcn/__pchip_deriv_

annotate libinterp/corefcn/__pchip_deriv__.cc @ 27923:bd51beb6205e

update formatting of copyright notices * Use <https://octave.org/copyright/> instead of <https://octave.org/COPYRIGHT.html/>. * For consistency with other comments in the Octave sources, use C++-style comments for copyright blocks in C and C++ files. * Use delimiters above and below copyright blocks that are appropriate for the language used in the file. * Eliminate extra spacing inside copyright blocks. * lex.ll (looks_like_copyright): Also allow newlines and carriage returns before the word "Copyright". * scripts/mk-doc.pl (gethelp): Also skip empty comment lines. * bp-table.cc, type.m: Adjust tests.

author	John W. Eaton <jwe@octave.org>
date	Wed, 08 Jan 2020 11:59:41 -0500
parents	1891570abac8
children	0a5b15007766

rev	line source
27923 bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	1 ////////////////////////////////////////////////////////////////////////
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	2 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	3 // Copyright (C) 2002-2020 The Octave Project Developers
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	4 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	5 // See the file COPYRIGHT.md in the top-level directory of this
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	6 // distribution or <https://octave.org/copyright/>.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	7 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	8 // This file is part of Octave.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	9 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	10 // Octave is free software: you can redistribute it and/or modify it
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	11 // under the terms of the GNU General Public License as published by
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	12 // the Free Software Foundation, either version 3 of the License, or
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	13 // (at your option) any later version.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	14 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	15 // Octave is distributed in the hope that it will be useful, but
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	16 // WITHOUT ANY WARRANTY; without even the implied warranty of
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	17 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	18 // GNU General Public License for more details.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	19 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	20 // You should have received a copy of the GNU General Public License
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	21 // along with Octave; see the file COPYING. If not, see
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	22 // <https://www.gnu.org/licenses/>.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	23 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	24 ////////////////////////////////////////////////////////////////////////
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	25
21724 aba2e6293dd8 use "#if ..." consistently instead of "#ifdef" and "#ifndef" John W. Eaton <jwe@octave.org> parents: 21301 diff changeset	26 #if defined (HAVE_CONFIG_H)
21301 40de9f8f23a6 Use '#include "config.h"' rather than <config.h>. Rik <rik@octave.org> parents: 21200 diff changeset	27 # include "config.h"
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	28 #endif
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	29
22322 93b3cdd36854 move most f77 function decls to separate header files John W. Eaton <jwe@octave.org> parents: 22133 diff changeset	30 #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	31
15039 e753177cde93 maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory Rik <rik@octave.org> parents: 14846 diff changeset	32 #include "defun.h"
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	33 #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	34 #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	35 #include "ovl.h"
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	36 #include "utils.h"
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	37 #include "f77-fcn.h"
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	38
6945 6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6549 diff changeset	39 // 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	40 // for piecewise polynomials.
6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe] jwe parents: 6549 diff changeset	41
15039 e753177cde93 maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory Rik <rik@octave.org> parents: 14846 diff changeset	42 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	43 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	44 @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	45 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	46 @end deftypefn */)
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	47 {
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	48 octave_value retval;
21078 49852ff04747 maint: Remove unnecessary declarations of retval. Rik <rik@octave.org> parents: 20940 diff changeset	49 int nargin = args.length ();
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	50
18112 b560bac0fca2 maint: Don't use space between 'args' and '(' when doing indexing. Rik <rik@octave.org> parents: 17787 diff changeset	51 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	52
010e15c7be9a support pchip method in interp2 Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	53 if (nargin >= 2)
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	54 {
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	55 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	56 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	57 FloatColumnVector xvec (args(0).float_vector_value ());
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	58 FloatMatrix ymat (args(1).float_matrix_value ());
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	59
22988 cd33c785e80e put to_f77_int inside octave namespace John W. Eaton <jwe@octave.org> parents: 22981 diff changeset	60 F77_INT nx = octave::to_f77_int (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	61
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 if (nx < 2)
20831 35241c4b696c eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20232 diff changeset	63 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	64
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	65 octave_idx_type nyr = ymat.rows ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	66 octave_idx_type nyc = ymat.columns ();
5838 376e02b2ce70 [project @ 2006-06-01 20:23:53 by jwe] jwe parents: 5837 diff changeset	67
8712 010e15c7be9a support pchip method in interp2 Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	68 if (nx != (rows ? nyc : nyr))
20831 35241c4b696c eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20232 diff changeset	69 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	70
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	71 FloatMatrix dmat (nyr, nyc);
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	72
22957 92c9f8ec4687 fix compilation of __pchip_deriv__.cc John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	73 F77_INT ierr;
23450 855122b993da maint: Wrap tertiary operator in parentheses "(COND ? x : y)". Rik <rik@octave.org> parents: 23220 diff changeset	74 const F77_INT incfd = (rows ? octave::to_f77_int (nyr) : 1);
855122b993da maint: Wrap tertiary operator in parentheses "(COND ? x : y)". Rik <rik@octave.org> parents: 23220 diff changeset	75 volatile const octave_idx_type inc = (rows ? 1 : nyr);
19787 4ae2987c5f66 avoid variable possibly clobbered by longjmp warning from GCC John W. Eaton <jwe@octave.org> parents: 19780 diff changeset	76 volatile octave_idx_type k = 0;
8712 010e15c7be9a support pchip method in interp2 Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	77
19873 5261186f1b00 __pchip_deriv__.cc: Stop variable might be clobbered by 'longjmp' warning. Rik <rik@octave.org> parents: 19787 diff changeset	78 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	79 {
18603 e7369c43ebe0 Fix segmentation fault for interp2 and pchip method (bug #41838). Stefan Mahr <dac922@gmx.de> parents: 17787 diff changeset	80 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	81 ymat.data () + k * inc,
4ae2987c5f66 avoid variable possibly clobbered by longjmp warning from GCC John W. Eaton <jwe@octave.org> parents: 19780 diff changeset	82 dmat.fortran_vec () + k * inc,
22957 92c9f8ec4687 fix compilation of __pchip_deriv__.cc John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	83 incfd, ierr));
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	84
19787 4ae2987c5f66 avoid variable possibly clobbered by longjmp warning from GCC John W. Eaton <jwe@octave.org> parents: 19780 diff changeset	85 k++;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	86
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	87 if (ierr < 0)
20831 35241c4b696c eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20232 diff changeset	88 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	89 }
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	90
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	91 retval = dmat;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	92 }
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	93 else
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	94 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	95 ColumnVector xvec (args(0).vector_value ());
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	96 Matrix ymat (args(1).matrix_value ());
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	97
22988 cd33c785e80e put to_f77_int inside octave namespace John W. Eaton <jwe@octave.org> parents: 22981 diff changeset	98 F77_INT nx = octave::to_f77_int (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	99
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 if (nx < 2)
20831 35241c4b696c eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20232 diff changeset	101 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	102
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	103 octave_idx_type nyr = ymat.rows ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	104 octave_idx_type nyc = ymat.columns ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	105
8712 010e15c7be9a support pchip method in interp2 Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	106 if (nx != (rows ? nyc : nyr))
20831 35241c4b696c eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20232 diff changeset	107 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	108
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	109 Matrix dmat (nyr, nyc);
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	110
22957 92c9f8ec4687 fix compilation of __pchip_deriv__.cc John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	111 F77_INT ierr;
23450 855122b993da maint: Wrap tertiary operator in parentheses "(COND ? x : y)". Rik <rik@octave.org> parents: 23220 diff changeset	112 const F77_INT incfd = (rows ? octave::to_f77_int (nyr) : 1);
855122b993da maint: Wrap tertiary operator in parentheses "(COND ? x : y)". Rik <rik@octave.org> parents: 23220 diff changeset	113 volatile const octave_idx_type inc = (rows ? 1 : nyr);
19787 4ae2987c5f66 avoid variable possibly clobbered by longjmp warning from GCC John W. Eaton <jwe@octave.org> parents: 19780 diff changeset	114 volatile octave_idx_type k = 0;
8712 010e15c7be9a support pchip method in interp2 Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	115
19873 5261186f1b00 __pchip_deriv__.cc: Stop variable might be clobbered by 'longjmp' warning. Rik <rik@octave.org> parents: 19787 diff changeset	116 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	117 {
18603 e7369c43ebe0 Fix segmentation fault for interp2 and pchip method (bug #41838). Stefan Mahr <dac922@gmx.de> parents: 17787 diff changeset	118 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	119 ymat.data () + k * inc,
4ae2987c5f66 avoid variable possibly clobbered by longjmp warning from GCC John W. Eaton <jwe@octave.org> parents: 19780 diff changeset	120 dmat.fortran_vec () + k * inc,
22957 92c9f8ec4687 fix compilation of __pchip_deriv__.cc John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	121 incfd, ierr));
19787 4ae2987c5f66 avoid variable possibly clobbered by longjmp warning from GCC John W. Eaton <jwe@octave.org> parents: 19780 diff changeset	122 k++;
5838 376e02b2ce70 [project @ 2006-06-01 20:23:53 by jwe] jwe parents: 5837 diff changeset	123
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	124 if (ierr < 0)
20831 35241c4b696c eliminate return statements after calls to error John W. Eaton <jwe@octave.org> parents: 20232 diff changeset	125 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	126 }
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	127
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	128 retval = dmat;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 8713 diff changeset	129 }
5837 55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	130 }
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	131
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	132 return retval;
55404f3b0da1 [project @ 2006-06-01 19:05:31 by jwe] jwe parents: diff changeset	133 }
12805 3641167e5b75 codesprint: .cc helper functions do not need tests Rik <octave@nomad.inbox5.com>* parents: 11586 diff changeset	134
3641167e5b75 codesprint: .cc helper functions do not need tests Rik <octave@nomad.inbox5.com>* parents: 11586 diff changeset	135 /*
3641167e5b75 codesprint: .cc helper functions do not need tests Rik <octave@nomad.inbox5.com>* parents: 11586 diff changeset	136 ## 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	137 %!assert (1)
3641167e5b75 codesprint: .cc helper functions do not need tests Rik <octave@nomad.inbox5.com>* parents: 11586 diff changeset	138 */

Mercurial > octave

annotate libinterp/corefcn/__pchip_deriv__.cc @ 27923:bd51beb6205e