home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 9 / 09.iso / l / l048 / 1.ddi / BEZIER.HGH < prev    next >
Encoding:
Text File  |  1986-03-20  |  2.3 KB  |  69 lines
  1. (***********************************************************)
  2. (*                                                         *)
  3. (*                TURBO GRAPHIX version 1.06A              *)
  4. (*                                                         *)
  5. (*                 Bezier polynomial module                *)
  6. (*                  Module version  1.06A                  *)
  7. (*                                                         *)
  8. (*                  Copyright (C) 1985 by                  *)
  9. (*                  BORLAND International                  *)
  10. (*                                                         *)
  11. (***********************************************************)
  12.  
  13. procedure Bezier(A : PlotArray; MaxContrPoints : integer;
  14.                  var B : PlotArray; MaxIntPoints : integer);
  15. const
  16.   MaxControlPoints = 25;
  17. type
  18.   CombiArray = array[0..MaxControlPoints] of real;
  19. var
  20.   N : integer;
  21.   ContrPoint, IntPoint : integer;
  22.   T, SumX, SumY, Prod, DeltaT, Quot : real;
  23.   Combi : CombiArray;
  24.  
  25. begin
  26.   MaxContrPoints := MaxContrPoints - 1;
  27.   DeltaT := 1.0 / (MaxIntPoints - 1);
  28.   Combi[0] := 1;
  29.   Combi[MaxContrPoints] := 1;
  30.   for N := 0 to MaxContrPoints - 2 do
  31.     Combi[N + 1] := Combi[N] * (MaxContrPoints - N) / (N + 1);
  32.   for IntPoint := 1 to MaxIntPoints do
  33.   begin
  34.     T := (IntPoint - 1) * DeltaT;
  35.     if T <= 0.5 then
  36.       begin
  37.         Prod := 1.0 - T;
  38.         Quot := Prod;
  39.         for N := 1 to MaxContrPoints - 1 do
  40.           Prod := Prod * Quot;
  41.         Quot := T / Quot;
  42.         SumX := A[MaxContrPoints + 1, 1];
  43.         SumY := A[MaxContrPoints + 1, 2];
  44.         for N := MaxContrPoints downto 1 do
  45.         begin
  46.           SumX := Combi[N - 1] * A[N, 1] + Quot * SumX;
  47.           SumY := Combi[N - 1] * A[N, 2] + Quot * SumY;
  48.         end;
  49.       end
  50.     else
  51.       begin
  52.         Prod := T;
  53.         Quot := Prod;
  54.         for N := 1 to MaxContrPoints - 1 do
  55.           Prod := Prod * Quot;
  56.         Quot := (1 - T) / Quot;
  57.         SumX := A[1, 1];
  58.         SumY := A[1, 2];
  59.         for N := 1 to MaxContrPoints do
  60.         begin
  61.           SumX := Combi[N] * A[N + 1, 1] + Quot * SumX;
  62.           SumY := Combi[N] * A[N + 1, 2] + Quot * SumY;
  63.         end;
  64.       end;
  65.     B[IntPoint, 1] := SumX * Prod;
  66.     B[IntPoint, 2] := SumY * Prod;
  67.   end;
  68. end; { Bezier }
  69.