/* bwrm99.c */
#include "bwrm.h"


int bwrm2ERROR(FILE *fp, char *s, double P, double T, double yv, double yl)
{
  
/*
  fprintf(fp, "%s at ", s);
  fprintf(fp, "P=%6.3fMPa, ", P/1e+6);
  fprintf(fp, "T=%6.1fC, ", T-T0bwrm);
  fprintf(fp, "YV=%6.2f, ", yv);
  fprintf(fp, "YL=%6.2f, ", yv);
  fprintf(fp, "\n");
*/
}

int 
bwrmSolveA22x_b(double A[2][2], double B[2], double X[2])
{
#define a A[0][0]
#define b A[0][1]
#define c A[1][0]
#define d A[1][1]
#define u B[0]
#define v B[1]
#define x X[0]
#define y X[1]
  double k=0;
  if ( a == 0.0){
	if (b == 0.0) return(1);
	y = u / b;
  }else{
	k = c/a*b - d;
	if (k == 0.0) return(1);
	y = (c/a*u - v) / k;
  }
  if ( a == 0.0 ){
	if (c == 0.0)  return(1);
	y = v / c;
  }else{
	x = (u - b * y ) / a;
  }
  return(0);
#undef a
#undef b
#undef c
#undef d
#undef u
#undef v
#undef x
#undef y
}