Mercurial > octave

comparison libinterp/corefcn/ellipj.cc @ 32660:f53ac65ffba6

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
maint: New method rwdata() as clearer alternative to fortran_vec(). * etc/NEWS.10.md: Announce new method and recommend its use in future code. * doc/interpreter/external.txi: Update Octave manual to use rwdata. * liboctave/array/Array.h (rwdata): New function prototype. * liboctave/array/Array.h (fortran_vec): Make inline function that calls rwdata(). * liboctave/array/DiagArray2.h (rwdata): New function prototype. * liboctave/array/DiagArray2.h (fortran_vec): Make inline function that calls rwdata(). * examples/code/fortrandemo.cc, libgui/graphics/QtHandlesUtils.cc, libinterp/corefcn/Cell.cc, libinterp/corefcn/__dsearchn__.cc, libinterp/corefcn/__ichol__.cc, libinterp/corefcn/__ilu__.cc, libinterp/corefcn/__lin_interpn__.cc, libinterp/corefcn/__magick_read__.cc, libinterp/corefcn/__pchip_deriv__.cc, libinterp/corefcn/amd.cc, libinterp/corefcn/data.cc, libinterp/corefcn/dot.cc, libinterp/corefcn/ellipj.cc, libinterp/corefcn/filter.cc, libinterp/corefcn/gcd.cc, libinterp/corefcn/gl-render.cc, libinterp/corefcn/graphics.cc, libinterp/corefcn/graphics.in.h, libinterp/corefcn/kron.cc, libinterp/corefcn/ls-mat4.cc, libinterp/corefcn/ls-mat5.cc, libinterp/corefcn/mex.cc, libinterp/corefcn/oct-map.cc, libinterp/corefcn/oct-stream.cc, libinterp/corefcn/ordqz.cc, libinterp/corefcn/ordschur.cc, libinterp/corefcn/perms.cc, libinterp/corefcn/psi.cc, libinterp/corefcn/quadcc.cc, libinterp/corefcn/qz.cc, libinterp/corefcn/rand.cc, libinterp/corefcn/sqrtm.cc, libinterp/corefcn/strfind.cc, libinterp/corefcn/symrcm.cc, libinterp/corefcn/tril.cc, libinterp/corefcn/typecast.cc, libinterp/corefcn/xdiv.cc, libinterp/dldfcn/__delaunayn__.cc, libinterp/dldfcn/__glpk__.cc, libinterp/dldfcn/__ode15__.cc, libinterp/dldfcn/__voronoi__.cc, libinterp/dldfcn/audioread.cc, libinterp/dldfcn/convhulln.cc, libinterp/octave-value/ov-base-int.cc, libinterp/octave-value/ov-bool-mat.cc, libinterp/octave-value/ov-ch-mat.cc, libinterp/octave-value/ov-cx-diag.cc, libinterp/octave-value/ov-cx-mat.cc, libinterp/octave-value/ov-flt-cx-diag.cc, libinterp/octave-value/ov-flt-cx-mat.cc, libinterp/octave-value/ov-flt-re-diag.cc, libinterp/octave-value/ov-flt-re-mat.cc, libinterp/octave-value/ov-intx.h, libinterp/octave-value/ov-java.cc, libinterp/octave-value/ov-perm.cc, libinterp/octave-value/ov-re-diag.cc, libinterp/octave-value/ov-re-mat.cc, libinterp/octave-value/ov-str-mat.cc, liboctave/array/Array-base.cc, liboctave/array/Array-util.cc, liboctave/array/CColVector.cc, liboctave/array/CDiagMatrix.cc, liboctave/array/CMatrix.cc, liboctave/array/CNDArray.cc, liboctave/array/CRowVector.cc, liboctave/array/CSparse.cc, liboctave/array/DiagArray2.h, liboctave/array/MArray.cc, liboctave/array/PermMatrix.cc, liboctave/array/Range.h, liboctave/array/Sparse.cc, liboctave/array/boolSparse.cc, liboctave/array/dColVector.cc, liboctave/array/dDiagMatrix.cc, liboctave/array/dMatrix.cc, liboctave/array/dNDArray.cc, liboctave/array/dRowVector.cc, liboctave/array/dSparse.cc, liboctave/array/fCColVector.cc, liboctave/array/fCDiagMatrix.cc, liboctave/array/fCMatrix.cc, liboctave/array/fCNDArray.cc, liboctave/array/fCRowVector.cc, liboctave/array/fColVector.cc, liboctave/array/fMatrix.cc, liboctave/array/fNDArray.cc, liboctave/array/fRowVector.cc, liboctave/array/idx-vector.cc, liboctave/numeric/CollocWt.cc, liboctave/numeric/DASPK.cc, liboctave/numeric/DASRT.cc, liboctave/numeric/DASSL.cc, liboctave/numeric/EIG.cc, liboctave/numeric/LSODE.cc, liboctave/numeric/Quad.cc, liboctave/numeric/aepbalance.cc, liboctave/numeric/bsxfun-defs.cc, liboctave/numeric/chol.cc, liboctave/numeric/eigs-base.cc, liboctave/numeric/fEIG.cc, liboctave/numeric/gepbalance.cc, liboctave/numeric/gsvd.cc, liboctave/numeric/hess.cc, liboctave/numeric/lu.cc, liboctave/numeric/oct-convn.cc, liboctave/numeric/oct-rand.cc, liboctave/numeric/qr.cc, liboctave/numeric/qrp.cc, liboctave/numeric/schur.cc, liboctave/numeric/sparse-dmsolve.cc, liboctave/numeric/sparse-lu.cc, liboctave/numeric/sparse-qr.cc, liboctave/numeric/svd.cc, liboctave/operators/mx-inlines.cc, liboctave/operators/mx-op-defs.h, liboctave/util/oct-base64.cc, liboctave/util/oct-binmap.h, liboctave/util/str-vec.cc: Replace calls to fortran_vec() with calls to rwdata().
author Rik <rik@octave.org>
date Wed, 27 Dec 2023 16:55:14 -0800
parents 4b601ca024d5
children
comparison
equal deleted inserted replaced
32659:71a23564ceca 32660:f53ac65ffba6
121 121
122 ComplexNDArray sn (sz_u), cn (sz_u), dn (sz_u); 122 ComplexNDArray sn (sz_u), cn (sz_u), dn (sz_u);
123 NDArray err (sz_u); 123 NDArray err (sz_u);
124 124
125 const Complex *pu = u.data (); 125 const Complex *pu = u.data ();
126 Complex *psn = sn.fortran_vec (); 126 Complex *psn = sn.rwdata ();
127 Complex *pcn = cn.fortran_vec (); 127 Complex *pcn = cn.rwdata ();
128 Complex *pdn = dn.fortran_vec (); 128 Complex *pdn = dn.rwdata ();
129 double *perr = err.fortran_vec (); 129 double *perr = err.rwdata ();
130 octave_idx_type nel = u.numel (); 130 octave_idx_type nel = u.numel ();
131 131
132 for (octave_idx_type i = 0; i < nel; i++) 132 for (octave_idx_type i = 0; i < nel; i++)
133 math::ellipj (pu[i], m, psn[i], pcn[i], pdn[i], perr[i]); 133 math::ellipj (pu[i], m, psn[i], pcn[i], pdn[i], perr[i]);
134 134
151 151
152 NDArray sn (sz_m), cn (sz_m), dn (sz_m); 152 NDArray sn (sz_m), cn (sz_m), dn (sz_m);
153 NDArray err (sz_m); 153 NDArray err (sz_m);
154 154
155 const double *pm = m.data (); 155 const double *pm = m.data ();
156 double *psn = sn.fortran_vec (); 156 double *psn = sn.rwdata ();
157 double *pcn = cn.fortran_vec (); 157 double *pcn = cn.rwdata ();
158 double *pdn = dn.fortran_vec (); 158 double *pdn = dn.rwdata ();
159 double *perr = err.fortran_vec (); 159 double *perr = err.rwdata ();
160 octave_idx_type nel = m.numel (); 160 octave_idx_type nel = m.numel ();
161 161
162 for (octave_idx_type i = 0; i < nel; i++) 162 for (octave_idx_type i = 0; i < nel; i++)
163 math::ellipj (u, pm[i], psn[i], pcn[i], pdn[i], perr[i]); 163 math::ellipj (u, pm[i], psn[i], pcn[i], pdn[i], perr[i]);
164 164
171 171
172 ComplexNDArray sn (sz_m), cn (sz_m), dn (sz_m); 172 ComplexNDArray sn (sz_m), cn (sz_m), dn (sz_m);
173 NDArray err (sz_m); 173 NDArray err (sz_m);
174 174
175 const double *pm = m.data (); 175 const double *pm = m.data ();
176 Complex *psn = sn.fortran_vec (); 176 Complex *psn = sn.rwdata ();
177 Complex *pcn = cn.fortran_vec (); 177 Complex *pcn = cn.rwdata ();
178 Complex *pdn = dn.fortran_vec (); 178 Complex *pdn = dn.rwdata ();
179 double *perr = err.fortran_vec (); 179 double *perr = err.rwdata ();
180 octave_idx_type nel = m.numel (); 180 octave_idx_type nel = m.numel ();
181 181
182 for (octave_idx_type i = 0; i < nel; i++) 182 for (octave_idx_type i = 0; i < nel; i++)
183 math::ellipj (u, pm[i], psn[i], pcn[i], pdn[i], perr[i]); 183 math::ellipj (u, pm[i], psn[i], pcn[i], pdn[i], perr[i]);
184 184
220 NDArray sn (sz_m), cn (sz_m), dn (sz_m); 220 NDArray sn (sz_m), cn (sz_m), dn (sz_m);
221 NDArray err (sz_m); 221 NDArray err (sz_m);
222 222
223 const double *pu = u.data (); 223 const double *pu = u.data ();
224 const double *pm = m.data (); 224 const double *pm = m.data ();
225 double *psn = sn.fortran_vec (); 225 double *psn = sn.rwdata ();
226 double *pcn = cn.fortran_vec (); 226 double *pcn = cn.rwdata ();
227 double *pdn = dn.fortran_vec (); 227 double *pdn = dn.rwdata ();
228 double *perr = err.fortran_vec (); 228 double *perr = err.rwdata ();
229 octave_idx_type nel = m.numel (); 229 octave_idx_type nel = m.numel ();
230 230
231 for (octave_idx_type i = 0; i < nel; i++) 231 for (octave_idx_type i = 0; i < nel; i++)
232 math::ellipj (pu[i], pm[i], psn[i], pcn[i], pdn[i], perr[i]); 232 math::ellipj (pu[i], pm[i], psn[i], pcn[i], pdn[i], perr[i]);
233 233
268 ComplexNDArray sn (sz_m), cn (sz_m), dn (sz_m); 268 ComplexNDArray sn (sz_m), cn (sz_m), dn (sz_m);
269 NDArray err (sz_m); 269 NDArray err (sz_m);
270 270
271 const Complex *pu = u.data (); 271 const Complex *pu = u.data ();
272 const double *pm = m.data (); 272 const double *pm = m.data ();
273 Complex *psn = sn.fortran_vec (); 273 Complex *psn = sn.rwdata ();
274 Complex *pcn = cn.fortran_vec (); 274 Complex *pcn = cn.rwdata ();
275 Complex *pdn = dn.fortran_vec (); 275 Complex *pdn = dn.rwdata ();
276 double *perr = err.fortran_vec (); 276 double *perr = err.rwdata ();
277 octave_idx_type nel = m.numel (); 277 octave_idx_type nel = m.numel ();
278 278
279 for (octave_idx_type i = 0; i < nel; i++) 279 for (octave_idx_type i = 0; i < nel; i++)
280 math::ellipj (pu[i], pm[i], psn[i], pcn[i], pdn[i], perr[i]); 280 math::ellipj (pu[i], pm[i], psn[i], pcn[i], pdn[i], perr[i]);
281 281