data
*dmn
, *dn1
, *dn2
, *dn3
, *dm
;
if(p
->type
!= DA
|| q
->type
!= DA
)
if((p
->rank
!= 1 && p
->rank
!= 2) || q
->rank
!= 2)
if(m
< n
|| m
!= p
->dim
[0])
a
= alloc(n
*(SINT
+ SDAT
*m
+ SDAT
*3) + SDAT
*m
);
a
= lsq(dmn
, dn1
, dn2
, dn3
, dm
, in
, m
, n
, o
, p
->datap
, q
->datap
);
lsq(dmn
, dn1
, dn2
, dn3
, dm
, in
, m
, n
, p
, d1
, d2
)
data
*dmn
, *dn1
, *dn2
, *dn3
, *dm
;
register data
*dp1
, *dp2
;
solve(m
, n
, dmn
, dn2
, in
, dm
, dn3
);
f1
=+ d2
[i
*n
+j
] * *dp2
++;
solve(m
, n
, dmn
, dn2
, in
, dm
, dn1
);
f1
=+ d2
[i
*n
+j
] * *dp2
++;
solve(m
, n
, dmn
, dn2
, in
, dm
, dn1
);
solve(m
, n
, dmn
, dn2
, in
, dm
, dn1
)
data
*dmn
, *dn1
, *dn2
, *dm
;
register data
*dp1
, *dp2
;
dn1
[in
[n
-1]] = *--dp1
/ *--dp2
;
f1
=+ dmn
[k
] * dn1
[in
[j
]];
dn1
[in
[i
]] = -f1
/ *--dp2
;