// test_inside.java   compile    javac -cp . test_inside.java
//                    execute    java  -cp . test_inside > test_inside_java.out

public class test_inside  // point between lines 
{
  public test_inside()
  {
    double p1x, p1y, p2x, p2y, p3x, p3y, p4x, p4y, pix, piy;
    double xmax, xmin, ymax, ymin, dx, dy, x, y;
    int nx = 9;
    int ny = 9;
    double ddx = 0.01;
    double ddy = 0.01;
    int nline = 3;
    int found = 0;
    double mat[][] = new double[2][3];
    double p[] = new double[2];
    double pinside[] = new double[2];
    double line[][] = new double[3][4];
    boolean up = false;
    boolean down = false;
    boolean left = false;
    boolean right = false;
    
    System.out.println("test_inside.java running");
    System.out.println("find if point inside lines");
    System.out.println("point cheched below, above, left, right of some line");
    line[0][0] =  4.0; // lp1 triangle 
    line[0][1] =  1.0;
    line[0][2] =  1.0; // lp2
    line[0][3] = 12.0;
    line[1][0] =  1.0; // lp2
    line[1][1] = 12.0;
    line[1][2] = 12.0; // lp3
    line[1][3] =  9.0;
    line[2][0] = 12.0; // lp3
    line[2][1] =  9.0;
    line[2][2] =  4.0; // lp1
    line[2][3] =  1.0;
    xmin =  2.0;
    xmax = 11.0;
    ymin =  2.0;
    ymax = 11.0;
    dx = (xmax-xmin)/(double)nx;
    dy = (ymax-ymin)/(double)ny;
    System.out.println("nx="+nx);
    System.out.println("xmax="+xmax+"  xmin="+xmin+"  dx="+dx);
    System.out.println("ny="+ny);
    System.out.println("ymax="+ymax+"  ymin="+ymin+"  dy="+dy);
    for(int i=0; i<nline; i++)
    {
      System.out.println("line["+i+"i][0]="+line[i][0]+"  line[i][1] ="+line[i][1]);
      System.out.println("line["+i+"i][2]="+line[i][2]+"  line[i][3] ="+line[i][3]);
      System.out.println(" ");
    } // end i

    x = xmin;
    for(int i=0; i<=nx; i++)
    {
      y = ymin;
      for(int j=0; j<=ny; j++)
      {
	up = false;
        down = false;
        left = false;
        right = false;

        System.out.println("point inside? x="+x+"  y="+y);
	found = 0; // need found == 4 for inside
        for(int m=0; m<nline; m++) // check x up 
	{
	  p1x = line[m][0];
	  p1y = line[m][1];
	  p2x = line[m][2];
	  p2y = line[m][3];
          // checking up
	  if(((x>=p1x && x<=p2x) || (x>=p2x && x<=p1x)) && (p1y>y || p2y>y))
	  {
	    piy = p1y+((p2y-p1y)*(x-p1x))/(p2x-p1x);
	    System.out.println("up at m="+m+ " piy="+piy);
            if(piy>y)
            {
              up = true;
              System.out.println("up at m="+m);
	    }
	  }
          // checking down
	  if(((x>=p1x && x<=p2x) || (x>=p2x && x<=p1x)) && (p1y<y || p2y<y))
	  {
	    piy = p1y+((p2y-p1y)*(x-p1x))/(p2x-p1x);
	    System.out.println("down at m="+m+ " piy="+piy);
            if(piy<y)
            {
              down = true;
              System.out.println("down at m="+m);
	    }
	  }
          // checking left
          if(((y>=p1y && y<=p2y) || (y>=p2y && y<=p1y)) && (p1x<x || p2x<x))
	  {
            pix = p1x+((p2x-p1x)*(y-p1y))/(p2y-p1y);
            System.out.println("left at m="+m+ " pix="+pix);
            if(pix<x)
	    {
	      left = true;
	      System.out.println("left at m="+m);
	    }
	  }
          // checking right
          if(((y>=p1y && y<=p2y) || (y>=p2y && y<=p1y)) && (p1x>x || p2x>x))
	  {
            pix = p1x+((p2x-p1x)*(y-p1y))/(p2y-p1y);
            System.out.println("right at m="+m+ " pix="+pix);
	    if(pix>x)
            {		
	      right = true;
	      System.out.println("right at m="+m);
	    }
	  }
	  if(up && down && left && right)
          {
            System.out.println("point is inside");
            System.out.println(" ");
            break; // exit m loop
	  }
	} // end m
        System.out.println(" ");
	y = y + dy;
      } // end j
      System.out.println(" ");
      x = x + dx;
    } // end y

    System.out.println("test_inside.java finished");
  } // end test_inside

  void inside(double p1x, double p1y, double p2x, double p2y,
	    double p3x, double p3y, double p4x, double p4y,
	    double pinside[]) // pix, piy
  {
    double d12x, d12y, d34x, d34y, pix, piy;
    double p[] = {0.0, 0.0, 0.0};
    double mat[][] = new double[2][3];
    pinside[0] = 0.0;
    pinside[1] = 0.0;
    d12x = p2x-p1x;
    d12y = p2y-p1y;
    d34x = p4x-p3x;
    d34y = p4y-p3y;
    // System.out.println("d12x="+d12x+"  d12y="+d12y);
    // System.out.println("d34x="+d34x+"  d34y="+d34y);
    // line p1 to p2  (x-p1x)/d12x = (y-p1y)/d12y 
    // x/d12x - y/d12y = p1x/d12x - p1y/d12y  
    // line p3 to p4  (x-p3x)/d34x = (y-p3y)/d34y
    // x/d34x - y/d34y = p3x/d34x - p3y/d34y
    // solve two equations in two unknowns
    if(Math.abs(d12x)<1.0e-6 || Math.abs(d12y)<1.0e-6) {return;}
    if(Math.abs(d34x)<1.0e-6 || Math.abs(d34y)<1.0e-6) {return;}
    mat[0][0] = 1.0/d12x;
    mat[0][1] = -1.0/d12y;
    mat[0][2] = p1x/d12x - p1y/d12y;
    mat[1][0] = 1.0/d34x;
    mat[1][1] = -1.0/d34y;
    mat[1][2] = p3x/d34x - p3y/d34y; 
    new simeq(2,mat,p);
    pix = p[0];
    piy = p[1];
    pinside[0] = pix;
    pinside[1] = piy;
  } // end inside

  public static void main (String[] args)
  {
    new test_inside();
  }
} // end test_inside