/*--------------------------------------------------------------------
  NAS Parallel Benchmarks 2.3 UPC version - BT source code tree
  2004 - GWU - HPCL -- check bt.c for complete information
  --------------------------------------------------------------------*/

#include <stdio.h>
#include <upc_relaxed.h>
#include "functions.h"

void binvrhs( double *lh, 
	      double *r )
{
  double pivot, coeff;
  
  /*--------------------------------------------------------------------
    c     
    c-------------------------------------------------------------------*/
  
  pivot = 1.00/lh[0];
  lh[5+0] = lh[5+0]*pivot;
  lh[10+0] = lh[10+0]*pivot;
  lh[15+0] = lh[15+0]*pivot;
  lh[20+0] = lh[20+0]*pivot;
  r[0]   = r[0]  *pivot;

  coeff = lh[1];
  lh[5+1]= lh[5+1] - coeff*lh[5+0];
  lh[10+1]= lh[10+1] - coeff*lh[10+0];
  lh[15+1]= lh[15+1] - coeff*lh[15+0];
  lh[20+1]= lh[20+1] - coeff*lh[20+0];
  r[1]   = r[1]   - coeff*r[0];

  coeff = lh[2];
  lh[5+2]= lh[5+2] - coeff*lh[5+0];
  lh[10+2]= lh[10+2] - coeff*lh[10+0];
  lh[15+2]= lh[15+2] - coeff*lh[15+0];
  lh[20+2]= lh[20+2] - coeff*lh[20+0];
  r[2]   = r[2]   - coeff*r[0];

  coeff = lh[3];
  lh[5+3]= lh[5+3] - coeff*lh[5+0];
  lh[10+3]= lh[10+3] - coeff*lh[10+0];
  lh[15+3]= lh[15+3] - coeff*lh[15+0];
  lh[20+3]= lh[20+3] - coeff*lh[20+0];
  r[3]   = r[3]   - coeff*r[0];

  coeff = lh[4];
  lh[5+4]= lh[5+4] - coeff*lh[5+0];
  lh[10+4]= lh[10+4] - coeff*lh[10+0];
  lh[15+4]= lh[15+4] - coeff*lh[15+0];
  lh[20+4]= lh[20+4] - coeff*lh[20+0];
  r[4]   = r[4]   - coeff*r[0];


  pivot = 1.00/lh[5+1];
  lh[10+1] = lh[10+1]*pivot;
  lh[15+1] = lh[15+1]*pivot;
  lh[20+1] = lh[20+1]*pivot;
  r[1]   = r[1]  *pivot;

  coeff = lh[5+0];
  lh[10+0]= lh[10+0] - coeff*lh[10+1];
  lh[15+0]= lh[15+0] - coeff*lh[15+1];
  lh[20+0]= lh[20+0] - coeff*lh[20+1];
  r[0]   = r[0]   - coeff*r[1];

  coeff = lh[5+2];
  lh[10+2]= lh[10+2] - coeff*lh[10+1];
  lh[15+2]= lh[15+2] - coeff*lh[15+1];
  lh[20+2]= lh[20+2] - coeff*lh[20+1];
  r[2]   = r[2]   - coeff*r[1];

  coeff = lh[5+3];
  lh[10+3]= lh[10+3] - coeff*lh[10+1];
  lh[15+3]= lh[15+3] - coeff*lh[15+1];
  lh[20+3]= lh[20+3] - coeff*lh[20+1];
  r[3]   = r[3]   - coeff*r[1];

  coeff = lh[5+4];
  lh[10+4]= lh[10+4] - coeff*lh[10+1];
  lh[15+4]= lh[15+4] - coeff*lh[15+1];
  lh[20+4]= lh[20+4] - coeff*lh[20+1];
  r[4]   = r[4]   - coeff*r[1];


  pivot = 1.00/lh[10+2];
  lh[15+2] = lh[15+2]*pivot;
  lh[20+2] = lh[20+2]*pivot;
  r[2]   = r[2]  *pivot;

  coeff = lh[10+0];
  lh[15+0]= lh[15+0] - coeff*lh[15+2];
  lh[20+0]= lh[20+0] - coeff*lh[20+2];
  r[0]   = r[0]   - coeff*r[2];

  coeff = lh[10+1];
  lh[15+1]= lh[15+1] - coeff*lh[15+2];
  lh[20+1]= lh[20+1] - coeff*lh[20+2];
  r[1]   = r[1]   - coeff*r[2];

  coeff = lh[10+3];
  lh[15+3]= lh[15+3] - coeff*lh[15+2];
  lh[20+3]= lh[20+3] - coeff*lh[20+2];
  r[3]   = r[3]   - coeff*r[2];

  coeff = lh[10+4];
  lh[15+4]= lh[15+4] - coeff*lh[15+2];
  lh[20+4]= lh[20+4] - coeff*lh[20+2];
  r[4]   = r[4]   - coeff*r[2];


  pivot = 1.00/lh[15+3];
  lh[20+3] = lh[20+3]*pivot;
  r[3]   = r[3]  *pivot;

  coeff = lh[15+0];
  lh[20+0]= lh[20+0] - coeff*lh[20+3];
  r[0]   = r[0]   - coeff*r[3];

  coeff = lh[15+1];
  lh[20+1]= lh[20+1] - coeff*lh[20+3];
  r[1]   = r[1]   - coeff*r[3];

  coeff = lh[15+2];
  lh[20+2]= lh[20+2] - coeff*lh[20+3];
  r[2]   = r[2]   - coeff*r[3];

  coeff = lh[15+4];
  lh[20+4]= lh[20+4] - coeff*lh[20+3];
  r[4]   = r[4]   - coeff*r[3];


  pivot = 1.00/lh[20+4];
  r[4]   = r[4]  *pivot;

  coeff = lh[20+0];
  r[0]   = r[0]   - coeff*r[4];

  coeff = lh[20+1];
  r[1]   = r[1]   - coeff*r[4];

  coeff = lh[20+2];
  r[2]   = r[2]   - coeff*r[4];

  coeff = lh[20+3];
  r[3]   = r[3]   - coeff*r[4];
}