#include "f2c.h"
void dlinrt(int n, double *a, int lda, int ipath, double *ainv, int ldainv)
#define A(I,J) a[(I) + (J)* ( lda)]
#define AINV(I,J) ainv[(I) + (J)* ( ldainv)]
{
  int i,j;
  char uplo, diag='n';
  int info;
  
  for (i=0; i<=n-1; i++) 
     for (j=0; j<=i-1; j++) 
        AINV(i,j)=A(i,j); 
  if (ipath==1) uplo='l'; else uplo='u';
  dtrtri_(&uplo, &diag, &n, ainv, &ldainv, &info);
}