Linear Algebra and the C Language/a0ae
Code study: dot_R();
/* ------------------------------------ */
/* ------- u.v = v^t u ---------------- */
/* ------------------------------------ */
double dot_R(
double **u,
double **v
)
{
double **vt = i_mR(R1, rsize_R(v));
double **vtu = i_mR(R1, C1);
double udotv;
transpose_mR(v,vt);
mul_mR(vt,u,vtu);
udotv = vtu[R1][C1];
f_mR(vt);
f_mR(vtu);
return (udotv);
}
/* ------------------------------------ */
/* ------------------------------------ */