/* tsimeq.h  solve A[n,n]*X[n]=Y[n] for vector X  threaded version */
/* np is number of worker processors, np+1 includes master */
/* A preserved */
void tsimeq(int np, int n, double A[], double Y[], double X[]);
/* last column of B is Y, B overwritten */
void tsimeqb(int np, int n, double B[], double X[]);