Mercurial > mxe-octave

comparison src/of-lssa-1-fixes.patch @ 4198:2f1d93109f43

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
temporary fixes for of-lssa package
author John W. Eaton <jwe@octave.org>
date Wed, 31 Aug 2016 15:01:40 -0400
parents
children
comparison
equal deleted inserted replaced
4197:667aa26d9dd4 4198:2f1d93109f43
1 # HG changeset patch
2 # User John Donoghue
3 # Date 1467900324 14400
4 # Thu Jul 07 10:05:24 2016 -0400
5 # Node ID 59af0bcd5cafce29162db3dfc42dcf326f56d8be
6 # Parent 3727b85ea021033ee96fb0b017eefb73fb2f7c35
7 Use complex real/imag set method to assign values (Bug #48424)
8
9 * src/fastlscomplex.cc
10 (flscomplex): Use record_current, record_next real/imag set method to assign +/- values
11
12 diff --git a/src/fastlscomplex.cc b/src/fastlscomplex.cc
13 --- a/src/fastlscomplex.cc
14 +++ b/src/fastlscomplex.cc
15 @@ -372,8 +372,10 @@ bool flscomplex (const RowVector & tvec,
16
17 --tpra;
18 h = *tpra * *exp_ptr;
19 - record_current->power_series[p].real() -= h.imag();
20 - record_current->power_series[p].imag() += h.real();
21 + record_current->power_series[p].real (
22 + record_current->power_series[p].real() - h.imag());
23 + record_current->power_series[p].imag (
24 + record_current->power_series[p].imag() + h.real());
25
26 if (++exp_ptr >= exp_pse_ptr )
27 break;
28 @@ -386,8 +388,10 @@ bool flscomplex (const RowVector & tvec,
29
30 --tpra;
31 h = -*tpra * *exp_ptr;
32 - record_current->power_series[p].real() -= h.imag();
33 - record_current->power_series[p].imag() += h.real();
34 + record_current->power_series[p].real (
35 + record_current->power_series[p].real() - h.imag());
36 + record_current->power_series[p].imag (
37 + record_current->power_series[p].imag() + h.real());
38
39 if (++exp_ptr >= exp_pse_ptr)
40 break;
41 @@ -422,8 +426,10 @@ bool flscomplex (const RowVector & tvec,
42
43 tprb -= 2;
44 h = *tprb * *exp_ptr;
45 - record_current->power_series[p].real() -= h.imag();
46 - record_current->power_series[p].imag() += h.real();
47 + record_current->power_series[p].real (
48 + record_current->power_series[p].real() - h.imag());
49 + record_current->power_series[p].imag (
50 + record_current->power_series[p].imag() + h.real());
51
52 if ( ++exp_ptr >= exp_pse_ptr )
53 break;
54 @@ -436,8 +442,10 @@ bool flscomplex (const RowVector & tvec,
55
56 tprb -= 2;
57 h = - *tprb * *exp_ptr;
58 - record_current->power_series[p].real() -= h.imag();
59 - record_current->power_series[p].imag() += h.real();
60 + record_current->power_series[p].real (
61 + record_current->power_series[p].real() - h.imag());
62 + record_current->power_series[p].imag (
63 + record_current->power_series[p].imag() + h.real());
64
65 if (++exp_ptr >= exp_pse_ptr)
66 break;
67 @@ -467,8 +475,10 @@ bool flscomplex (const RowVector & tvec,
68
69 --tpra;
70 h = *tpra * *exp_ptr;
71 - record_next->power_series[q].real() -= h.imag();
72 - record_next->power_series[q].imag() += h.real();
73 + record_next->power_series[q].real (
74 + record_next->power_series[q].real() - h.imag());
75 + record_next->power_series[q].imag (
76 + record_next->power_series[q].imag() + h.real());
77
78 if (++exp_ptr >= exp_pse_ptr)
79 break;
80 @@ -481,8 +491,10 @@ bool flscomplex (const RowVector & tvec,
81
82 --tpra;
83 h = -*tpra * *exp_ptr;
84 - record_next->power_series[q].real() -= h.imag();
85 - record_next->power_series[q].imag() += h.real();
86 + record_next->power_series[q].real (
87 + record_next->power_series[q].real() - h.imag());
88 + record_next->power_series[q].imag (
89 + record_next->power_series[q].imag() + h.real());
90
91 if (++exp_ptr >= exp_pse_ptr)
92 break;