Anonymous

checking for d
 penetratedp()
  a,b,c:vector3(0, 0, 0)vector3(1, 0, 0)vector3(1, 1, 0)
  X:vector3(0.5, 0.5, 6.12313e-17)
  u:vector3(0, 0, 1)
   range_in()
    before: a,b,c:vector3(0, 0, 0)vector3(1, 0, 0)vector3(1, 1, 0)
    X, u:vector3(0.5, 0.5, 6.12313e-17)vector3(0, 0, 1)
     calculating A..
    X, u:vector3(0.5, 0.5, 6.12313e-17)vector3(0, 0, 1)
     calculating B..
    X, u:vector3(0.5, 0.5, 6.12313e-17)vector3(0, 0, 1)
     calculating C..
    X, u:vector3(0.5, 0.5, 6.12313e-17)vector3(0, 0, 1)
     calculating D..
    X, u:vector3(0.5, 0.5, 6.12313e-17)vector3(0, 0, 1)
    after: a,b,c:vector3(0, 0, 0)vector3(1, 0, 0)vector3(1, 1, 0)
    A==1,B==1.22463e-16,C==-0.25,D==1
   range_in()
    before: a,b,c:vector3(0, 0, 0)vector3(1, 0, 0)vector3(1, 1, 0)
    X, u:vector3(0.5, 0.5, 6.12313e-17)vector3(0, 0, 1)
     calculating A..
    X, u:vector3(0.5, 0.5, 6.12313e-17)vector3(0, 0, 1)
     calculating B..
    X, u:vector3(0.5, 0.5, 6.12313e-17)vector3(0, 0, 1)
     calculating C..
    X, u:vector3(0.5, 0.5, 6.12313e-17)vector3(0, 0, 1)
     calculating D..
    X, u:vector3(0.5, 0.5, 6.12313e-17)vector3(0, 0, 1)
    after: a,b,c:vector3(0, 0, 0)vector3(1, 0, 0)vector3(1, 1, 0)
    A==1,B==1.22463e-16,C==-0.25,D==1
  a,b,c:vector3(0, 0, 0)vector3(1, 0, 0)vector3(1, 1, 0)
  r1:-0.5, 0.5
  r2:-0.5, 0.5
  t_on:-6.12313e-17
result:yes
checking for 7adouble
 penetratedp()
  a,b,c:vector3(0(a), 0(a), 0(a))vector3(1(a), 0(a), 0(a))vector3(1(a), 1(a), 0(a))
  X:vector3(0.5(a), 0.5(a), 6.12323e-17(a))
  u:vector3(0(a), 0(a), 1(a))
   range_in()
    before: a,b,c:vector3(0(a), 0(a), 0(a))vector3(1(a), 0(a), 0(a))vector3(1(a), 1(a), 0(a))
    X, u:vector3(0.5(a), 0.5(a), 6.12323e-17(a))vector3(0(a), 0(a), 1(a))
     calculating A..
    X, u:vector3(0.5(a), 0.5(a), 6.12323e-17(a))vector3(0(a), 0(a), 1(a))
     calculating B..
    X, u:vector3(0.5(a), 0.5(a), 6.12323e-17(a))vector3(0(a), 0(a), 1(a))
     calculating C..
    X, u:vector3(0.5(a), 0.5(a), 6.12323e-17(a))vector3(0(a), 0(a), 1(a))
     calculating D..
    X, u:vector3(0.5(a), 0.5(a), 6.12323e-17(a))vector3(0(a), 0(a), 1(a))
    after: a,b,c:vector3(0(a), 0(a), 0(a))vector3(1(a), 0(a), 0(a))vector3(1(a), 1(a), 0(a))
    A==1(a),B==1.22465e-16(a),C==-0.25(a),D==1(a)
   range_in()
    before: a,b,c:vector3(0(a), 0(a), 0(a))vector3(1(a), 0(a), 0(a))vector3(1(a), 1(a), 0(a))
    X, u:vector3(0.5(a), 0.5(a), 6.12323e-17(a))vector3(0(a), 0(a), 1(a))
     calculating A..
    X, u:vector3(0.5(a), 0.5(a), 6.12323e-17(a))vector3(0(a), 0(a), 0(a))
     calculating B..
    X, u:vector3(0.5(a), 0.5(a), 6.12323e-17(a))vector3(0(a), 0(a), 1(a))
     calculating C..
    X, u:vector3(0.5(a), 0.5(a), 6.12323e-17(a))vector3(0(a), 0(a), 0(a))
     calculating D..
    X, u:vector3(0.5(a), 0.5(a), 6.12323e-17(a))vector3(0(a), 0(a), 0(a))
    after: a,b,c:vector3(0(a), 0(a), 0(a))vector3(1(a), 0(a), 0(a))vector3(1(a), 1(a), 0(a))
    A==0(a),B==1.22465e-16(a),C==-0.25(a),D==1.49976e-32(a)
  a,b,c:vector3(0(a), 0(a), 0(a))vector3(1(a), 0(a), 0(a))vector3(1(a), 1(a), 0(a))
  r1:-0.5(a), 0.5(a)
  r2:-inf(a), -nan(a)
  t_on:-6.12323e-17(a)
result:no