Linear Algebra and the C Language/a017

/* ------------------------------------ */
/*  Save as :   vbtm.h                  */
/* ------------------------------------ */
double **transpose_mR(
double **A,
double **A_t
)
{
int r;
int c;

  canItranspose_mR(A,A_t,"transpose_mR();","(A or A_t)");

  for   ( r=R1; r<A[R_SIZE][C0]; r++)
    for ( c=C1; c<A[C_SIZE][C0]; c++)

        A_t[c][r]  = A[r][c];
        
return(A_t);
}
/* ------------------------------------ */
/* ------------------------------------ */
double trace_R(
double **A
)
{
int    r;
double trace = 0.0;

  isquare_mR(A,"trace_mR();","(A)"); 

  for( r=R1; r<A[R_SIZE][C0]; r++)

          trace += A[r][r];

  return(trace);
}
/* ------------------------------------ */
/* ------------------------------------ */