Mercurial > mxe-octave
changeset 4198:2f1d93109f43
temporary fixes for of-lssa package
| author | John W. Eaton <jwe@octave.org> |
|---|---|
| date | Wed, 31 Aug 2016 15:01:40 -0400 |
| parents | 667aa26d9dd4 |
| children | f7fce04ebafa |
| files | src/of-lssa-1-fixes.patch |
| diffstat | 1 files changed, 92 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/src/of-lssa-1-fixes.patch Wed Aug 31 15:01:40 2016 -0400 @@ -0,0 +1,92 @@ +# HG changeset patch +# User John Donoghue +# Date 1467900324 14400 +# Thu Jul 07 10:05:24 2016 -0400 +# Node ID 59af0bcd5cafce29162db3dfc42dcf326f56d8be +# Parent 3727b85ea021033ee96fb0b017eefb73fb2f7c35 +Use complex real/imag set method to assign values (Bug #48424) + +* src/fastlscomplex.cc + (flscomplex): Use record_current, record_next real/imag set method to assign +/- values + +diff --git a/src/fastlscomplex.cc b/src/fastlscomplex.cc +--- a/src/fastlscomplex.cc ++++ b/src/fastlscomplex.cc +@@ -372,8 +372,10 @@ bool flscomplex (const RowVector & tvec, + + --tpra; + h = *tpra * *exp_ptr; +- record_current->power_series[p].real() -= h.imag(); +- record_current->power_series[p].imag() += h.real(); ++ record_current->power_series[p].real ( ++ record_current->power_series[p].real() - h.imag()); ++ record_current->power_series[p].imag ( ++ record_current->power_series[p].imag() + h.real()); + + if (++exp_ptr >= exp_pse_ptr ) + break; +@@ -386,8 +388,10 @@ bool flscomplex (const RowVector & tvec, + + --tpra; + h = -*tpra * *exp_ptr; +- record_current->power_series[p].real() -= h.imag(); +- record_current->power_series[p].imag() += h.real(); ++ record_current->power_series[p].real ( ++ record_current->power_series[p].real() - h.imag()); ++ record_current->power_series[p].imag ( ++ record_current->power_series[p].imag() + h.real()); + + if (++exp_ptr >= exp_pse_ptr) + break; +@@ -422,8 +426,10 @@ bool flscomplex (const RowVector & tvec, + + tprb -= 2; + h = *tprb * *exp_ptr; +- record_current->power_series[p].real() -= h.imag(); +- record_current->power_series[p].imag() += h.real(); ++ record_current->power_series[p].real ( ++ record_current->power_series[p].real() - h.imag()); ++ record_current->power_series[p].imag ( ++ record_current->power_series[p].imag() + h.real()); + + if ( ++exp_ptr >= exp_pse_ptr ) + break; +@@ -436,8 +442,10 @@ bool flscomplex (const RowVector & tvec, + + tprb -= 2; + h = - *tprb * *exp_ptr; +- record_current->power_series[p].real() -= h.imag(); +- record_current->power_series[p].imag() += h.real(); ++ record_current->power_series[p].real ( ++ record_current->power_series[p].real() - h.imag()); ++ record_current->power_series[p].imag ( ++ record_current->power_series[p].imag() + h.real()); + + if (++exp_ptr >= exp_pse_ptr) + break; +@@ -467,8 +475,10 @@ bool flscomplex (const RowVector & tvec, + + --tpra; + h = *tpra * *exp_ptr; +- record_next->power_series[q].real() -= h.imag(); +- record_next->power_series[q].imag() += h.real(); ++ record_next->power_series[q].real ( ++ record_next->power_series[q].real() - h.imag()); ++ record_next->power_series[q].imag ( ++ record_next->power_series[q].imag() + h.real()); + + if (++exp_ptr >= exp_pse_ptr) + break; +@@ -481,8 +491,10 @@ bool flscomplex (const RowVector & tvec, + + --tpra; + h = -*tpra * *exp_ptr; +- record_next->power_series[q].real() -= h.imag(); +- record_next->power_series[q].imag() += h.real(); ++ record_next->power_series[q].real ( ++ record_next->power_series[q].real() - h.imag()); ++ record_next->power_series[q].imag ( ++ record_next->power_series[q].imag() + h.real()); + + if (++exp_ptr >= exp_pse_ptr) + break;