19 #define REPORT { static ExeCounter ExeCount(__LINE__,12); ++ExeCount; } 61 { sum += (
long_Real)*s * *s;
if (!(--i))
break; s +=
d; }
69 Real* s = store;
int d = spacing;
78 { sum += (
long_Real)*s * *s1;
if (!(--i))
break; s +=
d; s1 += d1; }
85 int i = n;
Real* s = store;
int d = spacing;
94 { *s += *s1 * r;
if (!(--i))
break; s +=
d; s1 += d1; }
100 int i = n;
Real* s = store;
int d = spacing;
108 if (i)
for (;;) { *s = *s1 / r;
if (!(--i))
break; s +=
d; s1 += d1; }
114 int i = n;
Real* s = store;
int d = spacing;
116 if (i)
for (;;) { *s /= r;
if (!(--i))
break; s +=
d; }
122 int i = n;
Real* s = store;
int d = spacing;
124 if (i)
for (;;) { *s = - *s;
if (!(--i))
break; s +=
d; }
130 int i = n;
Real* s = store;
int d = spacing;
132 if (i)
for (;;) { *s = 0.0;
if (!(--i))
break; s +=
d; }
153 Real z = *u * x - *v * y; *v = *u * y + *v * x; *u = z;
180 *u -= s * (zv + zu * tau); *v += s * (zu - zv * tau);
194 if (f==0 && g==0) { c=1.0; s=0.0;
return 0.0; }
195 Real af = f>=0 ? f : -f;
196 Real ag = g>=0 ? g : -g;
200 Real h = g/f;
Real sq = sqrt(1.0+h*h);
202 c = 1.0/sq; s = h/sq;
return sq*f;
207 Real h = f/g;
Real sq = sqrt(1.0+h*h);
209 s = 1.0/sq; c = h/sq;
return sq*g;
Real operator*(const RectMatrixRowCol &) const
void Reset(Real *st, int nx, int sp, int sh)
Internal newmat exception - shouldn't happen.
void Rotate(RectMatrixCol &U, RectMatrixCol &V, Real tau, Real s)
Real pythag(Real f, Real g, Real &c, Real &s)
void Reset(const Matrix &, int, int, int)
void ComplexScale(RectMatrixCol &U, RectMatrixCol &V, Real x, Real y)
void Divide(const RectMatrixRowCol &, Real)
The usual rectangular matrix.
Real sum(const BaseMatrix &B)
void AddScaled(const RectMatrixRowCol &, Real)
void Reset(const Matrix &, int, int, int)