#include <assert.h>
#include "TCernLib.h"
#include "TMath.h"
#include "TArrayD.h"
#include "TError.h"
Go to the source code of this file.
|
#define | TCL_MXMAD(n_, a, b, c, i, j, k) |
|
#define | TCL_MXMLRT(n__, a, b, c,ni, nj) |
|
#define | TCL_MXTRP(a, b, i, j) |
|
#define | TCL_TRAAT(a, s, m, n) |
|
#define | TCL_TRAL(a, u, b, m,n) |
|
#define | TCL_TRALT(a, u, b, m, n) |
|
#define | TCL_TRAS(a, s, b, m, n) |
|
#define | TCL_TRASAT(a, s, r__, m, n) |
|
|
| ClassImp (TCL) float *TCL |
|
*traat F translated by | f2c (version 19970219).**BEGIN_HTML<!--*/float *TCL |
| Symmetric Multiplication of Rectangular Matrices CERN PROGLIB# F112 TRAAT .VERSION KERNFOR 4.15 861204 ORIG. More...
|
|
#define TCL_MXMAD |
( |
|
n_, |
|
|
|
a, |
|
|
|
b, |
|
|
|
c, |
|
|
|
i, |
|
|
|
j, |
|
|
|
k |
|
) |
| |
#define TCL_MXMLRT |
( |
|
n__, |
|
|
|
a, |
|
|
|
b, |
|
|
|
c, |
|
|
|
ni, |
|
|
|
nj |
|
) |
| |
Value:if (ni <= 0 || nj <= 0)
return 0;
\
int ia, ib, ic, ja, kc,
ii, jj, kj, ki, ia1, ib1, ic1, ja1; \
int ipa = 1;
int jpa = nj;
\
if (n__ == 1) { ipa = ni; jpa = 1; } \
\
\
for (ii = 1; ii <= ni; ++
ii, ic1+=ni, ia1+=jpa) { \
for (kc = 1; kc <= ni; ++kc,ic++)
c[ic] = 0.; \
for (jj = 1; jj <= nj; ++jj,++ib1,ja1 += ipa) { \
for (kj = 1;kj <= nj;++kj,ia+=ipa,ib += nj) \
for (ki = 1; ki <= ni; ++ki,++ic,ja += jpa) \
} \
}
if(pyself &&pyself!=Py_None)
Definition at line 130 of file TCernLib.cxx.
Referenced by TCL::mxmlrt_0_().
#define TCL_MXTRP |
( |
|
a, |
|
|
|
b, |
|
|
|
i, |
|
|
|
j |
|
) |
| |
Value:if (i == 0 || j == 0) return 0; \
for (
int k = 1; k <= j; ++k) \
for (
int l = 1;
l <= i; ++
l,ia += j,++ib) b[ib] =
a[ia]; }
Definition at line 209 of file TCernLib.cxx.
Referenced by TCL::mxtrp().
#define TCL_TRAAT |
( |
|
a, |
|
|
|
s, |
|
|
|
m, |
|
|
|
n |
|
) |
| |
Value: \
int ipiv, i, j, ipivn, ia, is, iat;
\
for (i = 1; i <=
m; ++i) { \
ipiv = ia; \
for (j = 1; j <= i; ++j) { \
ia = ipiv; \
sum = 0.; \
do { \
++ia; ++iat; \
} while (ia < ipivn); \
++is; \
s[is] = sum; \
} \
} \
s++;
Definition at line 261 of file TCernLib.cxx.
Referenced by f2c().
#define TCL_TRAL |
( |
|
a, |
|
|
|
u, |
|
|
|
b, |
|
|
|
m, |
|
|
|
n |
|
) |
| |
Value:int indu, i, j, k, ia, ib, iu;
\
for (i = 1; i <=
m; ++i) { \
for (j = 1; j <=
n; ++j) { \
indu += j; \
ia = ib; \
iu = indu; \
for (k = j; k <=
n; ++k) {\
++ia; \
iu += k; \
} \
b[ib] = sum; \
++ib; \
} \
} \
b++;
Definition at line 324 of file TCernLib.cxx.
Referenced by TCL::tral().
#define TCL_TRALT |
( |
|
a, |
|
|
|
u, |
|
|
|
b, |
|
|
|
m, |
|
|
|
n |
|
) |
| |
Value:int indu, j, k, ia, ib, iu;
\
indu = (
n *
n +
n) / 2; \
do { \
for (j = 1; j <=
n; ++j) { \
ia = ib; \
for (k = j; k <=
n; ++k) {\
--ia; --iu; \
} \
b[ib] = sum; \
--ib; \
} \
} while (ib > 0); \
++b;
#define TCL_TRAS |
( |
|
a, |
|
|
|
s, |
|
|
|
b, |
|
|
|
m, |
|
|
|
n |
|
) |
| |
Value:int inds, i__, j, k, ia, ib, is;
\
ib = 0; inds = 0; i__ = 0; \
do { \
inds += i__; \
ia = 0; \
for (j = 1; j <=
m; ++j) { \
is = inds; \
sum = 0.; \
k = 0; \
++ia; \
++k; \
b[ib] = sum; \
ib += n; \
} \
++i__; \
++b;
if(pyself &&pyself!=Py_None)
#define TCL_TRASAT |
( |
|
a, |
|
|
|
s, |
|
|
|
r__, |
|
|
|
m, |
|
|
|
n |
|
) |
| |
*traat F translated by f2c |
( |
version |
19970219 | ) |
|
Symmetric Multiplication of Rectangular Matrices CERN PROGLIB# F112 TRAAT .VERSION KERNFOR 4.15 861204 ORIG.
18/12/74 WH
Definition at line 290 of file TCernLib.cxx.