octave: libinterp/corefcn/stream-euler.cc annotate

annotate libinterp/corefcn/stream-euler.cc @ 31605:e88a07dec498 stable

maint: Use macros to begin/end C++ namespaces. * oct-conf-post-public.in.h: Define two macros (OCTAVE_BEGIN_NAMESPACE, OCTAVE_END_NAMESPACE) that can be used to start/end a namespace. * mk-opts.pl, build-env.h, build-env.in.cc, __betainc__.cc, __contourc__.cc, __dsearchn__.cc, __eigs__.cc, __expint__.cc, __ftp__.cc, __gammainc__.cc, __ichol__.cc, __ilu__.cc, __isprimelarge__.cc, __lin_interpn__.cc, __magick_read__.cc, __pchip_deriv__.cc, __qp__.cc, amd.cc, auto-shlib.cc, auto-shlib.h, balance.cc, base-text-renderer.cc, base-text-renderer.h, besselj.cc, bitfcns.cc, bsxfun.cc, c-file-ptr-stream.cc, c-file-ptr-stream.h, call-stack.cc, call-stack.h, ccolamd.cc, cellfun.cc, chol.cc, colamd.cc, colloc.cc, conv2.cc, daspk.cc, dasrt.cc, dassl.cc, data.cc, data.h, debug.cc, defaults.cc, defaults.h, defun-int.h, defun.cc, det.cc, dirfns.cc, display.cc, display.h, dlmread.cc, dmperm.cc, dot.cc, dynamic-ld.cc, dynamic-ld.h, eig.cc, ellipj.cc, environment.cc, environment.h, error.cc, error.h, errwarn.h, event-manager.cc, event-manager.h, event-queue.cc, event-queue.h, fcn-info.cc, fcn-info.h, fft.cc, fft2.cc, fftn.cc, file-io.cc, filter.cc, find.cc, ft-text-renderer.cc, ft-text-renderer.h, gcd.cc, getgrent.cc, getpwent.cc, getrusage.cc, givens.cc, gl-render.cc, gl-render.h, gl2ps-print.cc, gl2ps-print.h, graphics-toolkit.cc, graphics-toolkit.h, graphics.cc, graphics.in.h, gsvd.cc, gtk-manager.cc, gtk-manager.h, hash.cc, help.cc, help.h, hess.cc, hex2num.cc, hook-fcn.cc, hook-fcn.h, input.cc, input.h, interpreter-private.cc, interpreter-private.h, interpreter.cc, interpreter.h, inv.cc, jsondecode.cc, jsonencode.cc, kron.cc, latex-text-renderer.cc, latex-text-renderer.h, load-path.cc, load-path.h, load-save.cc, load-save.h, lookup.cc, ls-ascii-helper.cc, ls-ascii-helper.h, ls-oct-text.cc, ls-utils.cc, ls-utils.h, lsode.cc, lu.cc, mappers.cc, matrix_type.cc, max.cc, mex-private.h, mex.cc, mgorth.cc, nproc.cc, oct-fstrm.cc, oct-fstrm.h, oct-hdf5-types.cc, oct-hdf5-types.h, oct-hist.cc, oct-hist.h, oct-iostrm.cc, oct-iostrm.h, oct-opengl.h, oct-prcstrm.cc, oct-prcstrm.h, oct-procbuf.cc, oct-procbuf.h, oct-process.cc, oct-process.h, oct-stdstrm.h, oct-stream.cc, oct-stream.h, oct-strstrm.cc, oct-strstrm.h, oct-tex-lexer.in.ll, oct-tex-parser.yy, ordqz.cc, ordschur.cc, pager.cc, pager.h, pinv.cc, pow2.cc, pr-flt-fmt.cc, pr-output.cc, procstream.cc, procstream.h, psi.cc, qr.cc, quad.cc, quadcc.cc, qz.cc, rand.cc, rcond.cc, regexp.cc, schur.cc, settings.cc, settings.h, sighandlers.cc, sighandlers.h, sparse-xdiv.cc, sparse-xdiv.h, sparse-xpow.cc, sparse-xpow.h, sparse.cc, spparms.cc, sqrtm.cc, stack-frame.cc, stack-frame.h, stream-euler.cc, strfind.cc, strfns.cc, sub2ind.cc, svd.cc, sylvester.cc, symbfact.cc, syminfo.cc, syminfo.h, symrcm.cc, symrec.cc, symrec.h, symscope.cc, symscope.h, symtab.cc, symtab.h, syscalls.cc, sysdep.cc, sysdep.h, text-engine.cc, text-engine.h, text-renderer.cc, text-renderer.h, time.cc, toplev.cc, tril.cc, tsearch.cc, typecast.cc, url-handle-manager.cc, url-handle-manager.h, urlwrite.cc, utils.cc, utils.h, variables.cc, variables.h, xdiv.cc, xdiv.h, xnorm.cc, xnorm.h, xpow.cc, xpow.h, __delaunayn__.cc, __fltk_uigetfile__.cc, __glpk__.cc, __init_fltk__.cc, __init_gnuplot__.cc, __ode15__.cc, __voronoi__.cc, audiodevinfo.cc, audioread.cc, convhulln.cc, fftw.cc, gzip.cc, mk-build-env-features.sh, mk-builtins.pl, cdef-class.cc, cdef-class.h, cdef-fwd.h, cdef-manager.cc, cdef-manager.h, cdef-method.cc, cdef-method.h, cdef-object.cc, cdef-object.h, cdef-package.cc, cdef-package.h, cdef-property.cc, cdef-property.h, cdef-utils.cc, cdef-utils.h, ov-base.cc, ov-base.h, ov-bool-mat.cc, ov-builtin.h, ov-cell.cc, ov-class.cc, ov-class.h, ov-classdef.cc, ov-classdef.h, ov-complex.cc, ov-fcn-handle.cc, ov-fcn-handle.h, ov-fcn.h, ov-java.cc, ov-java.h, ov-mex-fcn.h, ov-null-mat.cc, ov-oncleanup.cc, ov-struct.cc, ov-typeinfo.cc, ov-typeinfo.h, ov-usr-fcn.cc, ov-usr-fcn.h, ov.cc, ov.h, octave.cc, octave.h, mk-ops.sh, op-b-b.cc, op-b-bm.cc, op-b-sbm.cc, op-bm-b.cc, op-bm-bm.cc, op-bm-sbm.cc, op-cdm-cdm.cc, op-cell.cc, op-chm.cc, op-class.cc, op-cm-cm.cc, op-cm-cs.cc, op-cm-m.cc, op-cm-s.cc, op-cm-scm.cc, op-cm-sm.cc, op-cs-cm.cc, op-cs-cs.cc, op-cs-m.cc, op-cs-s.cc, op-cs-scm.cc, op-cs-sm.cc, op-dm-dm.cc, op-dm-scm.cc, op-dm-sm.cc, op-dm-template.cc, op-dms-template.cc, op-fcdm-fcdm.cc, op-fcm-fcm.cc, op-fcm-fcs.cc, op-fcm-fm.cc, op-fcm-fs.cc, op-fcn.cc, op-fcs-fcm.cc, op-fcs-fcs.cc, op-fcs-fm.cc, op-fcs-fs.cc, op-fdm-fdm.cc, op-fm-fcm.cc, op-fm-fcs.cc, op-fm-fm.cc, op-fm-fs.cc, op-fs-fcm.cc, op-fs-fcs.cc, op-fs-fm.cc, op-fs-fs.cc, op-i16-i16.cc, op-i32-i32.cc, op-i64-i64.cc, op-i8-i8.cc, op-int-concat.cc, op-m-cm.cc, op-m-cs.cc, op-m-m.cc, op-m-s.cc, op-m-scm.cc, op-m-sm.cc, op-mi.cc, op-pm-pm.cc, op-pm-scm.cc, op-pm-sm.cc, op-pm-template.cc, op-range.cc, op-s-cm.cc, op-s-cs.cc, op-s-m.cc, op-s-s.cc, op-s-scm.cc, op-s-sm.cc, op-sbm-b.cc, op-sbm-bm.cc, op-sbm-sbm.cc, op-scm-cm.cc, op-scm-cs.cc, op-scm-m.cc, op-scm-s.cc, op-scm-scm.cc, op-scm-sm.cc, op-sm-cm.cc, op-sm-cs.cc, op-sm-m.cc, op-sm-s.cc, op-sm-scm.cc, op-sm-sm.cc, op-str-m.cc, op-str-s.cc, op-str-str.cc, op-struct.cc, op-ui16-ui16.cc, op-ui32-ui32.cc, op-ui64-ui64.cc, op-ui8-ui8.cc, ops.h, anon-fcn-validator.cc, anon-fcn-validator.h, bp-table.cc, bp-table.h, comment-list.cc, comment-list.h, filepos.h, lex.h, lex.ll, oct-lvalue.cc, oct-lvalue.h, oct-parse.yy, parse.h, profiler.cc, profiler.h, pt-anon-scopes.cc, pt-anon-scopes.h, pt-arg-list.cc, pt-arg-list.h, pt-args-block.cc, pt-args-block.h, pt-array-list.cc, pt-array-list.h, pt-assign.cc, pt-assign.h, pt-binop.cc, pt-binop.h, pt-bp.cc, pt-bp.h, pt-cbinop.cc, pt-cbinop.h, pt-cell.cc, pt-cell.h, pt-check.cc, pt-check.h, pt-classdef.cc, pt-classdef.h, pt-cmd.h, pt-colon.cc, pt-colon.h, pt-const.cc, pt-const.h, pt-decl.cc, pt-decl.h, pt-eval.cc, pt-eval.h, pt-except.cc, pt-except.h, pt-exp.cc, pt-exp.h, pt-fcn-handle.cc, pt-fcn-handle.h, pt-id.cc, pt-id.h, pt-idx.cc, pt-idx.h, pt-jump.h, pt-loop.cc, pt-loop.h, pt-mat.cc, pt-mat.h, pt-misc.cc, pt-misc.h, pt-pr-code.cc, pt-pr-code.h, pt-select.cc, pt-select.h, pt-spmd.cc, pt-spmd.h, pt-stmt.cc, pt-stmt.h, pt-tm-const.cc, pt-tm-const.h, pt-unop.cc, pt-unop.h, pt-vm-eval.cc, pt-walk.cc, pt-walk.h, pt.cc, pt.h, token.cc, token.h, Range.cc, Range.h, idx-vector.cc, idx-vector.h, range-fwd.h, CollocWt.cc, CollocWt.h, aepbalance.cc, aepbalance.h, chol.cc, chol.h, gepbalance.cc, gepbalance.h, gsvd.cc, gsvd.h, hess.cc, hess.h, lo-mappers.cc, lo-mappers.h, lo-specfun.cc, lo-specfun.h, lu.cc, lu.h, oct-convn.cc, oct-convn.h, oct-fftw.cc, oct-fftw.h, oct-norm.cc, oct-norm.h, oct-rand.cc, oct-rand.h, oct-spparms.cc, oct-spparms.h, qr.cc, qr.h, qrp.cc, qrp.h, randgamma.cc, randgamma.h, randmtzig.cc, randmtzig.h, randpoisson.cc, randpoisson.h, schur.cc, schur.h, sparse-chol.cc, sparse-chol.h, sparse-lu.cc, sparse-lu.h, sparse-qr.cc, sparse-qr.h, svd.cc, svd.h, child-list.cc, child-list.h, dir-ops.cc, dir-ops.h, file-ops.cc, file-ops.h, file-stat.cc, file-stat.h, lo-sysdep.cc, lo-sysdep.h, lo-sysinfo.cc, lo-sysinfo.h, mach-info.cc, mach-info.h, oct-env.cc, oct-env.h, oct-group.cc, oct-group.h, oct-password.cc, oct-password.h, oct-syscalls.cc, oct-syscalls.h, oct-time.cc, oct-time.h, oct-uname.cc, oct-uname.h, action-container.cc, action-container.h, base-list.h, cmd-edit.cc, cmd-edit.h, cmd-hist.cc, cmd-hist.h, f77-fcn.h, file-info.cc, file-info.h, lo-array-errwarn.cc, lo-array-errwarn.h, lo-hash.cc, lo-hash.h, lo-ieee.h, lo-regexp.cc, lo-regexp.h, lo-utils.cc, lo-utils.h, oct-base64.cc, oct-base64.h, oct-glob.cc, oct-glob.h, oct-inttypes.h, oct-mutex.cc, oct-mutex.h, oct-refcount.h, oct-shlib.cc, oct-shlib.h, oct-sparse.cc, oct-sparse.h, oct-string.h, octave-preserve-stream-state.h, pathsearch.cc, pathsearch.h, quit.cc, quit.h, unwind-prot.cc, unwind-prot.h, url-transfer.cc, url-transfer.h : Use new macros to begin/end C++ namespaces.

author	Rik <rik@octave.org>
date	Thu, 01 Dec 2022 14:23:45 -0800
parents	a3d8915b9f87
children	aac27ad79be6

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 //
30564 796f54d4ddbf update Octave Project Developers copyright for the new year John W. Eaton <jwe@octave.org> parents: 30043 diff changeset	3 // Copyright (C) 2019-2022 The Octave Project Developers
27923 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 ////////////////////////////////////////////////////////////////////////
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	25
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	26 /*
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	27
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	28 References:
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	29
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	30 @article{
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	31 title = {Particle Tracing Algorithms for 3D Curvilinear Grids},
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	32 year = {2000},
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	33 author = {Nielson, Gregory and Uller, H. and Sadarjoen, I. and Walsum, Theo and Hin, Andrea and Post, Frits}
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	34 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	35
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	36 @article{
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	37 title = {Sources of error in the graphical analysis of CFD results},
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	38 publisher = {Journal of Scientific Computing},
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	39 year = {1988},
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	40 volume = {3},
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	41 number = {2},
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	42 pages = {149--164},
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	43 author = {Buning, Pieter G.},
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	44 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	45
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	46 */
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	47
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	48 #if defined (HAVE_CONFIG_H)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	49 # include "config.h"
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	50 #endif
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	51
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	52 #include "defun.h"
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	53 #include "error.h"
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	54 #include "ovl.h"
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	55
31605 e88a07dec498 maint: Use macros to begin/end C++ namespaces. Rik <rik@octave.org> parents: 31594 diff changeset	56 OCTAVE_BEGIN_NAMESPACE(octave)
29958 32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace John W. Eaton <jwe@octave.org> parents: 29358 diff changeset	57
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	58 // Coordinates of a point in C-Space (unit square mesh)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	59
30043 b8ffd22ee32c maint: use C++ style struct declaration instead of typedef struct in stream-euler.cc Rik <rik@octave.org> parents: 29961 diff changeset	60 struct Vector2
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	61 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	62 double x, y;
30043 b8ffd22ee32c maint: use C++ style struct declaration instead of typedef struct in stream-euler.cc Rik <rik@octave.org> parents: 29961 diff changeset	63 };
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	64
27810 f2b89a2e20b6 overhaul streamXXX.m series of functions. Rik <rik@octave.org> parents: 27808 diff changeset	65 // The integer value and the fractional value from a point in C-Space.
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	66 // Equivalent to the cell index the point is located in and the local
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	67 // coordinates of the point in the cell.
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	68
30043 b8ffd22ee32c maint: use C++ style struct declaration instead of typedef struct in stream-euler.cc Rik <rik@octave.org> parents: 29961 diff changeset	69 struct Cell2
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	70 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	71 double fcx, fcy;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	72 signed long idx, idy;
30043 b8ffd22ee32c maint: use C++ style struct declaration instead of typedef struct in stream-euler.cc Rik <rik@octave.org> parents: 29961 diff changeset	73 };
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	74
30043 b8ffd22ee32c maint: use C++ style struct declaration instead of typedef struct in stream-euler.cc Rik <rik@octave.org> parents: 29961 diff changeset	75 struct Vector3
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	76 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	77 double x, y, z;
30043 b8ffd22ee32c maint: use C++ style struct declaration instead of typedef struct in stream-euler.cc Rik <rik@octave.org> parents: 29961 diff changeset	78 };
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	79
30043 b8ffd22ee32c maint: use C++ style struct declaration instead of typedef struct in stream-euler.cc Rik <rik@octave.org> parents: 29961 diff changeset	80 struct Cell3
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	81 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	82 double fcx, fcy, fcz;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	83 signed long idx, idy, idz;
30043 b8ffd22ee32c maint: use C++ style struct declaration instead of typedef struct in stream-euler.cc Rik <rik@octave.org> parents: 29961 diff changeset	84 };
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	85
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	86 static inline void
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	87 number_to_fractional (signed long id, double fc, const double u)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	88 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	89 *id = floor (u);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	90 fc = u - id;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	91 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	92
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	93 static inline octave_idx_type
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	94 handle_border_index (const octave_idx_type id, const octave_idx_type N)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	95 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	96 return (id < N - 1 ? id : N - 2);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	97 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	98
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	99 static inline void
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	100 handle_border (octave_idx_type id2, double fc2, const octave_idx_type id1,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	101 const double fc1, const octave_idx_type N)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	102 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	103 if (id1 < N - 1)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	104 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	105 *id2 = id1;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	106 *fc2 = fc1;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	107 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	108 else
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	109 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	110 *id2 = N - 2;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	111 *fc2 = 1.0;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	112 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	113 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	114
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	115 static inline double
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	116 bilinear (const double u11, const double u21, const double u12,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	117 const double u22, const double x, const double y)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	118 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	119 return (u11 * (1-x) * (1-y) +
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	120 u21 * x * (1-y) +
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	121 u12 * (1-x) * y +
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	122 u22 * x * y);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	123 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	124
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	125 static inline Cell2
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	126 vector_to_cell2d (const Vector2 X)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	127 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	128 Cell2 Z;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	129
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	130 number_to_fractional (&Z.idx, &Z.fcx, X.x);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	131 number_to_fractional (&Z.idy, &Z.fcy, X.y);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	132
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	133 return (Z);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	134 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	135
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	136 static inline bool
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	137 is_in_definition_set2d (const Cell2 X, const octave_idx_type cols,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	138 const octave_idx_type rows)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	139 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	140 return ( (((X.idx >= 0) && (X.idx < cols-1)) \|\|
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	141 ((X.idx == cols-1) && (X.fcx == 0.0))) &&
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	142 (((X.idy >= 0) && (X.idy < rows-1)) \|\|
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	143 ((X.idy == rows-1) && (X.fcy == 0.0))) );
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	144 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	145
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	146 static inline Vector2
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	147 add2d (const Cell2 X, const Vector2 Y)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	148 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	149 Vector2 Z = {X.idx + X.fcx + Y.x,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	150 X.idy + X.fcy + Y.y};
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	151
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	152 return (Z);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	153 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	154
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	155 static inline Vector2
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	156 vector_interpolation2d (const Cell2 X, const Matrix& u, const Matrix& v,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	157 const octave_idx_type cols, const octave_idx_type rows)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	158 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	159 Vector2 V;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	160 double fcx, fcy;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	161 octave_idx_type idx, idy;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	162
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	163 handle_border (&idx, &fcx, X.idx, X.fcx, cols);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	164 handle_border (&idy, &fcy, X.idy, X.fcy, rows);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	165
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	166 V.x = bilinear (u(idy, idx), u(idy, idx+1), u(idy+1, idx),
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	167 u(idy+1, idx+1), fcx, fcy);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	168 V.y = bilinear (v(idy, idx), v(idy, idx+1), v(idy+1, idx),
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	169 v(idy+1, idx+1), fcx, fcy);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	170
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	171 return (V);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	172 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	173
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	174 // Apply the Jacobian matrix on the vector V.
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	175 // The step vector length is set to h.
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	176
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	177 static inline Vector2
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	178 calculate_step_vector2d (const Cell2 X, const Vector2 V,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	179 const RowVector& tx, const RowVector& ty,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	180 const octave_idx_type cols, const octave_idx_type rows,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	181 const double h)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	182 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	183 Vector2 S;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	184
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	185 const octave_idx_type idx = handle_border_index (X.idx, cols);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	186 const octave_idx_type idy = handle_border_index (X.idy, rows);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	187
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	188 const double x = V.x * tx(idx);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	189 const double y = V.y * ty(idy);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	190 const double n = 1.0 / sqrt (xx + yy);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	191 S.x = h * n * x;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	192 S.y = h * n * y;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	193
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	194 return (S);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	195 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	196
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	197 static inline bool
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	198 is_singular2d (const Vector2 V)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	199 {
29961 7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more John W. Eaton <jwe@octave.org> parents: 29958 diff changeset	200 return ((math::isnan (V.x) \|\| math::isnan (V.y)) \|\|
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	201 ((V.x == 0) && (V.y == 0)));
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	202 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	203
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	204 static void
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	205 euler2d (const octave_idx_type cols, const octave_idx_type rows,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	206 const Matrix& u, const Matrix& v,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	207 const RowVector& tx, const RowVector& ty,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	208 const double zeta, const double xi,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	209 const double h, const octave_idx_type maxnverts,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	210 Matrix& buffer, octave_idx_type *nverts)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	211 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	212 Vector2 V0, V1, S0, X1, Xnxt, S1;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	213 const Vector2 X0 = {zeta, xi};
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	214 Cell2 X0f, X1f;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	215
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	216 octave_idx_type i = 0;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	217
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	218 buffer(i, 0) = X0.x;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	219 buffer(i, 1) = X0.y;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	220
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	221 X0f = vector_to_cell2d (X0);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	222 while (true)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	223 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	224 if (! is_in_definition_set2d (X0f, cols, rows))
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	225 break;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	226
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	227 V0 = vector_interpolation2d (X0f, u, v, cols, rows);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	228 if (is_singular2d (V0))
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	229 break;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	230
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	231 S0 = calculate_step_vector2d (X0f, V0, tx, ty, cols, rows, h);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	232
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	233 X1 = add2d (X0f, S0);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	234 X1f = vector_to_cell2d (X1);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	235 if (! is_in_definition_set2d (X1f, cols, rows))
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	236 break;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	237
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	238 V1 = vector_interpolation2d (X1f, u, v, cols, rows);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	239 if (is_singular2d (V1))
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	240 break;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	241
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	242 S1 = calculate_step_vector2d (X1f, V1, tx, ty, cols, rows, h);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	243
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	244 // Runge Kutta - Heun's Scheme
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	245 const Vector2 S = {0.5 * (S0.x + S1.x),
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	246 0.5 * (S0.y + S1.y)};
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	247 Xnxt = add2d (X0f, S);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	248
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	249 X0f = vector_to_cell2d (Xnxt);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	250 if (! is_in_definition_set2d (X0f, cols, rows))
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	251 break;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	252
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	253 i++;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	254 buffer(i, 0) = Xnxt.x;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	255 buffer(i, 1) = Xnxt.y;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	256
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	257 if (i + 1 >= maxnverts)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	258 break;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	259 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	260
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	261 *nverts = i + 1;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	262 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	263
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	264 static inline double
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	265 trilinear (const double u111, const double u211, const double u121,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	266 const double u221, const double u112, const double u212,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	267 const double u122, const double u222,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	268 const double x, const double y, const double z)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	269 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	270 return (u111 * (1-x) * (1-y) * (1-z) +
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	271 u211 * x * (1-y) * (1-z) +
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	272 u121 * (1-x) * y * (1-z) +
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	273 u221 * x * y * (1-z) +
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	274 u112 * (1-x) * (1-y) * z +
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	275 u212 * x * (1-y) * z +
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	276 u122 * (1-x) * y * z +
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	277 u222 * x * y * z);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	278 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	279
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	280 static inline Cell3
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	281 vector_to_cell3d (const Vector3 X)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	282 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	283 Cell3 Z;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	284
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	285 number_to_fractional (&Z.idx, &Z.fcx, X.x);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	286 number_to_fractional (&Z.idy, &Z.fcy, X.y);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	287 number_to_fractional (&Z.idz, &Z.fcz, X.z);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	288
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	289 return (Z);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	290 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	291
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	292 static inline bool
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	293 is_in_definition_set3d (const Cell3 X, const octave_idx_type nx,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	294 const octave_idx_type ny, const octave_idx_type nz)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	295 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	296 return ( (((X.idx >= 0) && (X.idx < nx-1)) \|\|
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	297 ((X.idx == nx-1) && (X.fcx == 0.0))) &&
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	298 (((X.idy >= 0) && (X.idy < ny-1)) \|\|
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	299 ((X.idy == ny-1) && (X.fcy == 0.0))) &&
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	300 (((X.idz >= 0) && (X.idz < nz-1)) \|\|
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	301 ((X.idz == nz-1) && (X.fcz == 0.0))) );
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	302 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	303
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	304 static inline Vector3
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	305 add3d (const Cell3 X, const Vector3 Y)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	306 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	307 Vector3 Z = {X.idx + X.fcx + Y.x,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	308 X.idy + X.fcy + Y.y,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	309 X.idz + X.fcz + Y.z};
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	310
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	311 return (Z);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	312 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	313
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	314 static inline Vector3
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	315 vector_interpolation3d (const Cell3 X, const NDArray& u, const NDArray& v,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	316 const NDArray& w, const octave_idx_type nx,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	317 const octave_idx_type ny, const octave_idx_type nz)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	318 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	319 Vector3 V;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	320 double fcx, fcy, fcz;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	321 octave_idx_type idx, idy, idz;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	322
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	323 handle_border (&idx, &fcx, X.idx, X.fcx, nx);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	324 handle_border (&idy, &fcy, X.idy, X.fcy, ny);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	325 handle_border (&idz, &fcz, X.idz, X.fcz, nz);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	326
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	327 V.x = trilinear (u(idy, idx, idz), u(idy, idx+1, idz),
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	328 u(idy+1, idx, idz), u(idy+1, idx+1, idz),
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	329 u(idy, idx, idz+1), u(idy, idx+1, idz+1),
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	330 u(idy+1, idx, idz+1), u(idy+1, idx+1, idz+1),
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	331 fcx, fcy, fcz);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	332 V.y = trilinear (v(idy, idx, idz), v(idy, idx+1, idz),
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	333 v(idy+1, idx, idz), v(idy+1, idx+1, idz),
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	334 v(idy, idx, idz+1), v(idy, idx+1, idz+1),
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	335 v(idy+1, idx, idz+1), v(idy+1, idx+1, idz+1),
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	336 fcx, fcy, fcz);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	337 V.z = trilinear (w(idy, idx, idz), w(idy, idx+1, idz),
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	338 w(idy+1, idx, idz), w(idy+1, idx+1, idz),
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	339 w(idy, idx, idz+1), w(idy, idx+1, idz+1),
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	340 w(idy+1, idx, idz+1), w(idy+1, idx+1, idz+1),
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	341 fcx, fcy, fcz);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	342
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	343 return (V);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	344 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	345
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	346 static inline Vector3
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	347 calculate_step_vector3d (const Cell3 X, const Vector3 V,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	348 const RowVector& tx, const RowVector& ty, const RowVector& tz,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	349 const octave_idx_type nx, const octave_idx_type ny,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	350 const octave_idx_type nz, const double h)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	351 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	352 Vector3 S;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	353
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	354 const octave_idx_type idx = handle_border_index (X.idx, nx);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	355 const octave_idx_type idy = handle_border_index (X.idy, ny);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	356 const octave_idx_type idz = handle_border_index (X.idz, nz);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	357
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	358 const double x = V.x * tx(idx);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	359 const double y = V.y * ty(idy);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	360 const double z = V.z * tz(idz);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	361 const double n = 1.0 / sqrt (xx + yy + z*z);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	362 S.x = h * n * x;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	363 S.y = h * n * y;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	364 S.z = h * n * z;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	365
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	366 return (S);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	367 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	368
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	369 static inline bool
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	370 is_singular3d (const Vector3 V)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	371 {
29961 7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more John W. Eaton <jwe@octave.org> parents: 29958 diff changeset	372 return ((math::isnan (V.x) \|\| math::isnan (V.y) \|\| math::isnan (V.z)) \|\|
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	373 ((V.x == 0) && (V.y == 0) && (V.z == 0)));
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	374 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	375
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	376 static void
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	377 euler3d (const octave_idx_type nx, const octave_idx_type ny, const octave_idx_type nz,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	378 const NDArray& u, const NDArray& v, const NDArray& w,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	379 const RowVector& tx, const RowVector& ty, const RowVector& tz,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	380 const double zeta, const double xi, const double rho,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	381 const double h, const octave_idx_type maxnverts,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	382 Matrix& buffer, octave_idx_type *nverts)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	383 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	384 Vector3 V0, V1, S0, X1, Xnxt, S1;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	385 const Vector3 X0 = {zeta, xi, rho};
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	386 Cell3 X0f, X1f;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	387
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	388 octave_idx_type i = 0;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	389 buffer(i, 0) = X0.x;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	390 buffer(i, 1) = X0.y;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	391 buffer(i, 2) = X0.z;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	392
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	393 X0f = vector_to_cell3d (X0);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	394 while (true)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	395 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	396 if (! is_in_definition_set3d (X0f, nx, ny, nz))
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	397 break;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	398
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	399 V0 = vector_interpolation3d (X0f, u, v, w, nx, ny, nz);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	400 if (is_singular3d (V0))
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	401 break;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	402
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	403 S0 = calculate_step_vector3d (X0f, V0, tx, ty, tz, nx, ny, nz, h);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	404
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	405 X1 = add3d (X0f, S0);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	406
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	407 X1f = vector_to_cell3d (X1);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	408 if (! is_in_definition_set3d (X1f, nx, ny, nz))
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	409 break;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	410
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	411 V1 = vector_interpolation3d (X1f, u, v, w, nx, ny, nz);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	412 if (is_singular3d (V1))
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	413 break;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	414
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	415 S1 = calculate_step_vector3d (X1f, V1, tx, ty, tz, nx, ny, nz, h);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	416
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	417 // Runge Kutta - Heun's Scheme
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	418 const Vector3 S = {0.5 * (S0.x + S1.x),
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	419 0.5 * (S0.y + S1.y),
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	420 0.5 * (S0.z + S1.z)};
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	421 Xnxt = add3d (X0f, S);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	422
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	423 X0f = vector_to_cell3d (Xnxt);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	424 if (! is_in_definition_set3d (X0f, nx, ny, nz))
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	425 break;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	426
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	427 i++;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	428 buffer(i, 0) = Xnxt.x;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	429 buffer(i, 1) = Xnxt.y;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	430 buffer(i, 2) = Xnxt.z;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	431
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	432 if (i + 1 >= maxnverts)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	433 break;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	434
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	435 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	436
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	437 *nverts = i + 1;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	438 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	439
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	440 static octave_value
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	441 streameuler2d_internal (const octave_value_list& args)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	442 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	443
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	444 const int nargin = args.length ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	445 if (nargin != 8)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	446 print_usage ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	447
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	448 const Matrix U = args(0).matrix_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	449 const Matrix V = args(1).matrix_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	450 const RowVector TX = args(2).row_vector_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	451 const RowVector TY = args(3).row_vector_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	452 const double zeta = args(4).double_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	453 const double xi = args(5).double_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	454 const double h = args(6).double_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	455 const octave_idx_type maxnverts = args(7).idx_type_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	456
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	457 const octave_idx_type rows = U.rows ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	458 const octave_idx_type cols = U.columns ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	459
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	460 octave_idx_type nverts;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	461 Matrix buffer (maxnverts, 2);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	462
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	463 euler2d (cols, rows, U, V, TX, TY, zeta, xi, h, maxnverts,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	464 buffer, &nverts);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	465
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	466 Matrix xy = buffer.extract (0, 0, nverts-1, 1);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	467
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	468 return octave_value (xy);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	469 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	470
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	471 static octave_value
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	472 streameuler3d_internal (const octave_value_list& args, const char *fcn)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	473 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	474
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	475 const int nargin = args.length ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	476 if (nargin != 11)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	477 print_usage ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	478
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	479 const NDArray U = args(0).array_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	480 const NDArray V = args(1).array_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	481 const NDArray W = args(2).array_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	482 const RowVector TX = args(3).row_vector_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	483 const RowVector TY = args(4).row_vector_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	484 const RowVector TZ = args(5).row_vector_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	485 const double zeta = args(6).double_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	486 const double xi = args(7).double_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	487 const double rho = args(8).double_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	488 const double h = args(9).double_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	489 const octave_idx_type maxnverts = args(10).idx_type_value ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	490
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	491 const dim_vector dims = args(0).dims ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	492 const int ndims = dims.ndims ();
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	493 if (ndims != 3)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	494 error ("%s: dimension must be 3", fcn);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	495
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	496 octave_idx_type nverts;
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	497 Matrix buffer (maxnverts, 3);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	498
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	499 euler3d (dims(1), dims(0), dims(2), U, V, W, TX, TY, TZ, zeta, xi, rho,
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	500 h, maxnverts, buffer, &nverts);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	501
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	502 Matrix xyz = buffer.extract (0, 0, nverts-1, 2);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	503
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	504 return octave_value (xyz);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	505 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	506
27808 d7dfab7045d9 Make "streameuler2d" and "streameuler3d" internal functions. Markus Mützel <markus.muetzel@gmx.de> parents: 27795 diff changeset	507 DEFUN (__streameuler2d__, args, ,
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	508 doc: /* -- texinfo --
30888 32d2b6604a9f doc: Ensure documentation lists output argument when it exists for functions in libinterp/ Rik <rik@octave.org> parents: 30564 diff changeset	509 @deftypefn {} {@var{output} =} __streameuler2d__ (@var{U}, @var{V}, @var{TX}, @var{TY}, @var{ZETA}, @var{XI}, @var{H}, @var{MAXNVERTS})
31594 a3d8915b9f87 doc: more DOCSTRING changes ahead of 8.1 release. Rik <rik@octave.org> parents: 30888 diff changeset	510 Calculate the streamline in a vector field @code{[@var{U}, @var{V}]} starting
27810 f2b89a2e20b6 overhaul streamXXX.m series of functions. Rik <rik@octave.org> parents: 27808 diff changeset	511 from a seed point at position @code{[@var{ZETA}, @var{XI}]}. The integrator
f2b89a2e20b6 overhaul streamXXX.m series of functions. Rik <rik@octave.org> parents: 27808 diff changeset	512 used is Heun's Scheme. The step size can be controlled by @var{H}. The
f2b89a2e20b6 overhaul streamXXX.m series of functions. Rik <rik@octave.org> parents: 27808 diff changeset	513 Jacobian matrix can be defined for each grid cell by
f2b89a2e20b6 overhaul streamXXX.m series of functions. Rik <rik@octave.org> parents: 27808 diff changeset	514 @code{[@var{TX}, @var{TY}]}.
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	515
27808 d7dfab7045d9 Make "streameuler2d" and "streameuler3d" internal functions. Markus Mützel <markus.muetzel@gmx.de> parents: 27795 diff changeset	516 @seealso{streamline, stream2, stream3, __streameuler3d__}
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	517 @end deftypefn */)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	518 {
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	519 return streameuler2d_internal (args);
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	520 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	521
27808 d7dfab7045d9 Make "streameuler2d" and "streameuler3d" internal functions. Markus Mützel <markus.muetzel@gmx.de> parents: 27795 diff changeset	522 DEFUN (__streameuler3d__, args, ,
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	523 doc: /* -- texinfo --
30888 32d2b6604a9f doc: Ensure documentation lists output argument when it exists for functions in libinterp/ Rik <rik@octave.org> parents: 30564 diff changeset	524 @deftypefn {} {@var{output} =} __streameuler3d__ (@var{U}, @var{V}, @var{W}, @var{TX}, @var{TY}, @var{TZ}, @var{ZETA}, @var{XI}, @var{RHO}, @var{H}, @var{MAXNVERTS})
31594 a3d8915b9f87 doc: more DOCSTRING changes ahead of 8.1 release. Rik <rik@octave.org> parents: 30888 diff changeset	525 Calculate the streamline in a vector field @code{[@var{U}, @var{V}, @var{W}]}
27810 f2b89a2e20b6 overhaul streamXXX.m series of functions. Rik <rik@octave.org> parents: 27808 diff changeset	526 starting from a seed point at position
f2b89a2e20b6 overhaul streamXXX.m series of functions. Rik <rik@octave.org> parents: 27808 diff changeset	527 @code{[@var{ZETA}, @var{XI}, @var{RHO}]}. The integrator used is Heun's
f2b89a2e20b6 overhaul streamXXX.m series of functions. Rik <rik@octave.org> parents: 27808 diff changeset	528 Scheme. The step size can be controlled by @var{H}. The Jacobian matrix can
f2b89a2e20b6 overhaul streamXXX.m series of functions. Rik <rik@octave.org> parents: 27808 diff changeset	529 be defined for each grid cell by @code{[@var{TX}, @var{TY}, @var{TZ}]}.
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	530
27808 d7dfab7045d9 Make "streameuler2d" and "streameuler3d" internal functions. Markus Mützel <markus.muetzel@gmx.de> parents: 27795 diff changeset	531 @seealso{streamline, stream2, stream3, __streameuler2d__}
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	532 @end deftypefn */)
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	533 {
27808 d7dfab7045d9 Make "streameuler2d" and "streameuler3d" internal functions. Markus Mützel <markus.muetzel@gmx.de> parents: 27795 diff changeset	534 return streameuler3d_internal (args, "__streameuler3d__");
27759 992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	535 }
992e702ef0d7 Add stream* functions (patch #9859). Markus Meisinger <chloros2@gmx.de> parents: diff changeset	536
31605 e88a07dec498 maint: Use macros to begin/end C++ namespaces. Rik <rik@octave.org> parents: 31594 diff changeset	537 OCTAVE_END_NAMESPACE(octave)

Mercurial > octave

annotate libinterp/corefcn/stream-euler.cc @ 31605:e88a07dec498 stable