Mercurial > octave

comparison liboctave/numeric/lo-specfun.cc @ 22197:e43d83253e28

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
refill multi-line macro definitions Use the Emacs C++ mode style for line continuation markers in multi-line macro definitions. * make_int.cc, __dsearchn__.cc, __magick_read__.cc, besselj.cc, bitfcns.cc, bsxfun.cc, cellfun.cc, data.cc, defun-dld.h, defun-int.h, defun.h, det.cc, error.h, find.cc, gcd.cc, graphics.cc, interpreter.h, jit-ir.h, jit-typeinfo.h, lookup.cc, ls-mat5.cc, max.cc, mexproto.h, mxarray.in.h, oct-stream.cc, ordschur.cc, pr-output.cc, profiler.h, psi.cc, regexp.cc, sparse-xdiv.cc, sparse-xpow.cc, tril.cc, txt-eng.h, utils.cc, variables.cc, variables.h, xdiv.cc, xpow.cc, __glpk__.cc, ov-base.cc, ov-base.h, ov-cell.cc, ov-ch-mat.cc, ov-classdef.cc, ov-complex.cc, ov-cx-mat.cc, ov-cx-sparse.cc, ov-float.cc, ov-float.h, ov-flt-complex.cc, ov-flt-cx-mat.cc, ov-flt-re-mat.cc, ov-int-traits.h, ov-lazy-idx.h, ov-perm.cc, ov-re-mat.cc, ov-re-sparse.cc, ov-scalar.cc, ov-scalar.h, ov-str-mat.cc, ov-type-conv.h, ov.cc, ov.h, op-class.cc, op-int-conv.cc, op-int.h, op-str-str.cc, ops.h, lex.ll, Array.cc, CMatrix.cc, CSparse.cc, MArray.cc, MArray.h, MDiagArray2.cc, MDiagArray2.h, MSparse.h, Sparse.cc, dMatrix.cc, dSparse.cc, fCMatrix.cc, fMatrix.cc, idx-vector.cc, f77-fcn.h, quit.h, bsxfun-decl.h, bsxfun-defs.cc, lo-specfun.cc, oct-convn.cc, oct-convn.h, oct-norm.cc, oct-norm.h, oct-rand.cc, Sparse-op-decls.h, Sparse-op-defs.h, mx-inlines.cc, mx-op-decl.h, mx-op-defs.h, mach-info.cc, oct-group.cc, oct-passwd.cc, oct-syscalls.cc, oct-time.cc, data-conv.cc, kpse.cc, lo-ieee.h, lo-macros.h, oct-cmplx.h, oct-glob.cc, oct-inttypes.cc, oct-inttypes.h, oct-locbuf.h, oct-sparse.h, url-transfer.cc, oct-conf-post.in.h, shared-fcns.h: Refill macro definitions.
author John W. Eaton <jwe@octave.org>
date Mon, 01 Aug 2016 12:40:18 -0400
parents 407c66ae1e20
children edcad35c364f
comparison
equal deleted inserted replaced
22196:dd992fd74fce 22197:e43d83253e28
1400 retval(i,j) = f (x(i), alpha(j), (scaled ? 2 : 1), ierr(i,j)); 1400 retval(i,j) = f (x(i), alpha(j), (scaled ? 2 : 1), ierr(i,j));
1401 1401
1402 return retval; 1402 return retval;
1403 } 1403 }
1404 1404
1405 #define SS_BESSEL(name, fcn) \ 1405 #define SS_BESSEL(name, fcn) \
1406 Complex \ 1406 Complex \
1407 name (double alpha, const Complex& x, bool scaled, octave_idx_type& ierr) \ 1407 name (double alpha, const Complex& x, bool scaled, octave_idx_type& ierr) \
1408 { \ 1408 { \
1409 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \ 1409 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \
1410 } 1410 }
1411 1411
1412 #define SM_BESSEL(name, fcn) \ 1412 #define SM_BESSEL(name, fcn) \
1413 ComplexMatrix \ 1413 ComplexMatrix \
1414 name (double alpha, const ComplexMatrix& x, bool scaled, \ 1414 name (double alpha, const ComplexMatrix& x, bool scaled, \
1415 Array<octave_idx_type>& ierr) \ 1415 Array<octave_idx_type>& ierr) \
1416 { \ 1416 { \
1417 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \ 1417 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \
1418 } 1418 }
1419 1419
1420 #define MS_BESSEL(name, fcn) \ 1420 #define MS_BESSEL(name, fcn) \
1421 ComplexMatrix \ 1421 ComplexMatrix \
1422 name (const Matrix& alpha, const Complex& x, bool scaled, \ 1422 name (const Matrix& alpha, const Complex& x, bool scaled, \
1423 Array<octave_idx_type>& ierr) \ 1423 Array<octave_idx_type>& ierr) \
1424 { \ 1424 { \
1425 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \ 1425 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \
1426 } 1426 }
1427 1427
1428 #define MM_BESSEL(name, fcn) \ 1428 #define MM_BESSEL(name, fcn) \
1429 ComplexMatrix \ 1429 ComplexMatrix \
1430 name (const Matrix& alpha, const ComplexMatrix& x, bool scaled, \ 1430 name (const Matrix& alpha, const ComplexMatrix& x, bool scaled, \
1431 Array<octave_idx_type>& ierr) \ 1431 Array<octave_idx_type>& ierr) \
1432 { \ 1432 { \
1433 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \ 1433 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \
1434 } 1434 }
1435 1435
1436 #define SN_BESSEL(name, fcn) \ 1436 #define SN_BESSEL(name, fcn) \
1437 ComplexNDArray \ 1437 ComplexNDArray \
1438 name (double alpha, const ComplexNDArray& x, bool scaled, \ 1438 name (double alpha, const ComplexNDArray& x, bool scaled, \
1439 Array<octave_idx_type>& ierr) \ 1439 Array<octave_idx_type>& ierr) \
1440 { \ 1440 { \
1441 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \ 1441 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \
1442 } 1442 }
1443 1443
1444 #define NS_BESSEL(name, fcn) \ 1444 #define NS_BESSEL(name, fcn) \
1445 ComplexNDArray \ 1445 ComplexNDArray \
1446 name (const NDArray& alpha, const Complex& x, bool scaled, \ 1446 name (const NDArray& alpha, const Complex& x, bool scaled, \
1447 Array<octave_idx_type>& ierr) \ 1447 Array<octave_idx_type>& ierr) \
1448 { \ 1448 { \
1449 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \ 1449 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \
1450 } 1450 }
1451 1451
1452 #define NN_BESSEL(name, fcn) \ 1452 #define NN_BESSEL(name, fcn) \
1453 ComplexNDArray \ 1453 ComplexNDArray \
1454 name (const NDArray& alpha, const ComplexNDArray& x, bool scaled, \ 1454 name (const NDArray& alpha, const ComplexNDArray& x, bool scaled, \
1455 Array<octave_idx_type>& ierr) \ 1455 Array<octave_idx_type>& ierr) \
1456 { \ 1456 { \
1457 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \ 1457 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \
1458 } 1458 }
1459 1459
1460 #define RC_BESSEL(name, fcn) \ 1460 #define RC_BESSEL(name, fcn) \
1461 ComplexMatrix \ 1461 ComplexMatrix \
1462 name (const RowVector& alpha, const ComplexColumnVector& x, bool scaled, \ 1462 name (const RowVector& alpha, const ComplexColumnVector& x, bool scaled, \
1463 Array<octave_idx_type>& ierr) \ 1463 Array<octave_idx_type>& ierr) \
1464 { \ 1464 { \
1465 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \ 1465 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \
1466 } 1466 }
1467 1467
1468 #define ALL_BESSEL(name, fcn) \ 1468 #define ALL_BESSEL(name, fcn) \
1469 SS_BESSEL (name, fcn) \ 1469 SS_BESSEL (name, fcn) \
1470 SM_BESSEL (name, fcn) \ 1470 SM_BESSEL (name, fcn) \
1471 MS_BESSEL (name, fcn) \ 1471 MS_BESSEL (name, fcn) \
1472 MM_BESSEL (name, fcn) \ 1472 MM_BESSEL (name, fcn) \
1473 SN_BESSEL (name, fcn) \ 1473 SN_BESSEL (name, fcn) \
1474 NS_BESSEL (name, fcn) \ 1474 NS_BESSEL (name, fcn) \
1475 NN_BESSEL (name, fcn) \ 1475 NN_BESSEL (name, fcn) \
1476 RC_BESSEL (name, fcn) 1476 RC_BESSEL (name, fcn)
1477 1477
1478 ALL_BESSEL (besselj, zbesj) 1478 ALL_BESSEL (besselj, zbesj)
1479 ALL_BESSEL (bessely, zbesy) 1479 ALL_BESSEL (bessely, zbesy)
1480 ALL_BESSEL (besseli, zbesi) 1480 ALL_BESSEL (besseli, zbesi)
2013 retval(i,j) = f (x(i), alpha(j), (scaled ? 2 : 1), ierr(i,j)); 2013 retval(i,j) = f (x(i), alpha(j), (scaled ? 2 : 1), ierr(i,j));
2014 2014
2015 return retval; 2015 return retval;
2016 } 2016 }
2017 2017
2018 #define SS_BESSEL(name, fcn) \ 2018 #define SS_BESSEL(name, fcn) \
2019 FloatComplex \ 2019 FloatComplex \
2020 name (float alpha, const FloatComplex& x, bool scaled, octave_idx_type& ierr) \ 2020 name (float alpha, const FloatComplex& x, bool scaled, octave_idx_type& ierr) \
2021 { \ 2021 { \
2022 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \ 2022 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \
2023 } 2023 }
2024 2024
2025 #define SM_BESSEL(name, fcn) \ 2025 #define SM_BESSEL(name, fcn) \
2026 FloatComplexMatrix \ 2026 FloatComplexMatrix \
2027 name (float alpha, const FloatComplexMatrix& x, bool scaled, \ 2027 name (float alpha, const FloatComplexMatrix& x, bool scaled, \
2028 Array<octave_idx_type>& ierr) \ 2028 Array<octave_idx_type>& ierr) \
2029 { \ 2029 { \
2030 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \ 2030 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \
2031 } 2031 }
2032 2032
2033 #define MS_BESSEL(name, fcn) \ 2033 #define MS_BESSEL(name, fcn) \
2034 FloatComplexMatrix \ 2034 FloatComplexMatrix \
2035 name (const FloatMatrix& alpha, const FloatComplex& x, bool scaled, \ 2035 name (const FloatMatrix& alpha, const FloatComplex& x, bool scaled, \
2036 Array<octave_idx_type>& ierr) \ 2036 Array<octave_idx_type>& ierr) \
2037 { \ 2037 { \
2038 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \ 2038 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \
2039 } 2039 }
2040 2040
2041 #define MM_BESSEL(name, fcn) \ 2041 #define MM_BESSEL(name, fcn) \
2042 FloatComplexMatrix \ 2042 FloatComplexMatrix \
2043 name (const FloatMatrix& alpha, const FloatComplexMatrix& x, bool scaled, \ 2043 name (const FloatMatrix& alpha, const FloatComplexMatrix& x, bool scaled, \
2044 Array<octave_idx_type>& ierr) \ 2044 Array<octave_idx_type>& ierr) \
2045 { \ 2045 { \
2046 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \ 2046 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \
2047 } 2047 }
2048 2048
2049 #define SN_BESSEL(name, fcn) \ 2049 #define SN_BESSEL(name, fcn) \
2050 FloatComplexNDArray \ 2050 FloatComplexNDArray \
2051 name (float alpha, const FloatComplexNDArray& x, bool scaled, \ 2051 name (float alpha, const FloatComplexNDArray& x, bool scaled, \
2052 Array<octave_idx_type>& ierr) \ 2052 Array<octave_idx_type>& ierr) \
2053 { \ 2053 { \
2054 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \ 2054 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \
2055 } 2055 }
2056 2056
2057 #define NS_BESSEL(name, fcn) \ 2057 #define NS_BESSEL(name, fcn) \
2058 FloatComplexNDArray \ 2058 FloatComplexNDArray \
2059 name (const FloatNDArray& alpha, const FloatComplex& x, bool scaled, \ 2059 name (const FloatNDArray& alpha, const FloatComplex& x, bool scaled, \
2060 Array<octave_idx_type>& ierr) \ 2060 Array<octave_idx_type>& ierr) \
2061 { \ 2061 { \
2062 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \ 2062 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \
2063 } 2063 }
2064 2064
2065 #define NN_BESSEL(name, fcn) \ 2065 #define NN_BESSEL(name, fcn) \
2066 FloatComplexNDArray \ 2066 FloatComplexNDArray \
2067 name (const FloatNDArray& alpha, const FloatComplexNDArray& x, bool scaled, \ 2067 name (const FloatNDArray& alpha, const FloatComplexNDArray& x, bool scaled, \
2068 Array<octave_idx_type>& ierr) \ 2068 Array<octave_idx_type>& ierr) \
2069 { \ 2069 { \
2070 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \ 2070 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \
2071 } 2071 }
2072 2072
2073 #define RC_BESSEL(name, fcn) \ 2073 #define RC_BESSEL(name, fcn) \
2074 FloatComplexMatrix \ 2074 FloatComplexMatrix \
2075 name (const FloatRowVector& alpha, const FloatComplexColumnVector& x, bool scaled, \ 2075 name (const FloatRowVector& alpha, const FloatComplexColumnVector& x, bool scaled, \
2076 Array<octave_idx_type>& ierr) \ 2076 Array<octave_idx_type>& ierr) \
2077 { \ 2077 { \
2078 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \ 2078 return do_bessel (fcn, #name, alpha, x, scaled, ierr); \
2079 } 2079 }
2080 2080
2081 #define ALL_BESSEL(name, fcn) \ 2081 #define ALL_BESSEL(name, fcn) \
2082 SS_BESSEL (name, fcn) \ 2082 SS_BESSEL (name, fcn) \
2083 SM_BESSEL (name, fcn) \ 2083 SM_BESSEL (name, fcn) \
2084 MS_BESSEL (name, fcn) \ 2084 MS_BESSEL (name, fcn) \
2085 MM_BESSEL (name, fcn) \ 2085 MM_BESSEL (name, fcn) \
2086 SN_BESSEL (name, fcn) \ 2086 SN_BESSEL (name, fcn) \
2087 NS_BESSEL (name, fcn) \ 2087 NS_BESSEL (name, fcn) \
2088 NN_BESSEL (name, fcn) \ 2088 NN_BESSEL (name, fcn) \
2089 RC_BESSEL (name, fcn) 2089 RC_BESSEL (name, fcn)
2090 2090
2091 ALL_BESSEL (besselj, cbesj) 2091 ALL_BESSEL (besselj, cbesj)
2092 ALL_BESSEL (bessely, cbesy) 2092 ALL_BESSEL (bessely, cbesy)
2093 ALL_BESSEL (besseli, cbesi) 2093 ALL_BESSEL (besseli, cbesi)