Vectors

call.psl

deffun vector(x, y, z) { x: x, y: y, z: z } in

deffun sub(v1, v2)
  vector( -(v1.x, v2.x), -(v1.y, v2.y), -(v1.z, v2.z))
in

deffun abs(x)
  if <(x, 0) then -(0, x) else x
in

deffun mult(x, y)
  defvar yi = 0 in
  defvar r = 0 in {
    for(yi = 0; <(yi, abs(y)); yi++) {
      r = +(r, x);
    };
    abs(r);
  }
in

deffun dotP(v1, v2)
  +(mult(v1.x, v2.x), mult(v1.y, v2.y), mult(v1.z, v2.z));
in

deffun dist2(v1, v2)
  defvar td = sub(v1, v2) in {
    dotP(td, td)
  }
in {
  print(dotP(vector(1,2,3), vector(4,5,6)));
  print(" ");
  print(dist2(vector(1,2,3), vector(4,5,6)));
  "";
}