9 void otrsm(
double *sm,
double *a,
double *
b,
int n)
11 double z, *q0, *p, *s, *
t;
15 q0 = (
double *)calloc(n,
sizeof(
double));
16 for (i = 0; i < n; ++i) {
17 for (j = 0,
t =
b; j < n; ++j) {
18 for (k = 0, s = a + i * n, z = 0.; k < n; ++k)
22 for (j = 0, p = sm + i,
t = a; j <= i; ++j, p += n) {
23 for (k = 0, s = q0, z = 0.; k < n; ++k)
void otrsm(double *sm, double *a, double *b, int n)