Sorry your browser is not supported!

Rem Project: Matrix Mouse Selection
Rem Created: 22.01.2003 20:50:47

`*******************************************************************************
`*                    By: DMiTR0S   dmitr0s@newmail.ru                         *
`*******************************************************************************

`*******************************************************************************
`*                            Free License Agreement.                          *
`*                                                                             *
`*   The user/developer agrees to include a visible credit to "Dmitri Kuschev" *
`*   within your programs documentation and or web site.                       *
`*                                                                             *
`*   The user/developer agrees to NOT offer the this source code for Sale.     *
`*                                                                             *
`*   Providing the user/Developer agrees in full to the previous Free          *
`*   License Agreement clauses, the user/developer may use this source         *
`*   code FREELY in all their products, commercial or otherwise.               *
`*                                                                             *
`*******************************************************************************
Set Display Mode 1024, 768, 32

Autocam Off

Type TPoint
   X# As Float
   Y# As Float
   Z# As Float
EndType

Type MatrixPlane
   mX As Float
   mY As Float
   mZ As Float
   mD As Float

   R1X As Float
   R1Y As Float
   R1Z As Float
   R1D As Integer

   R2X As Float
   R2Y As Float
   R2Z As Float
   R2D As Float

   R3X As Float
   R3Y As Float
   R3Z As Float
   R3D As Float

   mXP As Integer
   mZP As Integer

   pX As Float
   pY As Float
   pZ As Float
EndType

Global MatView=3
Global Vector=1
Global DestPoint=2
Global VecA=3
Global VecB=4
Global VecC=5
Global VecD=6
Global Vec1=7
Global Vec2=8
Global MaxY# As Float
Global SelPoly As Integer
Global Epsilon#=0.000001
Global IntersectionPoint As TPoint
Global MNearestX
Global MNearestZ
Global Selected

FOV#=61.930722089872112754955664756038
Set Camera FOV FOV#
MaxY#=Tan(FOV#/2.0)

r=Make Matrix4(MatView)
r=Make Vector3(Vector)
r=Make Vector3(DestPoint)
r=Make Vector3(VecA)
r=Make Vector3(VecB)
r=Make Vector3(VecC)
r=Make Vector3(VecD)
r=Make Vector3(Vec1)
r=Make Vector3(Vec2)


MatrixWidth#=200.0
MatrixHeight#=200.0
MatrixXSegments=50
MatrixZSegments=50
MatrixXSegmentSize#=MatrixWidth#/(MatrixXSegments*1.0)
MatrixZSegmentSize#=MatrixHeight#/(MatrixZSegments*1.0)
MVCount=(MatrixXSegments+1)*(MatrixZSegments+1)

Make Matrix 1, MatrixWidth#, MatrixHeight#, MatrixXSegments, MatrixZSegments
Set Matrix 1, 1, 1, 1, 1, 1, 1, 1
Randomize Matrix 1, 5
Update Matrix 1

Dim MatrixPlaneList(MatrixXSegments*MatrixZSegments*2) As MatrixPlane
Dim MatrixVerticles(MVCount) As TPoint
UpdatePlanesOfMatrix(1, MatrixXSegments, MatrixZSegments, MatrixXSegmentSize#, MatrixZSegmentSize#, MatrixXSegments*MatrixZSegments)

Pointer=1

Make Object Sphere Pointer, 0.5
Color Object Pointer, RGB(255,0,0)

Make Object Sphere 12, 1
Color Object 12, RGB(255,255,128)


Sync Rate 0
Sync On
Position Camera 50,20,-10
Do

   Mouse2DTo3D(Pointer)

   If Mouseclick()=1
      Mouse2DTo3D(Pointer)
      FindIntersection(Pointer)
   EndIf

   Control Camera Using Arrowkeys 0, 2, 2

   Text 10, 10, "FPS="+Str$(Screen FPS())
   If Object Position X(1)>Matrix Position X(1) And Object Position X(1)<Matrix Position X(1)+MatrixWidth# And Object Position Z(1)>Matrix Position Z(1) And Object Position Z(1)<Matrix Position Z(1)+MatrixHeight#
      Text 10, 70, "Matrix tile X="+Str$(Int((Object Position X(1)-Matrix Position X(1))/MatrixXSegmentSize#))+" Z="+Str$(Int((Object Position Z(1)-Matrix Position Z(1))/MatrixZSegmentSize#))
      Text 10, 90, "Position of a pointer on the matrix:"
      Text 10, 110, "X="+Str$(Object Position X(1)-Matrix Position X(1))
      Text 10, 130, "Y="+Str$(Object Position Y(1)-Matrix Position Y(1))
      Text 10, 150, "Z="+Str$(Object Position Z(1)-Matrix Position Z(1))
      FindNearestPoint(1, MatrixXSegments, MatrixZSegments, MatrixXSegmentSize#, MatrixZSegmentSize#)
   EndIf
   If Inkey$()="+"
      Set Matrix Height 1, MNearestX, MNearestZ, Get Matrix Height(1, MNearestX, MNearestZ)+1.0
      Update Matrix 1
      UpdatePlanesOfMatrix(1, MatrixXSegments, MatrixZSegments, MatrixXSegmentSize#, MatrixZSegmentSize#, MatrixXSegments*MatrixZSegments)
      Position Object 12, Matrix Position X(1)+MNearestX*MatrixXSegmentSize#, Matrix Position Y(1)+Get Matrix Height(1, MNearestX, MNearestZ), Matrix Position Z(1)+MNearestZ*MatrixZSegmentSize#
   EndIf
   If Inkey$()="-"
      Set Matrix Height 1, MNearestX, MNearestZ, Get Matrix Height(1, MNearestX, MNearestZ)-1.0
      Update Matrix 1
      UpdatePlanesOfMatrix(1, MatrixXSegments, MatrixZSegments, MatrixXSegmentSize#, MatrixZSegmentSize#, MatrixXSegments*MatrixZSegments)
      Position Object 12, Matrix Position X(1)+MNearestX*MatrixXSegmentSize#, Matrix Position Y(1)+Get Matrix Height(1, MNearestX, MNearestZ), Matrix Position Z(1)+MNearestZ*MatrixZSegmentSize#
   EndIf
   Sync
Loop

End

Function UpdatePlanesOfMatrix(Matrix, MXSegments, MZSegments, MXSegmentSize#, MZSegmentSize#, Offset)
   For SX=1 to MXSegments
      For SZ=1 To MZSegments
         Num=MXSegments*(SZ-1)+SX

         PolyY#=Matrix Position Y(Matrix)

         Poly1X1#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
         Poly1Y1#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX-1, SZ-1)
         Poly1Z1#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#

         Poly1X2#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
         Poly1Y2#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX-1, SZ)
         Poly1Z2#=Matrix Position Z(Matrix)+SZ*MZSegmentSize#

         Poly1X3#=Matrix Position X(Matrix)+SX*MXSegmentSize#
         Poly1Y3#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX, SZ)
         Poly1Z3#=Matrix Position Z(Matrix)+SZ*MZSegmentSize#

         Poly2X1#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
         Poly2Y1#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX-1, SZ-1)
         Poly2Z1#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#

         Poly2X2#=Matrix Position X(Matrix)+SX*MXSegmentSize#
         Poly2Y2#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX, SZ)
         Poly2Z2#=Matrix Position Z(Matrix)+SZ*MZSegmentSize#

         Poly2X3#=Matrix Position X(Matrix)+SX*MXSegmentSize#
         Poly2Y3#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX, SZ-1)
         Poly2Z3#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#

         `Find matrix plane
         Set Vector3 VecA, Poly1X1#, Poly1Y#, Poly1Z1#
         Set Vector3 VecB, Poly1X2#, Poly1Y#, Poly1Z2#
         Set Vector3 VecC, Poly1X3#, Poly1Y#, Poly1Z3#

         Subtract Vector3 Vec1, VecB, VecA
         Subtract Vector3 Vec2, VecC, VecA
         Cross Product Vector3 Vec1, Vec1, Vec2
         Normalize Vector3 Vec1, Vec1
         Set Vector3 VecD, X Vector3(Vec1), Y Vector3(Vec1), Z Vector3(Vec1)

         `Poly1 plane
         Set Vector3 VecA, Poly1X1#, Poly1Y1#, Poly1Z1#
         Set Vector3 VecB, Poly1X2#, Poly1Y2#, Poly1Z2#
         Set Vector3 VecC, Poly1X3#, Poly1Y3#, Poly1Z3#

         Subtract Vector3 Vec1, VecB, VecA
         Subtract Vector3 Vec2, VecC, VecA
         Cross Product Vector3 Vec1, Vec1, Vec2
         Normalize Vector3 Vec1, Vec1

         MatrixPlaneList(Num).mX=X Vector3(Vec1)
         MatrixPlaneList(Num).mY=Y Vector3(Vec1)
         MatrixPlaneList(Num).mZ=Z Vector3(Vec1)
         MatrixPlaneList(Num).mD=Dot Product Vector3(Vec1, VecA)*-1.0

         `Poly1 Edge1 plane
         Set Vector3 VecA, Poly1X1#, PolyY#, Poly1Z1#
         Set Vector3 VecB, Poly1X2#, PolyY#, Poly1Z2#
         Set Vector3 VecC, Poly1X2#+X Vector3(VecD), PolyY#+Y Vector3(VecD), Poly1Z2#+Z Vector3(VecD)

         Subtract Vector3 Vec1, VecB, VecA
         Subtract Vector3 Vec2, VecC, VecA
         Cross Product Vector3 Vec1, Vec1, Vec2

         MatrixPlaneList(Num).R1X=X Vector3(Vec1)
         MatrixPlaneList(Num).R1Y=Y Vector3(Vec1)
         MatrixPlaneList(Num).R1Z=Z Vector3(Vec1)
         MatrixPlaneList(Num).R1D=Dot Product Vector3(Vec1, VecA)*-1.0


         `Poly1 Edge2 plane
         Set Vector3 VecA, Poly1X2#, PolyY#, Poly1Z2#
         Set Vector3 VecB, Poly1X3#, PolyY#, Poly1Z3#
         Set Vector3 VecC, Poly1X3#+X Vector3(VecD), PolyY#+Y Vector3(VecD), Poly1Z3#+Z Vector3(VecD)

         Subtract Vector3 Vec1, VecB, VecA
         Subtract Vector3 Vec2, VecC, VecA
         Cross Product Vector3 Vec1, Vec1, Vec2

         MatrixPlaneList(Num).R2X=X Vector3(Vec1)
         MatrixPlaneList(Num).R2Y=Y Vector3(Vec1)
         MatrixPlaneList(Num).R2Z=Z Vector3(Vec1)
         MatrixPlaneList(Num).R2D=Dot Product Vector3(Vec1, VecA)*-1.0

         `Poly1 Edge3 plane
         Set Vector3 VecA, Poly1X3#, PolyY#, Poly1Z3#
         Set Vector3 VecB, Poly1X1#, PolyY#, Poly1Z1#
         Set Vector3 VecC, Poly1X1#+X Vector3(VecD), PolyY#+Y Vector3(VecD), Poly1Z1#+Z Vector3(VecD)

         Subtract Vector3 Vec1, VecB, VecA
         Subtract Vector3 Vec2, VecC, VecA
         Cross Product Vector3 Vec1, Vec1, Vec2

         MatrixPlaneList(Num).R3X=X Vector3(Vec1)
         MatrixPlaneList(Num).R3Y=Y Vector3(Vec1)
         MatrixPlaneList(Num).R3Z=Z Vector3(Vec1)
         MatrixPlaneList(Num).R3D=Dot Product Vector3(Vec1, VecA)*-1.0

         MatrixPlaneList(Num).mXP=SX
         MatrixPlaneList(Num).mZP=SZ

         MatrixPlaneList(Num).pX=Poly1X1#
         MatrixPlaneList(Num).pY=Poly1Y1#
         MatrixPlaneList(Num).pZ=Poly1Z1#



         `Find Poly2 plane
         Set Vector3 VecA, Poly2X1#, Poly2Y1#, Poly2Z1#
         Set Vector3 VecB, Poly2X2#, Poly2Y2#, Poly2Z2#
         Set Vector3 VecC, Poly2X3#, Poly2Y3#, Poly2Z3#

         Subtract Vector3 Vec1, VecB, VecA
         Subtract Vector3 Vec2, VecC, VecA
         Cross Product Vector3 Vec1, Vec1, Vec2
         Normalize Vector3 Vec1, Vec1

         MatrixPlaneList(Num+Offset).mX=X Vector3(Vec1)
         MatrixPlaneList(Num+Offset).mY=Y Vector3(Vec1)
         MatrixPlaneList(Num+Offset).mZ=Z Vector3(Vec1)
         MatrixPlaneList(Num+Offset).mD=Dot Product Vector3(Vec1, VecA)*-1.0

         `Poly2 Edge1 plane
         Set Vector3 VecA, Poly2X1#, PolyY#, Poly2Z1#
         Set Vector3 VecB, Poly2X2#, PolyY#, Poly2Z2#
         Set Vector3 VecC, Poly2X2#+X Vector3(VecD), PolyY#+Y Vector3(VecD), Poly2Z2#+Z Vector3(VecD)

         Subtract Vector3 Vec1, VecB, VecA
         Subtract Vector3 Vec2, VecC, VecA
         Cross Product Vector3 Vec1, Vec1, Vec2

         MatrixPlaneList(Num+Offset).R1X=X Vector3(Vec1)
         MatrixPlaneList(Num+Offset).R1Y=Y Vector3(Vec1)
         MatrixPlaneList(Num+Offset).R1Z=Z Vector3(Vec1)
         MatrixPlaneList(Num+Offset).R1D=Dot Product Vector3(VecA, Vec1)*-1.0

         `Poly2 Edge2 plane
         Set Vector3 VecA, Poly2X2#, PolyY#, Poly2Z2#
         Set Vector3 VecB, Poly2X3#, PolyY#, Poly2Z3#
         Set Vector3 VecC, Poly2X3#+X Vector3(VecD), PolyY#+Y Vector3(VecD), Poly2Z3#+Z Vector3(VecD)

         Subtract Vector3 Vec1, VecB, VecA
         Subtract Vector3 Vec2, VecC, VecA
         Cross Product Vector3 Vec1, Vec1, Vec2

         MatrixPlaneList(Num+Offset).R2X=X Vector3(Vec1)
         MatrixPlaneList(Num+Offset).R2Y=Y Vector3(Vec1)
         MatrixPlaneList(Num+Offset).R2Z=Z Vector3(Vec1)
         MatrixPlaneList(Num+Offset).R2D=Dot Product Vector3(Vec1, VecA)*-1.0

         `Poly2 Edge3 plane
         Set Vector3 VecA, Poly2X3#, PolyY#, Poly2Z3#
         Set Vector3 VecB, Poly2X1#, PolyY#, Poly2Z1#
         Set Vector3 VecC, Poly2X1#+X Vector3(VecD), PolyY#+Y Vector3(VecD), Poly2Z1#+Z Vector3(VecD)

         Subtract Vector3 Vec1, VecB, VecA
         Subtract Vector3 Vec2, VecC, VecA
         Cross Product Vector3 Vec1, Vec1, Vec2

         MatrixPlaneList(Num+Offset).R3X=X Vector3(Vec1)
         MatrixPlaneList(Num+Offset).R3Y=Y Vector3(Vec1)
         MatrixPlaneList(Num+Offset).R3Z=Z Vector3(Vec1)
         MatrixPlaneList(Num+Offset).R3D=Dot Product Vector3(Vec1, VecA)*-1.0

         MatrixPlaneList(Num+Offset).mXP=SX
         MatrixPlaneList(Num+Offset).mZP=SZ

         MatrixPlaneList(Num+Offset).pX=Poly2X1#
         MatrixPlaneList(Num+Offset).pY=Poly2Y1#
         MatrixPlaneList(Num+Offset).pZ=Poly2Z1#


      Next SZ

   Next SX

EndFunction


Function Mouse2DTo3D(Obj)

   NMouseX#=1.0-2.0*(Mousex()*1.0)/(Screen Width()*1.0)
   NMouseY#=1.0-2.0*(MouseY()*1.0)/(Screen Height()*1.0)

   Aspectratio#=1.3333333333333333
   CameraspaceY#=NMouseY#*MaxY#
   CameraspaceX#=(-1.0*NMouseX#*Aspectratio#)*MaxY#
   CameraspaceZ#=1.0

   View Matrix4 MatView
   r#=Inverse Matrix4(MatView, MatView)

   Set Vector3 Vector, 1000.0*CameraspaceX#, 1000.0*CameraspaceY#, 1000.0*CameraspaceZ#
   Transform Coords Vector3 DestPoint, Vector, MatView

   Position Object Obj, X Vector3(DestPoint), Y Vector3(DestPoint), Z Vector3(DestPoint)


EndFunction

Function FindIntersection(Obj)
spX#=Camera Position X()
spY#=Camera Position Y()
spZ#=Camera Position Z()
dpX#=Object Position X(Obj)
dpY#=Object Position Y(Obj)
dpZ#=Object Position Z(Obj)

NearDist#=100000.0

For i=1  To Array Count(MatrixPlaneList())
   If GetIntersect(spX#, spY#, spZ#, dpX#, dpY#, dpZ#, MatrixPlaneList(i).pX, MatrixPlaneList(i).pY, MatrixPlaneList(i).pZ, MatrixPlaneList(i).mX, MatrixPlaneList(i).mY, MatrixPlaneList(i).mZ, MatrixPlaneList(i).mD)
      If InPoly(i)
         Set Vector3 VecA, spX#-IntersectionPoint.X#, spY#-IntersectionPoint.Y#, spZ#-IntersectionPoint.Z#
         Dist#=Length Vector3(VecA)
         If Dist#<NearDist#
            NearDist#=Dist#
            SelPoly=i
            Position Object Obj, IntersectionPoint.X#, IntersectionPoint.Y#, IntersectionPoint.Z#
         EndIf
      EndIf
   EndIf
Next i


EndFunction

Function GetIntersect(V1X#, V1Y#, V1Z#, V2X#, V2Y#, V2Z#, PVX#, PVY#, PVZ#, PX#, PY#, PZ#, PD#)
   Set Vector3 VecA, V2X#-V1X#, V2Y#-V1Y#, V2Z#-V1Z#
   Set Vector3 VecB, PX#, PY#, PZ#
   LineLength#=Dot Product Vector3(VecA, VecB)

   If ABS(LineLength#)<Epsilon# Then ExitFunction 0

   Set Vector3 VecC, PVX#-V1X#, PVY#-V1Y#, PVZ#-V1Z#

   Dist#=Dot Product Vector3(VecC, VecB)
   Percentage#=Dist#/LineLength#

   If Percentage#<0.0 Or Percentage#>1.0 Then ExitFunction 0

   Multiply Vector3 VecA, Percentage#

   IntersectionPoint.X#=V1X#+X Vector3(VecA)
   IntersectionPoint.Y#=V1Y#+Y Vector3(VecA)
   IntersectionPoint.Z#=V1Z#+Z Vector3(VecA)

EndFunction 1

Function DistanceToPlane(PlaneX#, PlaneY#, PlaneZ#, PlaneD#, VertexX#, VertexY#, VertexZ#)
   Set Vector3 Vec1, PlaneX#, PlaneY#, PlaneZ#
   Set Vector3 Vec2, VertexX#, VertexY#, VertexZ#
   Dist#=Dot Product Vector3(Vec1, Vec2)+PlaneD#
EndFunction Dist#

Function ClassifyPoint(PlaneX#, PlaneY#, PlaneZ#, PlaneD#, VertexX#, VertexY#, VertexZ#)
   Dist#=DistanceToPlane(PlaneX#, PlaneY#, PlaneZ#, PlaneD#, VertexX#, VertexY#, VertexZ#)
   If Dist#>Epsilon#
      ExitFunction "BACK"
   Else
      ExitFunction "FRONT"
   EndIf
EndFunction "ONPLANE"

Function InPoly(Poly)
   If ClassifyPoint(MatrixPlaneList(Poly).R1X, MatrixPlaneList(Poly).R1Y, MatrixPlaneList(Poly).R1Z, MatrixPlaneList(Poly).R1D, IntersectionPoint.X#, IntersectionPoint.Y#, IntersectionPoint.Z#)="BACK" Then ExitFunction 0
   If ClassifyPoint(MatrixPlaneList(Poly).R2X, MatrixPlaneList(Poly).R2Y, MatrixPlaneList(Poly).R2Z, MatrixPlaneList(Poly).R2D, IntersectionPoint.X#, IntersectionPoint.Y#, IntersectionPoint.Z#)="BACK" Then ExitFunction 0
   If ClassifyPoint(MatrixPlaneList(Poly).R3X, MatrixPlaneList(Poly).R3Y, MatrixPlaneList(Poly).R3Z, MatrixPlaneList(Poly).R3D, IntersectionPoint.X#, IntersectionPoint.Y#, IntersectionPoint.Z#)="BACK" Then ExitFunction 0
EndFunction 1

Function FindNearestPoint(Matrix, MXSegments, MZSegments, MXSegmentSize#, MZSegmentSize#)
   NearestDist#=1000.0
   Selected=0
   For SX=1 to MXSegments+1
      For SZ=1 To MZSegments+1
        X#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
        Z#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#
        Dist#=Sqrt((X#-Object Position X(1))*(X#-Object Position X(1))+(Z#-Object Position Z(1))*(Z#-Object Position Z(1)))
        If Dist#<NearestDist#
           NearestDist#=Dist#
           MatrixNearestX=SX-1
           MatrixNearestZ=SZ-1
           Selected=1
        EndIf
      Next SZ
   Next SX
   If Selected=1
      Position Object 12, Matrix Position X(Matrix)+(MatrixNearestX)*MXSegmentSize#, Matrix Position Y(Matrix)+Get Matrix Height(Matrix, MatrixNearestX, MatrixNearestZ), Matrix Position Z(Matrix)+(MatrixNearestZ)*MZSegmentSize#
      MNearestX=MatrixNearestX
      MNearestX=MatrixNearestX
      MNearestZ=MatrixNearestZ
   Else
      MNearestX=0
      MNearestZ=0
   EndIf
EndFunction

`*******************************************************************************
`*                      Mouse selection of matrix tile                         *
`*                           Created: 01/25/2003                               *
`*                     By: DMiTR0S   dmitr0s@newmail.ru                        *
`*******************************************************************************

`*******************************************************************************
`*                            Free License Agreement.                          *
`*                                                                             *
`*   The user/developer agrees to include a visible credit to "Dmitri Kuschev" *
`*   within your programs documentation and or web site.                       *
`*                                                                             *
`*   The user/developer agrees to NOT offer the this source code for Sale.     *
`*                                                                             *
`*   Providing the user/Developer agrees in full to the previous Free          *
`*   License Agreement clauses, the user/developer may use this source         *
`*   code FREELY in all their products, commercial or otherwise.               *
`*                                                                             *
`*******************************************************************************
Set Display Mode 1024, 768, 32

Autocam off
Position Camera 50.0, 25.0, -60.0

Global ViewMatrix
Global ProjectionMatrix
Global Ray As Integer
Global Vector As Integer
Global DestPoint As Integer
Global NearestPoint As Integer
Global MaxY# As Float
Global Distance# As Float
Global IntersectionPointX# As Float
Global IntersectionPointY# As Float
Global IntersectionPointZ# As Float
Global MNearestX
Global MNearestZ
Global Selected

Vector=1
DestPoint=2
NearestPoint=3
ViewMatrix=101
ProjectionMatrix=102



r=Make Matrix4(ViewMatrix)
r=Make Matrix4(ProjectionMatrix)
r=Make Vector3(Vector)
r=Make Vector3(DestPoint)
r=Make Vector3(NearestPoint)

MatrixClone=1
Pointer=2
Pointer1=3
Pointer2=4
Pointer3=5
Pointer4=6
Matrix=1

MatrixWidth#=100.0
MatrixHeight#=100.0
MatrixXSegments=20
MatrixZSegments=13
MatrixXSegmentSize#=MatrixWidth#/(MatrixXSegments*1.0)
MatrixZSegmentSize#=MatrixHeight#/(MatrixZSegments*1.0)
MVCount=(MatrixXSegments+1)*(MatrixZSegments+1)

Make Matrix Matrix, MatrixWidth#, MatrixHeight#, MatrixXSegments, MatrixZSegments
Randomize Matrix Matrix, 3
Update Matrix Matrix

Make Object Sphere Pointer1, 0.5
Color Object Pointer1, RGB(255,0,0)
Make Object Sphere Pointer2, 0.5
Color Object Pointer2, RGB(255,0,0)
Make Object Sphere Pointer3, 0.5
Color Object Pointer3, RGB(255,0,0)
Make Object Sphere Pointer4, 0.5
Color Object Pointer4, RGB(255,0,0)

Make Object Sphere 10, 1.0
Color Object 10, RGB(0,255,0)

Make Object Sphere 11, 0.5
Color Object 11, RGB(0,0,255)

MakeObjectFromMatrix(Matrix, MatrixClone, MatrixXSegments, MatrixZSegments, MatrixXSegmentSize#, MatrixZSegmentSize#, MatrixXSegments+1)

Sync Rate 0
Sync On
Do

   If Inkey$()="w" Then Move Camera 1

   If Inkey$()="s" Then Move Camera -1

   If Upkey()Then Rotate Camera WrapValue(Camera Angle X()-2), Camera Angle Y(), Camera Angle Z()

   If Downkey() Then Rotate Camera WrapValue(Camera Angle X()+2), Camera Angle Y(), Camera Angle Z()

   If Rightkey() Then Rotate Camera Camera Angle X(), WrapValue(Camera Angle Y()+3), Camera Angle Z()

   If Leftkey() Then Rotate Camera Camera Angle X(), WrapValue(Camera Angle Y()-3), Camera Angle Z()

   If MouseClick()=1
      Set Matrix Height 1, MNearestX, MNearestZ, Get Matrix Height(1, MNearestX, MNearestZ)+1.0
      Update Matrix 1
      MakeObjectFromMatrix(Matrix, MatrixClone, MatrixXSegments, MatrixZSegments, MatrixXSegmentSize#, MatrixZSegmentSize#, MatrixXSegments+1)
      Position Object 10, Matrix Position X(1)+MNearestX*MatrixXSegmentSize#, Matrix Position Y(1)+Get Matrix Height(1, MNearestX, MNearestZ), Matrix Position Z(1)+MNearestZ*MatrixZSegmentSize#
   EndIf

   If MouseClick()=2
      Set Matrix Height 1, MNearestX, MNearestZ, Get Matrix Height(1, MNearestX, MNearestZ)-1.0
      Update Matrix 1
      MakeObjectFromMatrix(Matrix, MatrixClone, MatrixXSegments, MatrixZSegments, MatrixXSegmentSize#, MatrixZSegmentSize#, MatrixXSegments+1)
      Position Object 10, Matrix Position X(1)+MNearestX*MatrixXSegmentSize#, Matrix Position Y(1)+Get Matrix Height(1, MNearestX, MNearestZ), Matrix Position Z(1)+MNearestZ*MatrixZSegmentSize#
   EndIf

   GetMatrixCloneIntersection(MatrixClone)
   If GetDistanceToObject()>0.0 Then Position Object 11, GetIntersectionPointX(), GetIntersectionPointY(), GetIntersectionPointZ()


   MatrixX=Int((GetIntersectionPointX()-Matrix Position X(Matrix))/MatrixXSegmentSize#)
   MatrixZ=Int((GetIntersectionPointZ()-Matrix Position Z(Matrix))/MatrixZSegmentSize#)
   MatrixYHeight#=GetIntersectionPointY()-Matrix Position Y(Matrix)


   If MatrixX=>0 And MatrixX<MatrixXSegments And MatrixZ=>0 And MatrixZ<MatrixZSegments And GetDistanceToObject()>0
      Position Object Pointer1, Matrix Position X(Matrix)+MatrixX*MatrixXSegmentSize#, Matrix Position Y(Matrix)+Get Matrix Height(Matrix, MatrixX, MatrixZ), Matrix Position Z(Matrix)+MatrixZ*MatrixZSegmentSize#
      Position Object Pointer2, Matrix Position X(Matrix)+(MatrixX+1)*MatrixXSegmentSize#, Matrix Position Y(Matrix)+Get Matrix Height(Matrix, MatrixX+1, MatrixZ), Matrix Position Z(Matrix)+MatrixZ*MatrixZSegmentSize#
      Position Object Pointer3, Matrix Position X(Matrix)+(MatrixX+1)*MatrixXSegmentSize#, Matrix Position Y(Matrix)+Get Matrix Height(Matrix, MatrixX+1, MatrixZ+1), Matrix Position Z(Matrix)+(MatrixZ+1)*MatrixZSegmentSize#
      Position Object Pointer4, Matrix Position X(Matrix)+MatrixX*MatrixXSegmentSize#, Matrix Position Y(Matrix)+Get Matrix Height(Matrix, MatrixX, MatrixZ+1), Matrix Position Z(Matrix)+(MatrixZ+1)*MatrixZSegmentSize#
   EndIf

   Text 10, 10, "Code by: Dmitri Kuschev"
   Text 10, 30, "E-mail: dmitr0s@newmail.ru"
   Text 10, 50, "FPS "+Str$(Screen FPS())
   Text 10, 70, "Keys 'w'/'s' - Forward/Backward.  Arrowkeys - Rotate camera"

   Text 10, 110, "Distance to the matrix "+Str$(GetDistanceToObject())
   Text 10, 130, "World intersection point:"
   Text 10, 150, "X = "+Str$(GetIntersectionPointX())
   Text 10, 170, "Y = "+Str$(GetIntersectionPointY())
   Text 10, 190, "Z = "+Str$(GetIntersectionPointZ())
   Text 10, 210, "Matrix tile: X="+Str$(MatrixX)+" Z="+Str$(MatrixZ)
   Text 10, 230, "Intersection point height = "+Str$(MatrixYHeight#)
   FindNearestPoint(1, MatrixXSegments, MatrixZSegments, MatrixXSegmentSize#, MatrixZSegmentSize#)
   Position Object 10, X Vector3(NearestPoint),  Y Vector3(NearestPoint),  Z Vector3(NearestPoint)
   Sync
Loop

Function MakeObjectFromMatrix(Matrix, Object, ,MXSegments, MZSegments, MXSegmentSize#, MZSegmentSize#, VIR)
   VerticlesCount=(MXSegments)*(MZSegments)*6
   FVFFormat=338
   FVFSize=36
   If Memblock Exist(1) Then Delete Memblock 1
   Make Memblock 1, 12+VerticlesCount*FVFSize
   Write Memblock Dword 1, 0, FVFFormat
   Write Memblock Dword 1, 4, FVFSize
   Write Memblock Dword 1, 8, VerticlesCount
   For SX=1 to MXSegments
      For SZ=1 To MZSegments
         Num=MXSegments*(SZ-1)+SX
         Poly1X1#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
         Poly1Y1#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX-1, SZ-1)
         Poly1Z1#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#

         Poly1X2#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
         Poly1Y2#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX-1, SZ)
         Poly1Z2#=Matrix Position Z(Matrix)+SZ*MZSegmentSize#

         Poly1X3#=Matrix Position X(Matrix)+SX*MXSegmentSize#
         Poly1Y3#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX, SZ)
         Poly1Z3#=Matrix Position Z(Matrix)+SZ*MZSegmentSize#

         Poly2X1#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
         Poly2Y1#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX-1, SZ-1)
         Poly2Z1#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#

         Poly2X2#=Matrix Position X(Matrix)+SX*MXSegmentSize#
         Poly2Y2#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX, SZ)
         Poly2Z2#=Matrix Position Z(Matrix)+SZ*MZSegmentSize#

         Poly2X3#=Matrix Position X(Matrix)+SX*MXSegmentSize#
         Poly2Y3#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX, SZ-1)
         Poly2Z3#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#

         Write Memblock Float 1, 12+216*(Num-1), Poly1X1#
         Write Memblock Float 1, 16+216*(Num-1), Poly1Y1#
         Write Memblock Float 1, 20+216*(Num-1), Poly1Z1#
         Write Memblock Float 1, 24+216*(Num-1), 0.0
         Write Memblock Float 1, 28+216*(Num-1), 0.0
         Write Memblock Float 1, 32+216*(Num-1), -1.0
         Write Memblock Dword 1, 36+216*(Num-1), -1
         Write Memblock Float 1, 40+216*(Num-1), 0.0
         Write Memblock Float 1, 44+216*(Num-1), 1.0

         Write Memblock Float 1, 48+216*(Num-1), Poly1X2#
         Write Memblock Float 1, 52+216*(Num-1), Poly1Y2#
         Write Memblock Float 1, 56+216*(Num-1), Poly1Z2#
         Write Memblock Float 1, 60+216*(Num-1), 0.0
         Write Memblock Float 1, 64+216*(Num-1), 0.0
         Write Memblock Float 1, 68+216*(Num-1), -1.0
         Write Memblock Dword 1, 72+216*(Num-1), -1
         Write Memblock Float 1, 76+216*(Num-1), 0.0
         Write Memblock Float 1, 80+216*(Num-1), 0.0

         Write Memblock Float 1, 84+216*(Num-1), Poly1X3#
         Write Memblock Float 1, 88+216*(Num-1), Poly1Y3#
         Write Memblock Float 1, 92+216*(Num-1), Poly1Z3#
         Write Memblock Float 1, 96+216*(Num-1), 0.0
         Write Memblock Float 1, 100+216*(Num-1), 0.0
         Write Memblock Float 1, 104+216*(Num-1), -1.0
         Write Memblock Dword 1, 108+216*(Num-1), -1
         Write Memblock Float 1, 112+216*(Num-1), 1.0
         Write Memblock Float 1, 116+216*(Num-1), 1.0

         Write Memblock Float 1, 120+216*(Num-1), Poly2X1#
         Write Memblock Float 1, 124+216*(Num-1), Poly2Y1#
         Write Memblock Float 1, 128+216*(Num-1), Poly2Z1#
         Write Memblock Float 1, 132+216*(Num-1), 0.0
         Write Memblock Float 1, 136+216*(Num-1), 0.0
         Write Memblock Float 1, 140+216*(Num-1), -1.0
         Write Memblock Dword 1, 144+216*(Num-1), -1
         Write Memblock Float 1, 148+216*(Num-1), 1.0
         Write Memblock Float 1, 152+216*(Num-1), 0.0

         Write Memblock Float 1, 156+216*(Num-1), Poly2X2#
         Write Memblock Float 1, 160+216*(Num-1), Poly2Y2#
         Write Memblock Float 1, 164+216*(Num-1), Poly2Z2#
         Write Memblock Float 1, 168+216*(Num-1), 0.0
         Write Memblock Float 1, 172+216*(Num-1), 0.0
         Write Memblock Float 1, 176+216*(Num-1), -1.0
         Write Memblock Dword 1, 180+216*(Num-1), -1
         Write Memblock Float 1, 184+216*(Num-1), 1.0
         Write Memblock Float 1, 188+216*(Num-1), 1.0

         Write Memblock Float 1, 192+216*(Num-1), Poly2X3#
         Write Memblock Float 1, 196+216*(Num-1), Poly2Y3#
         Write Memblock Float 1, 200+216*(Num-1), Poly2Z3#
         Write Memblock Float 1, 204+216*(Num-1), 0.0
         Write Memblock Float 1, 208+216*(Num-1), 0.0
         Write Memblock Float 1, 212+216*(Num-1), -1.0
         Write Memblock Dword 1, 216+216*(Num-1), -1
         Write Memblock Float 1, 220+216*(Num-1), 0.0
         Write Memblock Float 1, 224+216*(Num-1), 1.0

      Next SZ

   Next SX

  Make Mesh From Memblock 1, 1
  If Object Exist(Object) Then Delete Object Object
  Make Object Object, 1, 22
  Set Object Light 1, 0
EndFunction

Function GetMatrixCloneIntersection(Obj)

   X#=2.0*(Mousex()*1.0)/(Screen Width()*1.0)-1.0
   Y#=1.0-2.0*(MouseY()*1.0)/(Screen Height()*1.0)
   Set Vector3 DestPoint, X#, Y#, 1.0
   View Matrix4 ViewMatrix
   Projection Matrix4 ProjectionMatrix
   r=Inverse Matrix4(ViewMatrix, ViewMatrix)
   r=Inverse Matrix4(ProjectionMatrix, ProjectionMatrix)
   Transform Coords Vector3 DestPoint, DestPoint, ProjectionMatrix
   Normalize Vector3 DestPoint, DestPoint
   Multiply Vector3 DestPoint, 1000.0
   Transform Coords Vector3 DestPoint, DestPoint, ViewMatrix

   Distance#=Intersect Object(Obj, Camera Position X(), Camera Position Y(), Camera Position Z(), X Vector3(DestPoint), Y Vector3(DestPoint), Z Vector3(DestPoint))

   Set Vector3 DestPoint, (X Vector3(DestPoint)-Camera Position X()), (Y Vector3(DestPoint)-Camera Position Y()), (Z Vector3(DestPoint)-Camera Position Z())
   Normalize Vector3 DestPoint, DestPoint
   Multiply Vector3 DestPoint, Distance#
   Set Vector3 DestPoint, X Vector3(DestPoint)+Camera Position X(), Y Vector3(DestPoint)+Camera Position Y(), Z Vector3(DestPoint)+Camera Position Z()

   If Distance#>0.0
      IntersectionPointX#=X Vector3(DestPoint)
      IntersectionPointY#=Y Vector3(DestPoint)
      IntersectionPointZ#=Z Vector3(DestPoint)
   Else
      IntersectionPointX#=0.0
      IntersectionPointY#=0.0
      IntersectionPointZ#=0.0
   EndIf

EndFunction

Function FindNearestPoint(Matrix, MXSegments, MZSegments, MXSegmentSize#, MZSegmentSize#)
   NearestDist#=1000.0
   Selected=0
   For SX=1 to MXSegments+1
      For SZ=1 To MZSegments+1
        X#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
        Z#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#
        Dist#=Sqrt((X#-GetIntersectionPointX())*(X#-GetIntersectionPointX())+(Z#-GetIntersectionPointZ())*(Z#-GetIntersectionPointZ()))
        If Dist#<NearestDist#
           NearestDist#=Dist#
           MatrixNearestX=SX-1
           MatrixNearestZ=SZ-1
           Selected=1
        EndIf
      Next SZ
   Next SX
   If Selected=1
      Set Vector3 NearestPoint, Matrix Position X(Matrix)+(MatrixNearestX)*MXSegmentSize#, Matrix Position Y(Matrix)+Get Matrix Height(Matrix, MatrixNearestX, MatrixNearestZ), Matrix Position Z(Matrix)+(MatrixNearestZ)*MZSegmentSize#
      MNearestX=MatrixNearestX
      MNearestZ=MatrixNearestZ
   Else
      MNearestX=0
      MNearestZ=0
   EndIf
EndFunction

Function GetIntersectionPointX()
EndFunction IntersectionPointX#

Function GetIntersectionPointY()
EndFunction IntersectionPointY#

Function GetIntersectionPointZ()
EndFunction IntersectionPointZ#

Function GetDistanceToObject()
EndFunction Distance#


Delete Memblock 1

`*******************************************************************************
`*                      Mouse selection of matrix tile                         *
`*                           Created: 01/25/2003                               *
`*                     By: DMiTR0S   dmitr0s@newmail.ru                        *
`*******************************************************************************

`*******************************************************************************
`*                            Free License Agreement.                          *
`*                                                                             *
`*   The user/developer agrees to include a visible credit to "Dmitri Kuschev" *
`*   within your programs documentation and or web site.                       *
`*                                                                             *
`*   The user/developer agrees to NOT offer the this source code for Sale.     *
`*                                                                             *
`*   Providing the user/Developer agrees in full to the previous Free          *
`*   License Agreement clauses, the user/developer may use this source         *
`*   code FREELY in all their products, commercial or otherwise.               *
`*                                                                             *
`*******************************************************************************
Set Display Mode 1024, 768, 32

Autocam off
Position Camera 50.0, 25.0, -60.0
Randomize Timer()

Global Distance# As Float
Global IntersectionPointX# As Float
Global IntersectionPointY# As Float
Global IntersectionPointZ# As Float

MatrixClone=1
Matrix=1

MatrixWidth#=100.0
MatrixHeight#=100.0
MatrixXSegments=20
MatrixZSegments=13
MatrixXSegmentSize#=MatrixWidth#/(MatrixXSegments*1.0)
MatrixZSegmentSize#=MatrixHeight#/(MatrixZSegments*1.0)

Make Matrix Matrix, MatrixWidth#, MatrixHeight#, MatrixXSegments, MatrixZSegments
Randomize Matrix Matrix, 5
Update Matrix Matrix
MakeObjectFromMatrix(Matrix, MatrixClone, MatrixXSegments, MatrixZSegments, MatrixXSegmentSize#, MatrixZSegmentSize#, MatrixXSegments+1)

Make Object Sphere 10, 0.5
Color Object 10, RGB(0,255,0)

Sync Rate 0
Sync On
Do

   If Inkey$()="w" Then Move Camera 1

   If Inkey$()="s" Then Move Camera -1

   If Upkey()Then Rotate Camera WrapValue(Camera Angle X()-2), Camera Angle Y(), Camera Angle Z()

   If Downkey() Then Rotate Camera WrapValue(Camera Angle X()+2), Camera Angle Y(), Camera Angle Z()

   If Rightkey() Then Rotate Camera Camera Angle X(), WrapValue(Camera Angle Y()+3), Camera Angle Z()

   If Leftkey() Then Rotate Camera Camera Angle X(), WrapValue(Camera Angle Y()-3), Camera Angle Z()

   If MouseClick()=1
      GetMatrixCloneIntersection(MatrixClone)
      If GetDistanceToObject()>0.0
         Show Object 10
         Position Object 10, GetIntersectionPointX(), GetIntersectionPointY(), GetIntersectionPointZ()
      EndIf
      Text 10, 110, "Distance to the matrix "+Str$(GetDistanceToObject())
      Text 10, 130, "Intersection point:"
      Text 10, 150, "X = "+Str$(GetIntersectionPointX())
      Text 10, 170, "Y = "+Str$(GetIntersectionPointY())
      Text 10, 190, "Z = "+Str$(GetIntersectionPointZ())
   Else
      Hide Object 10
   EndIf

   MatrixX=Int((GetIntersectionPointX()-Matrix Position X(Matrix))/MatrixXSegmentSize#)
   MatrixZ=Int((GetIntersectionPointZ()-Matrix Position Z(Matrix))/MatrixZSegmentSize#)
   MatrixYHeight#=GetIntersectionPointY()-Matrix Position Y(Matrix)

   Text 10, 10, "Code by: Dmitri Kuschev"
   Text 10, 30, "E-mail: dmitr0s@newmail.ru"
   Text 10, 50, "FPS "+Str$(Screen FPS())
   Text 10, 70, "Keys 'w'/'s' - Forward/Backward.  Arrowkeys - Rotate camera"
   Sync
Loop

Function MakeObjectFromMatrix(Matrix, Object, ,MXSegments, MZSegments, MXSegmentSize#, MZSegmentSize#, VIR)
   VerticlesCount=(MXSegments)*(MZSegments)*6
   FVFFormat=338
   FVFSize=36
   If Memblock Exist(1) Then Delete Memblock 1
   Make Memblock 1, 12+VerticlesCount*FVFSize
   Write Memblock Dword 1, 0, FVFFormat
   Write Memblock Dword 1, 4, FVFSize
   Write Memblock Dword 1, 8, VerticlesCount
   For SX=1 to MXSegments
      For SZ=1 To MZSegments
         Num=MXSegments*(SZ-1)+SX
         Poly1X1#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
         Poly1Y1#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX-1, SZ-1)
         Poly1Z1#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#

         Poly1X2#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
         Poly1Y2#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX-1, SZ)
         Poly1Z2#=Matrix Position Z(Matrix)+SZ*MZSegmentSize#

         Poly1X3#=Matrix Position X(Matrix)+SX*MXSegmentSize#
         Poly1Y3#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX, SZ)
         Poly1Z3#=Matrix Position Z(Matrix)+SZ*MZSegmentSize#

         Poly2X1#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
         Poly2Y1#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX-1, SZ-1)
         Poly2Z1#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#

         Poly2X2#=Matrix Position X(Matrix)+SX*MXSegmentSize#
         Poly2Y2#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX, SZ)
         Poly2Z2#=Matrix Position Z(Matrix)+SZ*MZSegmentSize#

         Poly2X3#=Matrix Position X(Matrix)+SX*MXSegmentSize#
         Poly2Y3#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX, SZ-1)
         Poly2Z3#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#

         Write Memblock Float 1, 12+216*(Num-1), Poly1X1#
         Write Memblock Float 1, 16+216*(Num-1), Poly1Y1#
         Write Memblock Float 1, 20+216*(Num-1), Poly1Z1#
         Write Memblock Float 1, 24+216*(Num-1), 0.0
         Write Memblock Float 1, 28+216*(Num-1), 0.0
         Write Memblock Float 1, 32+216*(Num-1), -1.0
         Write Memblock Dword 1, 36+216*(Num-1), -1
         Write Memblock Float 1, 40+216*(Num-1), 0.0
         Write Memblock Float 1, 44+216*(Num-1), 1.0

         Write Memblock Float 1, 48+216*(Num-1), Poly1X2#
         Write Memblock Float 1, 52+216*(Num-1), Poly1Y2#
         Write Memblock Float 1, 56+216*(Num-1), Poly1Z2#
         Write Memblock Float 1, 60+216*(Num-1), 0.0
         Write Memblock Float 1, 64+216*(Num-1), 0.0
         Write Memblock Float 1, 68+216*(Num-1), -1.0
         Write Memblock Dword 1, 72+216*(Num-1), -1
         Write Memblock Float 1, 76+216*(Num-1), 0.0
         Write Memblock Float 1, 80+216*(Num-1), 0.0

         Write Memblock Float 1, 84+216*(Num-1), Poly1X3#
         Write Memblock Float 1, 88+216*(Num-1), Poly1Y3#
         Write Memblock Float 1, 92+216*(Num-1), Poly1Z3#
         Write Memblock Float 1, 96+216*(Num-1), 0.0
         Write Memblock Float 1, 100+216*(Num-1), 0.0
         Write Memblock Float 1, 104+216*(Num-1), -1.0
         Write Memblock Dword 1, 108+216*(Num-1), -1
         Write Memblock Float 1, 112+216*(Num-1), 1.0
         Write Memblock Float 1, 116+216*(Num-1), 1.0

         Write Memblock Float 1, 120+216*(Num-1), Poly2X1#
         Write Memblock Float 1, 124+216*(Num-1), Poly2Y1#
         Write Memblock Float 1, 128+216*(Num-1), Poly2Z1#
         Write Memblock Float 1, 132+216*(Num-1), 0.0
         Write Memblock Float 1, 136+216*(Num-1), 0.0
         Write Memblock Float 1, 140+216*(Num-1), -1.0
         Write Memblock Dword 1, 144+216*(Num-1), -1
         Write Memblock Float 1, 148+216*(Num-1), 1.0
         Write Memblock Float 1, 152+216*(Num-1), 0.0

         Write Memblock Float 1, 156+216*(Num-1), Poly2X2#
         Write Memblock Float 1, 160+216*(Num-1), Poly2Y2#
         Write Memblock Float 1, 164+216*(Num-1), Poly2Z2#
         Write Memblock Float 1, 168+216*(Num-1), 0.0
         Write Memblock Float 1, 172+216*(Num-1), 0.0
         Write Memblock Float 1, 176+216*(Num-1), -1.0
         Write Memblock Dword 1, 180+216*(Num-1), -1
         Write Memblock Float 1, 184+216*(Num-1), 1.0
         Write Memblock Float 1, 188+216*(Num-1), 1.0

         Write Memblock Float 1, 192+216*(Num-1), Poly2X3#
         Write Memblock Float 1, 196+216*(Num-1), Poly2Y3#
         Write Memblock Float 1, 200+216*(Num-1), Poly2Z3#
         Write Memblock Float 1, 204+216*(Num-1), 0.0
         Write Memblock Float 1, 208+216*(Num-1), 0.0
         Write Memblock Float 1, 212+216*(Num-1), -1.0
         Write Memblock Dword 1, 216+216*(Num-1), -1
         Write Memblock Float 1, 220+216*(Num-1), 0.0
         Write Memblock Float 1, 224+216*(Num-1), 1.0

      Next SZ

   Next SX

  Make Mesh From Memblock 1, 1
  If Object Exist(Object) Then Delete Object Object
  Make Object Object, 1, 22
  Set Object Light 1, 0
EndFunction

Function GetMatrixCloneIntersection(Obj)
   DestPoint=2
   ViewMatrix=101
   ProjectionMatrix=102
   r=Make Matrix4(ViewMatrix)
   r=Make Matrix4(ProjectionMatrix)
   r=Make Vector3(DestPoint)
   X#=2.0*(Mousex()*1.0)/(Screen Width()*1.0)-1.0
   Y#=1.0-2.0*(MouseY()*1.0)/(Screen Height()*1.0)
   Set Vector3 DestPoint, X#, Y#, 1.0
   View Matrix4 ViewMatrix
   Projection Matrix4 ProjectionMatrix
   r=Inverse Matrix4(ViewMatrix, ViewMatrix)
   r=Inverse Matrix4(ProjectionMatrix, ProjectionMatrix)
   Transform Coords Vector3 DestPoint, DestPoint, ProjectionMatrix
   Normalize Vector3 DestPoint, DestPoint
   Multiply Vector3 DestPoint, 1000.0
   Transform Coords Vector3 DestPoint, DestPoint, ViewMatrix

   Distance#=Intersect Object(Obj, Camera Position X(), Camera Position Y(), Camera Position Z(), X Vector3(DestPoint), Y Vector3(DestPoint), Z Vector3(DestPoint))

   Set Vector3 DestPoint, (X Vector3(DestPoint)-Camera Position X()), (Y Vector3(DestPoint)-Camera Position Y()), (Z Vector3(DestPoint)-Camera Position Z())
   Normalize Vector3 DestPoint, DestPoint
   Multiply Vector3 DestPoint, Distance#
   Set Vector3 DestPoint, X Vector3(DestPoint)+Camera Position X(), Y Vector3(DestPoint)+Camera Position Y(), Z Vector3(DestPoint)+Camera Position Z()

   If Distance#>0.0
      IntersectionPointX#=X Vector3(DestPoint)
      IntersectionPointY#=Y Vector3(DestPoint)
      IntersectionPointZ#=Z Vector3(DestPoint)
   Else
      IntersectionPointX#=0.0
      IntersectionPointY#=0.0
      IntersectionPointZ#=0.0
   EndIf
EndFunction

Function GetIntersectionPointX()
EndFunction IntersectionPointX#

Function GetIntersectionPointY()
EndFunction IntersectionPointY#

Function GetIntersectionPointZ()
EndFunction IntersectionPointZ#

Function GetDistanceToObject()
EndFunction Distance#


Delete Memblock 1

sync on
sync rate
make matrix
randomize matrix
position camera
point camera
etc....etc...

do
if upkey()=1 then move camera blah blah
if mouseclick()=1 then convert2dmouseto3d()
print x#
print y#
print z#

sync
loop

function convert2dmouseto3d()
complex maths stuff here
blah blah blah
blah blah blah
blah blah blah
x# = .....
y# = .....
z# = .....
endfunction

`*******************************************************************************
`*                      Mouse selection of matrix tile                         *
`*                           Created: 01/25/2003                               *
`*                     By: DMiTR0S   dmitr0s@newmail.ru                        *
`*******************************************************************************

`*******************************************************************************
`*                            Free License Agreement.                          *
`*                                                                             *
`*   The user/developer agrees to include a visible credit to "Dmitri Kuschev" *
`*   within your programs documentation and or web site.                       *
`*                                                                             *
`*   The user/developer agrees to NOT offer the this source code for Sale.     *
`*                                                                             *
`*   Providing the user/Developer agrees in full to the previous Free          *
`*   License Agreement clauses, the user/developer may use this source         *
`*   code FREELY in all their products, commercial or otherwise.               *
`*                                                                             *
`*******************************************************************************
Set Display Mode 1024, 768, 32

Autocam off
Position Camera 50.0, 25.0, -60.0
Randomize Timer()

Global Distance# As Float
Global IntersectionPointX# As Float
Global IntersectionPointY# As Float
Global IntersectionPointZ# As Float

MatrixClone=1
Matrix=1

MatrixWidth#=100.0
MatrixHeight#=100.0
MatrixXSegments=20
MatrixZSegments=13
MatrixXSegmentSize#=MatrixWidth#/(MatrixXSegments*1.0)
MatrixZSegmentSize#=MatrixHeight#/(MatrixZSegments*1.0)

Make Matrix Matrix, MatrixWidth#, MatrixHeight#, MatrixXSegments, MatrixZSegments
Position Matrix Matrix, RND(10), RND(10), RND(10)
Randomize Matrix Matrix, 5
Update Matrix Matrix
MakeObjectFromMatrix(Matrix, MatrixClone, MatrixXSegments, MatrixZSegments, MatrixXSegmentSize#, MatrixZSegmentSize#, MatrixXSegments+1)

Make Object Sphere 10, 0.5
Color Object 10, RGB(0,255,0)

Sync Rate 0
Sync On
Do

   If Inkey$()="w" Then Move Camera 1

   If Inkey$()="s" Then Move Camera -1

   If Upkey()Then Rotate Camera WrapValue(Camera Angle X()-2), Camera Angle Y(), Camera Angle Z()

   If Downkey() Then Rotate Camera WrapValue(Camera Angle X()+2), Camera Angle Y(), Camera Angle Z()

   If Rightkey() Then Rotate Camera Camera Angle X(), WrapValue(Camera Angle Y()+3), Camera Angle Z()

   If Leftkey() Then Rotate Camera Camera Angle X(), WrapValue(Camera Angle Y()-3), Camera Angle Z()

   If MouseClick()=1
      GetMatrixCloneIntersection(MatrixClone)
      If GetDistanceToObject()>0.0
         Show Object 10
         Position Object 10, GetIntersectionPointX(), GetIntersectionPointY(), GetIntersectionPointZ()
      EndIf
      Text 10, 110, "Distance to the matrix "+Str$(GetDistanceToObject())
      Text 10, 130, "Position of the cursor on the matrix in world coordinates:"
      Text 10, 150, "X = "+Str$(GetIntersectionPointX())
      Text 10, 170, "Y = "+Str$(GetIntersectionPointY())
      Text 10, 190, "Z = "+Str$(GetIntersectionPointZ())
      Text 10, 210, "Position of the cursor on the matrix in matrix coordinates:"
      Text 10, 230, "X="+Str$(GetIntersectionPointX()-Matrix Position X(1))
      Text 10, 250, "Y="+Str$(GetIntersectionPointY()-Matrix Position Y(1))
      Text 10, 270, "Z="+Str$(GetIntersectionPointZ()-Matrix Position Z(1))
   Else
      Hide Object 10
   EndIf

   MatrixX=Int((GetIntersectionPointX()-Matrix Position X(Matrix))/MatrixXSegmentSize#)
   MatrixZ=Int((GetIntersectionPointZ()-Matrix Position Z(Matrix))/MatrixZSegmentSize#)
   MatrixYHeight#=GetIntersectionPointY()-Matrix Position Y(Matrix)

   Text 10, 10, "Code by: Dmitri Kuschev"
   Text 10, 30, "E-mail: dmitr0s@newmail.ru"
   Text 10, 50, "FPS "+Str$(Screen FPS())
   Text 10, 70, "Keys 'w'/'s' - Forward/Backward.  Arrowkeys - Rotate camera"
   Sync
Loop

Function MakeObjectFromMatrix(Matrix, Object, ,MXSegments, MZSegments, MXSegmentSize#, MZSegmentSize#, VIR)
   VerticlesCount=(MXSegments)*(MZSegments)*6
   FVFFormat=338
   FVFSize=36
   If Memblock Exist(1) Then Delete Memblock 1
   Make Memblock 1, 12+VerticlesCount*FVFSize
   Write Memblock Dword 1, 0, FVFFormat
   Write Memblock Dword 1, 4, FVFSize
   Write Memblock Dword 1, 8, VerticlesCount
   For SX=1 to MXSegments
      For SZ=1 To MZSegments
         Num=MXSegments*(SZ-1)+SX
         Poly1X1#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
         Poly1Y1#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX-1, SZ-1)
         Poly1Z1#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#

         Poly1X2#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
         Poly1Y2#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX-1, SZ)
         Poly1Z2#=Matrix Position Z(Matrix)+SZ*MZSegmentSize#

         Poly1X3#=Matrix Position X(Matrix)+SX*MXSegmentSize#
         Poly1Y3#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX, SZ)
         Poly1Z3#=Matrix Position Z(Matrix)+SZ*MZSegmentSize#

         Poly2X1#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
         Poly2Y1#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX-1, SZ-1)
         Poly2Z1#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#

         Poly2X2#=Matrix Position X(Matrix)+SX*MXSegmentSize#
         Poly2Y2#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX, SZ)
         Poly2Z2#=Matrix Position Z(Matrix)+SZ*MZSegmentSize#

         Poly2X3#=Matrix Position X(Matrix)+SX*MXSegmentSize#
         Poly2Y3#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX, SZ-1)
         Poly2Z3#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#

         Write Memblock Float 1, 12+216*(Num-1), Poly1X1#
         Write Memblock Float 1, 16+216*(Num-1), Poly1Y1#
         Write Memblock Float 1, 20+216*(Num-1), Poly1Z1#
         Write Memblock Float 1, 24+216*(Num-1), 0.0
         Write Memblock Float 1, 28+216*(Num-1), 0.0
         Write Memblock Float 1, 32+216*(Num-1), -1.0
         Write Memblock Dword 1, 36+216*(Num-1), -1
         Write Memblock Float 1, 40+216*(Num-1), 0.0
         Write Memblock Float 1, 44+216*(Num-1), 1.0

         Write Memblock Float 1, 48+216*(Num-1), Poly1X2#
         Write Memblock Float 1, 52+216*(Num-1), Poly1Y2#
         Write Memblock Float 1, 56+216*(Num-1), Poly1Z2#
         Write Memblock Float 1, 60+216*(Num-1), 0.0
         Write Memblock Float 1, 64+216*(Num-1), 0.0
         Write Memblock Float 1, 68+216*(Num-1), -1.0
         Write Memblock Dword 1, 72+216*(Num-1), -1
         Write Memblock Float 1, 76+216*(Num-1), 0.0
         Write Memblock Float 1, 80+216*(Num-1), 0.0

         Write Memblock Float 1, 84+216*(Num-1), Poly1X3#
         Write Memblock Float 1, 88+216*(Num-1), Poly1Y3#
         Write Memblock Float 1, 92+216*(Num-1), Poly1Z3#
         Write Memblock Float 1, 96+216*(Num-1), 0.0
         Write Memblock Float 1, 100+216*(Num-1), 0.0
         Write Memblock Float 1, 104+216*(Num-1), -1.0
         Write Memblock Dword 1, 108+216*(Num-1), -1
         Write Memblock Float 1, 112+216*(Num-1), 1.0
         Write Memblock Float 1, 116+216*(Num-1), 1.0

         Write Memblock Float 1, 120+216*(Num-1), Poly2X1#
         Write Memblock Float 1, 124+216*(Num-1), Poly2Y1#
         Write Memblock Float 1, 128+216*(Num-1), Poly2Z1#
         Write Memblock Float 1, 132+216*(Num-1), 0.0
         Write Memblock Float 1, 136+216*(Num-1), 0.0
         Write Memblock Float 1, 140+216*(Num-1), -1.0
         Write Memblock Dword 1, 144+216*(Num-1), -1
         Write Memblock Float 1, 148+216*(Num-1), 1.0
         Write Memblock Float 1, 152+216*(Num-1), 0.0

         Write Memblock Float 1, 156+216*(Num-1), Poly2X2#
         Write Memblock Float 1, 160+216*(Num-1), Poly2Y2#
         Write Memblock Float 1, 164+216*(Num-1), Poly2Z2#
         Write Memblock Float 1, 168+216*(Num-1), 0.0
         Write Memblock Float 1, 172+216*(Num-1), 0.0
         Write Memblock Float 1, 176+216*(Num-1), -1.0
         Write Memblock Dword 1, 180+216*(Num-1), -1
         Write Memblock Float 1, 184+216*(Num-1), 1.0
         Write Memblock Float 1, 188+216*(Num-1), 1.0

         Write Memblock Float 1, 192+216*(Num-1), Poly2X3#
         Write Memblock Float 1, 196+216*(Num-1), Poly2Y3#
         Write Memblock Float 1, 200+216*(Num-1), Poly2Z3#
         Write Memblock Float 1, 204+216*(Num-1), 0.0
         Write Memblock Float 1, 208+216*(Num-1), 0.0
         Write Memblock Float 1, 212+216*(Num-1), -1.0
         Write Memblock Dword 1, 216+216*(Num-1), -1
         Write Memblock Float 1, 220+216*(Num-1), 0.0
         Write Memblock Float 1, 224+216*(Num-1), 1.0

      Next SZ

   Next SX

  Make Mesh From Memblock 1, 1
  If Object Exist(Object) Then Delete Object Object
  Make Object Object, 1, 22
  Hide Object Object
EndFunction

Function GetMatrixCloneIntersection(Obj)
   DestPoint=2
   ViewMatrix=101
   ProjectionMatrix=102
   r=Make Matrix4(ViewMatrix)
   r=Make Matrix4(ProjectionMatrix)
   r=Make Vector3(DestPoint)
   X#=2.0*(Mousex()*1.0)/(Screen Width()*1.0)-1.0
   Y#=1.0-2.0*(MouseY()*1.0)/(Screen Height()*1.0)
   Set Vector3 DestPoint, X#, Y#, 1.0
   View Matrix4 ViewMatrix
   Projection Matrix4 ProjectionMatrix
   r=Inverse Matrix4(ViewMatrix, ViewMatrix)
   r=Inverse Matrix4(ProjectionMatrix, ProjectionMatrix)
   Transform Coords Vector3 DestPoint, DestPoint, ProjectionMatrix
   Normalize Vector3 DestPoint, DestPoint
   Multiply Vector3 DestPoint, 1000.0
   Transform Coords Vector3 DestPoint, DestPoint, ViewMatrix

   Distance#=Intersect Object(Obj, Camera Position X(), Camera Position Y(), Camera Position Z(), X Vector3(DestPoint), Y Vector3(DestPoint), Z Vector3(DestPoint))

   Set Vector3 DestPoint, (X Vector3(DestPoint)-Camera Position X()), (Y Vector3(DestPoint)-Camera Position Y()), (Z Vector3(DestPoint)-Camera Position Z())
   Normalize Vector3 DestPoint, DestPoint
   Multiply Vector3 DestPoint, Distance#
   Set Vector3 DestPoint, X Vector3(DestPoint)+Camera Position X(), Y Vector3(DestPoint)+Camera Position Y(), Z Vector3(DestPoint)+Camera Position Z()

   If Distance#>0.0
      IntersectionPointX#=X Vector3(DestPoint)
      IntersectionPointY#=Y Vector3(DestPoint)
      IntersectionPointZ#=Z Vector3(DestPoint)
   Else
      IntersectionPointX#=0.0
      IntersectionPointY#=0.0
      IntersectionPointZ#=0.0
   EndIf
EndFunction

Function GetIntersectionPointX()
EndFunction IntersectionPointX#

Function GetIntersectionPointY()
EndFunction IntersectionPointY#

Function GetIntersectionPointZ()
EndFunction IntersectionPointZ#

Function GetDistanceToObject()
EndFunction Distance#


Delete Memblock 1

View_Matrix=2
Projection_Matrix=3
VectorResult=1
r=Make Matrix4(View_Matrix)
r=Make Matrix4(Projection_Matrix)
r=Make Vector3(VectorResult)

2DTo3D(VectorResult, View_Matrix, Projection_Matrix, 1000.0)

Function 2DTo3D(Vector, ViewMatrix, ProjectionMatrix, Range#)
   X#=2.0*(Mousex()*1.0)/(Screen Width()*1.0)-1.0
   Y#=1.0-2.0*(MouseY()*1.0)/(Screen Height()*1.0)
   Set Vector3 Vector, X#, Y#, 1.0
   View Matrix4 ViewMatrix
   Projection Matrix4 ProjectionMatrix
   r=Inverse Matrix4(ViewMatrix, ViewMatrix)
   r=Inverse Matrix4(ProjectionMatrix, ProjectionMatrix)
   Transform Coords Vector3 Vector, Vector, ProjectionMatrix
   Normalize Vector3 Vector, Vector
   Multiply Vector3 Vector, Range#
   Transform Coords Vector3 Vector, Vector, ViewMatrix
   Projection Matrix4 ProjectionMatrix
EndFunction

View_Matrix=2
Projection_Matrix=3
VectorResult=1
r=Make Matrix4(View_Matrix)
r=Make Matrix4(Projection_Matrix)
r=Make Vector3(VectorResult)
Make Object Sphere 1, 1.0
Sync Rate 0
Sync On
Do
   Control Camera Using ArrowKeys 0, 5.0, 5.0

   2DTo3D(VectorResult, View_Matrix, Projection_Matrix, 100.0)

   Position Object 1, X Vector3(VectorResult), Y Vector3(VectorResult), Z Vector3(VectorResult)

   Text 10, 10, " X (X Vector3(VectorResult)) = "+Str$(X Vector3(VectorResult))
   Text 10, 30, " Y (Y Vector3(VectorResult)) = "+Str$(Y Vector3(VectorResult))
   Text 10, 50, " Z (Z Vector3(VectorResult)) = "+Str$(Z Vector3(VectorResult))

   Sync
Loop

Function 2DTo3D(Vector, ViewMatrix, ProjectionMatrix, Range#)
   X#=2.0*(Mousex()*1.0)/(Screen Width()*1.0)-1.0
   Y#=1.0-2.0*(MouseY()*1.0)/(Screen Height()*1.0)
   Set Vector3 Vector, X#, Y#, 1.0
   View Matrix4 ViewMatrix
   Projection Matrix4 ProjectionMatrix
   r=Inverse Matrix4(ViewMatrix, ViewMatrix)
   r=Inverse Matrix4(ProjectionMatrix, ProjectionMatrix)
   Transform Coords Vector3 Vector, Vector, ProjectionMatrix
   Normalize Vector3 Vector, Vector
   Multiply Vector3 Vector, Range#
   Transform Coords Vector3 Vector, Vector, ViewMatrix
   Projection Matrix4 ProjectionMatrix
EndFunction

sync on
sync rate 30

width=5000*2
depth=5000*2
xTiles=100
zTiles=100
ax#=45 : ay#=0 :ox#=width/2*tilex : oy#=0 : oz#=depth/2*tilez : r#=500 :FOV#=60

make matrix 1,width,depth,xTiles,zTiles
randomize matrix 1,150
update matrix 1

View_Matrix=2
Projection_Matrix=3
VectorResult=1
r=Make Matrix4(View_Matrix)
r=Make Matrix4(Projection_Matrix)
r=Make Vector3(VectorResult)

make object sphere 1,50

do
gosub camera
2DTo3D(VectorResult, View_Matrix, Projection_Matrix, 100.0)
xpos# = X Vector3(VectorResult)
ypos# = Y Vector3(VectorResult)
zpos# = Z Vector3(VectorResult)
 Text 10, 10, " X (X Vector3(VectorResult)) = "+Str$(xpos#)
 Text 10, 30, " Y (Y Vector3(VectorResult)) = "+Str$(ypos#)
 Text 10, 50, " Z (Z Vector3(VectorResult)) = "+Str$(zpos#)
 if upkey()=1 then point object 1,xpos#,ypos#,zpos#:move object 1,2
sync
loop

Function 2DTo3D(Vector, ViewMatrix, ProjectionMatrix, Range#)
   X#=2.0*(Mousex()*1.0)/(Screen Width()*1.0)-1.0
   Y#=1.0-2.0*(MouseY()*1.0)/(Screen Height()*1.0)
   Set Vector3 Vector, X#, Y#, 1.0
   View Matrix4 ViewMatrix
   Projection Matrix4 ProjectionMatrix
   r=Inverse Matrix4(ViewMatrix, ViewMatrix)
   r=Inverse Matrix4(ProjectionMatrix, ProjectionMatrix)
   Transform Coords Vector3 Vector, Vector, ProjectionMatrix
   Normalize Vector3 Vector, Vector
   Multiply Vector3 Vector, Range#
   Transform Coords Vector3 Vector, Vector, ViewMatrix
   Projection Matrix4 ProjectionMatrix
EndFunction


camera:

if mouseclick()=2

    ax#=ax#+(my-mousey())
    ay#=wrapvalue(ay#+(mx-mousex()))
    if ax#>85 then ax#=85
    if ax#<5 then ax#=5
  endif

  mx=mousex()
  my=mousey()

  x#=ox#+r#*sin(ay#)*cos(ax#)
  y#=oy#+r#*sin(ax#)
  z#=oz#+r#*cos(ay#)*cos(ax#)

  position camera x#,y#,z#
  point camera ox#,oy#,oz#


return

`*******************************************************************************
`*                      Mouse selection of matrix tile                         *
`*                           Created: 01/25/2003                               *
`*                     By: DMiTR0S   dmitr0s@newmail.ru                        *
`*******************************************************************************

`*******************************************************************************
`*                            Free License Agreement.                          *
`*                                                                             *
`*   The user/developer agrees to include a visible credit to "Dmitri Kuschev" *
`*   within your programs documentation and or web site.                       *
`*                                                                             *
`*   The user/developer agrees to NOT offer the this source code for Sale.     *
`*                                                                             *
`*   Providing the user/Developer agrees in full to the previous Free          *
`*   License Agreement clauses, the user/developer may use this source         *
`*   code FREELY in all their products, commercial or otherwise.               *
`*                                                                             *
`*******************************************************************************
Set Display Mode 1024, 768, 32

Autocam off
Randomize Timer()

Global Distance# As Float
Global IntersectionPointX# As Float
Global IntersectionPointY# As Float
Global IntersectionPointZ# As Float

MoveVector=1
MatrixClone=1
Matrix=1
Step#=5.0
width=500*2.0
depth=500*2.0
xTiles=30
zTiles=30
r=Make Vector3(MoveVector)

Make Matrix Matrix ,width, depth, xTiles, zTiles
Randomize Matrix Matrix, 50
Update Matrix Matrix
Position Matrix Matrix, RND(10), RND(10), RND(10)

Position Camera width/2, 200.0, -500.0

MakeObjectFromMatrix(Matrix, MatrixClone, xTiles, zTiles, width/(xTiles*1.0), depth/(zTiles*1.0), xTiles+1)

Make Object Sphere 10, 5.0
Color Object 10, RGB(0,255,0)

Make Object Sphere 11, 10.0
Color Object 11, RGB(255,0,0)
Position Object 11, Matrix Position X(Matrix)+RND(width), 0.0, Matrix Position Z(Matrix)+RND(height)

Sync Rate 0
Sync On
Do
   Gosub Camera
   Gosub NewDestinationPoint
   Gosub MoveSphere
   Gosub PrintInfo

   Sync
Loop

NewDestinationPoint:
   GetMatrixCloneIntersection(MatrixClone)
   If GetDistanceToObject()>0.0
      Show Object 10
      Position Object 10, GetIntersectionPointX(), GetIntersectionPointY(), GetIntersectionPointZ()
   EndIf
 Return

 MoveSphere:
   If UpKey()
      Set Vector3 MoveVector, GetIntersectionPointX()-Object Position X(11), GetIntersectionPointY()-Object Position Y(11), GetIntersectionPointZ()-Object Position Z(11)
      If Length Vector3(MoveVector) > Step#
         Normalize Vector3 MoveVector, MoveVector
         Multiply Vector3 MoveVector, Step#
         NewX#=Object Position X(11)+X Vector3(MoveVector)
         NewZ#=Object Position Z(11)+Z Vector3(MoveVector)
         NewY#=Matrix Position Y(Matrix)+Get Ground Height(Matrix, NewX#-Matrix Position X(Matrix), NewZ#-Matrix Position Z(Matrix))
         Position Object 11, NewX#, NewY#, NewZ#
      Else
         Position Object 11, GetIntersectionPointX(), GetIntersectionPointY(), GetIntersectionPointZ()
      EndIf
   EndIf
Return

Camera:
   If MouseClick()=1
      cay#=cay#+(Mousemovex()/10.0)
      cax#=cax#+(Mousemovey()/10.0)
      If cax#<-90.0 Then cax#=-90.0
      If cax#>90.0 Then cax#=90.0
      ax#=Curveangle(cax#,ax#,2.0)
      ay#=Curveangle(cay#,ay#,2.0)
      Rotate Camera Wrapvalue(ax#*-1),Wrapvalue(ay#),Wrapvalue(caz)
   EndIf

   If Mouseclick()=2
      cy#=cy#+MousemoveY()
      Move Camera MousemoveX()
      Position Camera Camera Position X(), cy#, Camera Position Z()
   EndIf
return

PrintInfo:
   Text 10, 10, "Code by: Dmitri Kuschev"
   Text 10, 30, "E-mail: dmitr0s@newmail.ru"
   Text 10, 50, "FPS "+Str$(Screen FPS())
   Text 10, 70, "Press Up to move the red sphere towards the destination point"
   Text 10, 90, "Step: "+Str$(Step#)
   Text 10, 110, "Distance from the camera to the destination point: "+Str$(GetDistanceToObject())
   Text 10, 130, "Position of the red sphere on the matrix in world coordinates:"
   Text 10, 150, "X = "+Str$(Object Position X(11))
   Text 10, 170, "Y = "+Str$(Object Position Y(11))
   Text 10, 190, "Z = "+Str$(Object Position Z(11))
   Text 10, 210, "Position of the red sphere on the matrix in matrix coordinates:"
   Text 10, 230, "X = "+Str$(Object Position X(11)-Matrix Position X(1))
   Text 10, 250, "Y = "+Str$(Object Position Y(11)-Matrix Position Y(1))
   Text 10, 270, "Get Ground Height( Matrix, X, Z ) = "+Str$(Get Ground Height(Matrix, Object Position X(11)-Matrix Position X(1), Object Position Z(11)-Matrix Position Z(1)))
   Text 10, 290, "Z = "+Str$(Object Position Z(11)-Matrix Position Z(1))
   Text 10, 310, "Destination point in world coordinates:"
   Text 10, 330, "GetIntersectionPointX() X = "+Str$(GetIntersectionPointX())
   Text 10, 350, "GetIntersectionPointY() Y = "+Str$(GetIntersectionPointY())
   Text 10, 370, "GetIntersectionPointZ() Z = "+Str$(GetIntersectionPointZ())
   Text 10, 390, "Destination point in matrix coordinates:"
   Text 10, 410, "X = "+Str$(GetIntersectionPointX()-Matrix Position X(1))
   Text 10, 430, "Y = "+Str$(GetIntersectionPointY()-Matrix Position Y(1))
   Text 10, 450, "Z = "+Str$(GetIntersectionPointZ()-Matrix Position Z(1))
Return


Function MakeObjectFromMatrix(Matrix, Object, ,MXSegments, MZSegments, MXSegmentSize#, MZSegmentSize#, VIR)
   VerticlesCount=(MXSegments)*(MZSegments)*6
   FVFFormat=338
   FVFSize=36
   If Memblock Exist(1) Then Delete Memblock 1
   Make Memblock 1, 12+VerticlesCount*FVFSize
   Write Memblock Dword 1, 0, FVFFormat
   Write Memblock Dword 1, 4, FVFSize
   Write Memblock Dword 1, 8, VerticlesCount
   For SX=1 to MXSegments
      For SZ=1 To MZSegments
         Num=MXSegments*(SZ-1)+SX
         Poly1X1#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
         Poly1Y1#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX-1, SZ-1)
         Poly1Z1#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#

         Poly1X2#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
         Poly1Y2#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX-1, SZ)
         Poly1Z2#=Matrix Position Z(Matrix)+SZ*MZSegmentSize#

         Poly1X3#=Matrix Position X(Matrix)+SX*MXSegmentSize#
         Poly1Y3#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX, SZ)
         Poly1Z3#=Matrix Position Z(Matrix)+SZ*MZSegmentSize#

         Poly2X1#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
         Poly2Y1#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX-1, SZ-1)
         Poly2Z1#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#

         Poly2X2#=Matrix Position X(Matrix)+SX*MXSegmentSize#
         Poly2Y2#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX, SZ)
         Poly2Z2#=Matrix Position Z(Matrix)+SZ*MZSegmentSize#

         Poly2X3#=Matrix Position X(Matrix)+SX*MXSegmentSize#
         Poly2Y3#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX, SZ-1)
         Poly2Z3#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#

         Write Memblock Float 1, 12+216*(Num-1), Poly1X1#
         Write Memblock Float 1, 16+216*(Num-1), Poly1Y1#
         Write Memblock Float 1, 20+216*(Num-1), Poly1Z1#
         Write Memblock Float 1, 24+216*(Num-1), 0.0
         Write Memblock Float 1, 28+216*(Num-1), 0.0
         Write Memblock Float 1, 32+216*(Num-1), -1.0
         Write Memblock Dword 1, 36+216*(Num-1), -1
         Write Memblock Float 1, 40+216*(Num-1), 0.0
         Write Memblock Float 1, 44+216*(Num-1), 1.0

         Write Memblock Float 1, 48+216*(Num-1), Poly1X2#
         Write Memblock Float 1, 52+216*(Num-1), Poly1Y2#
         Write Memblock Float 1, 56+216*(Num-1), Poly1Z2#
         Write Memblock Float 1, 60+216*(Num-1), 0.0
         Write Memblock Float 1, 64+216*(Num-1), 0.0
         Write Memblock Float 1, 68+216*(Num-1), -1.0
         Write Memblock Dword 1, 72+216*(Num-1), -1
         Write Memblock Float 1, 76+216*(Num-1), 0.0
         Write Memblock Float 1, 80+216*(Num-1), 0.0

         Write Memblock Float 1, 84+216*(Num-1), Poly1X3#
         Write Memblock Float 1, 88+216*(Num-1), Poly1Y3#
         Write Memblock Float 1, 92+216*(Num-1), Poly1Z3#
         Write Memblock Float 1, 96+216*(Num-1), 0.0
         Write Memblock Float 1, 100+216*(Num-1), 0.0
         Write Memblock Float 1, 104+216*(Num-1), -1.0
         Write Memblock Dword 1, 108+216*(Num-1), -1
         Write Memblock Float 1, 112+216*(Num-1), 1.0
         Write Memblock Float 1, 116+216*(Num-1), 1.0

         Write Memblock Float 1, 120+216*(Num-1), Poly2X1#
         Write Memblock Float 1, 124+216*(Num-1), Poly2Y1#
         Write Memblock Float 1, 128+216*(Num-1), Poly2Z1#
         Write Memblock Float 1, 132+216*(Num-1), 0.0
         Write Memblock Float 1, 136+216*(Num-1), 0.0
         Write Memblock Float 1, 140+216*(Num-1), -1.0
         Write Memblock Dword 1, 144+216*(Num-1), -1
         Write Memblock Float 1, 148+216*(Num-1), 1.0
         Write Memblock Float 1, 152+216*(Num-1), 0.0

         Write Memblock Float 1, 156+216*(Num-1), Poly2X2#
         Write Memblock Float 1, 160+216*(Num-1), Poly2Y2#
         Write Memblock Float 1, 164+216*(Num-1), Poly2Z2#
         Write Memblock Float 1, 168+216*(Num-1), 0.0
         Write Memblock Float 1, 172+216*(Num-1), 0.0
         Write Memblock Float 1, 176+216*(Num-1), -1.0
         Write Memblock Dword 1, 180+216*(Num-1), -1
         Write Memblock Float 1, 184+216*(Num-1), 1.0
         Write Memblock Float 1, 188+216*(Num-1), 1.0

         Write Memblock Float 1, 192+216*(Num-1), Poly2X3#
         Write Memblock Float 1, 196+216*(Num-1), Poly2Y3#
         Write Memblock Float 1, 200+216*(Num-1), Poly2Z3#
         Write Memblock Float 1, 204+216*(Num-1), 0.0
         Write Memblock Float 1, 208+216*(Num-1), 0.0
         Write Memblock Float 1, 212+216*(Num-1), -1.0
         Write Memblock Dword 1, 216+216*(Num-1), -1
         Write Memblock Float 1, 220+216*(Num-1), 0.0
         Write Memblock Float 1, 224+216*(Num-1), 1.0

      Next SZ

   Next SX

  Make Mesh From Memblock 1, 1
  If Object Exist(Object) Then Delete Object Object
  Make Object Object, 1, 22
  Hide Object Object
EndFunction

Function GetMatrixCloneIntersection(Obj)
   DestPoint=2
   ViewMatrix=101
   ProjectionMatrix=102
   r=Make Matrix4(ViewMatrix)
   r=Make Matrix4(ProjectionMatrix)
   r=Make Vector3(DestPoint)
   X#=2.0*(Mousex()*1.0)/(Screen Width()*1.0)-1.0
   Y#=1.0-2.0*(MouseY()*1.0)/(Screen Height()*1.0)
   Set Vector3 DestPoint, X#, Y#, 1.0
   View Matrix4 ViewMatrix
   Projection Matrix4 ProjectionMatrix
   r=Inverse Matrix4(ViewMatrix, ViewMatrix)
   r=Inverse Matrix4(ProjectionMatrix, ProjectionMatrix)
   Transform Coords Vector3 DestPoint, DestPoint, ProjectionMatrix
   Normalize Vector3 DestPoint, DestPoint
   Multiply Vector3 DestPoint, 1000.0
   Transform Coords Vector3 DestPoint, DestPoint, ViewMatrix

   Distance#=Intersect Object(Obj, Camera Position X(), Camera Position Y(), Camera Position Z(), X Vector3(DestPoint), Y Vector3(DestPoint), Z Vector3(DestPoint))

   Set Vector3 DestPoint, (X Vector3(DestPoint)-Camera Position X()), (Y Vector3(DestPoint)-Camera Position Y()), (Z Vector3(DestPoint)-Camera Position Z())
   Normalize Vector3 DestPoint, DestPoint
   Multiply Vector3 DestPoint, Distance#
   Set Vector3 DestPoint, X Vector3(DestPoint)+Camera Position X(), Y Vector3(DestPoint)+Camera Position Y(), Z Vector3(DestPoint)+Camera Position Z()

   If Distance#>0.0
      IntersectionPointX#=X Vector3(DestPoint)
      IntersectionPointY#=Y Vector3(DestPoint)
      IntersectionPointZ#=Z Vector3(DestPoint)
   Else
      IntersectionPointX#=0.0
      IntersectionPointY#=0.0
      IntersectionPointZ#=0.0
   EndIf
EndFunction

Function GetIntersectionPointX()
EndFunction IntersectionPointX#

Function GetIntersectionPointY()
EndFunction IntersectionPointY#

Function GetIntersectionPointZ()
EndFunction IntersectionPointZ#

Function GetDistanceToObject()
EndFunction Distance#

Delete Memblock 1

`*******************************************************************************
`*                      Mouse selection of matrix tile                         *
`*                           Created: 01/25/2003                               *
`*                     By: DMiTR0S   dmitr0s@newmail.ru                        *
`*******************************************************************************

`*******************************************************************************
`*                            Free License Agreement.                          *
`*                                                                             *
`*   The user/developer agrees to include a visible credit to "Dmitri Kuschev" *
`*   within your programs documentation and or web site.                       *
`*                                                                             *
`*   The user/developer agrees to NOT offer the this source code for Sale.     *
`*                                                                             *
`*   Providing the user/Developer agrees in full to the previous Free          *
`*   License Agreement clauses, the user/developer may use this source         *
`*   code FREELY in all their products, commercial or otherwise.               *
`*                                                                             *
`*******************************************************************************
Set Display Mode 1024, 768, 32

Autocam off
Randomize Timer()

Global Distance# As Float
Global IntersectionPointX# As Float
Global IntersectionPointY# As Float
Global IntersectionPointZ# As Float

MoveVector=1
MatrixClone=1
Matrix=1
Step#=5.0
width=500*2.0
depth=500*2.0
xTiles=30
zTiles=30
r=Make Vector3(MoveVector)

Make Matrix Matrix ,width, depth, xTiles, zTiles
Randomize Matrix Matrix, 50
Update Matrix Matrix
Position Matrix Matrix, RND(10), RND(10), RND(10)

Position Camera width/2, 200.0, -500.0

MakeObjectFromMatrix(Matrix, MatrixClone, xTiles, zTiles, width/(xTiles*1.0), depth/(zTiles*1.0), xTiles+1)

Make Object Sphere 10, 5.0
Color Object 10, RGB(0,255,0)

Make Object Sphere 11, 10.0
Color Object 11, RGB(255,0,0)
Position Object 11, Matrix Position X(Matrix)+RND(width), 0.0, Matrix Position Z(Matrix)+RND(height)

Sync Rate 0
Sync On
Do
   Gosub Camera
   Gosub NewDestinationPoint
   Gosub MoveSphere
   Gosub PrintInfo

   Sync
Loop

NewDestinationPoint:
   GetMatrixCloneIntersection(MatrixClone)
   If GetDistanceToObject()>0.0
      Show Object 10
      Position Object 10, GetIntersectionPointX(), GetIntersectionPointY(), GetIntersectionPointZ()
   EndIf
 Return

 MoveSphere:
   If UpKey()
      Set Vector3 MoveVector, GetIntersectionPointX()-Object Position X(11), GetIntersectionPointY()-Object Position Y(11), GetIntersectionPointZ()-Object Position Z(11)
      If Length Vector3(MoveVector) > Step#
         Normalize Vector3 MoveVector, MoveVector
         Multiply Vector3 MoveVector, Step#
         NewX#=Object Position X(11)+X Vector3(MoveVector)
         NewY#=Object Position Y(11)+Y Vector3(MoveVector)
         NewZ#=Object Position Z(11)+Z Vector3(MoveVector)
         Position Object 11, NewX#, NewY#, NewZ#
      Else
         Position Object 11, GetIntersectionPointX(), GetIntersectionPointY(), GetIntersectionPointZ()
      EndIf
   EndIf
Return

Camera:
   If MouseClick()=1
      cay#=cay#+(Mousemovex()/10.0)
      cax#=cax#+(Mousemovey()/10.0)
      If cax#<-90.0 Then cax#=-90.0
      If cax#>90.0 Then cax#=90.0
      ax#=Curveangle(cax#,ax#,2.0)
      ay#=Curveangle(cay#,ay#,2.0)
      Rotate Camera Wrapvalue(ax#*-1),Wrapvalue(ay#),Wrapvalue(caz)
   EndIf

   If Mouseclick()=2
      cy#=cy#+MousemoveY()
      Move Camera MousemoveX()
      Position Camera Camera Position X(), cy#, Camera Position Z()
   EndIf
return

PrintInfo:
   Text 10, 10, "Code by: Dmitri Kuschev"
   Text 10, 30, "E-mail: dmitr0s@newmail.ru"
   Text 10, 50, "FPS "+Str$(Screen FPS())
   Text 10, 70, "Press Up to move the red sphere towards the destination point"
   Text 10, 90, "Step: "+Str$(Step#)
   Text 10, 110, "Distance from the camera to the destination point: "+Str$(GetDistanceToObject())
   Text 10, 130, "Position of the red sphere on the matrix in world coordinates:"
   Text 10, 150, "X = "+Str$(Object Position X(11))
   Text 10, 170, "Y = "+Str$(Object Position Y(11))
   Text 10, 190, "Z = "+Str$(Object Position Z(11))
   Text 10, 210, "Position of the red sphere on the matrix in matrix coordinates:"
   Text 10, 230, "X = "+Str$(Object Position X(11)-Matrix Position X(1))
   Text 10, 250, "Y = "+Str$(Object Position Y(11)-Matrix Position Y(1))
   Text 10, 270, "Get Ground Height( Matrix, X, Z ) = "+Str$(Get Ground Height(Matrix, Object Position X(11)-Matrix Position X(1), Object Position Z(11)-Matrix Position Z(1)))
   Text 10, 290, "Z = "+Str$(Object Position Z(11)-Matrix Position Z(1))
   Text 10, 310, "Destination point in world coordinates:"
   Text 10, 330, "GetIntersectionPointX() X = "+Str$(GetIntersectionPointX())
   Text 10, 350, "GetIntersectionPointY() Y = "+Str$(GetIntersectionPointY())
   Text 10, 370, "GetIntersectionPointZ() Z = "+Str$(GetIntersectionPointZ())
   Text 10, 390, "Destination point in matrix coordinates:"
   Text 10, 410, "X = "+Str$(GetIntersectionPointX()-Matrix Position X(1))
   Text 10, 430, "Y = "+Str$(GetIntersectionPointY()-Matrix Position Y(1))
   Text 10, 450, "Z = "+Str$(GetIntersectionPointZ()-Matrix Position Z(1))
Return


Function MakeObjectFromMatrix(Matrix, Object, ,MXSegments, MZSegments, MXSegmentSize#, MZSegmentSize#, VIR)
   VerticlesCount=(MXSegments)*(MZSegments)*6
   FVFFormat=338
   FVFSize=36
   If Memblock Exist(1) Then Delete Memblock 1
   Make Memblock 1, 12+VerticlesCount*FVFSize
   Write Memblock Dword 1, 0, FVFFormat
   Write Memblock Dword 1, 4, FVFSize
   Write Memblock Dword 1, 8, VerticlesCount
   For SX=1 to MXSegments
      For SZ=1 To MZSegments
         Num=MXSegments*(SZ-1)+SX
         Poly1X1#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
         Poly1Y1#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX-1, SZ-1)
         Poly1Z1#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#

         Poly1X2#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
         Poly1Y2#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX-1, SZ)
         Poly1Z2#=Matrix Position Z(Matrix)+SZ*MZSegmentSize#

         Poly1X3#=Matrix Position X(Matrix)+SX*MXSegmentSize#
         Poly1Y3#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX, SZ)
         Poly1Z3#=Matrix Position Z(Matrix)+SZ*MZSegmentSize#

         Poly2X1#=Matrix Position X(Matrix)+(SX-1)*MXSegmentSize#
         Poly2Y1#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX-1, SZ-1)
         Poly2Z1#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#

         Poly2X2#=Matrix Position X(Matrix)+SX*MXSegmentSize#
         Poly2Y2#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX, SZ)
         Poly2Z2#=Matrix Position Z(Matrix)+SZ*MZSegmentSize#

         Poly2X3#=Matrix Position X(Matrix)+SX*MXSegmentSize#
         Poly2Y3#=Matrix Position Y(Matrix)+Get Matrix Height(1, SX, SZ-1)
         Poly2Z3#=Matrix Position Z(Matrix)+(SZ-1)*MZSegmentSize#

         Write Memblock Float 1, 12+216*(Num-1), Poly1X1#
         Write Memblock Float 1, 16+216*(Num-1), Poly1Y1#
         Write Memblock Float 1, 20+216*(Num-1), Poly1Z1#
         Write Memblock Float 1, 24+216*(Num-1), 0.0
         Write Memblock Float 1, 28+216*(Num-1), 0.0
         Write Memblock Float 1, 32+216*(Num-1), -1.0
         Write Memblock Dword 1, 36+216*(Num-1), -1
         Write Memblock Float 1, 40+216*(Num-1), 0.0
         Write Memblock Float 1, 44+216*(Num-1), 1.0

         Write Memblock Float 1, 48+216*(Num-1), Poly1X2#
         Write Memblock Float 1, 52+216*(Num-1), Poly1Y2#
         Write Memblock Float 1, 56+216*(Num-1), Poly1Z2#
         Write Memblock Float 1, 60+216*(Num-1), 0.0
         Write Memblock Float 1, 64+216*(Num-1), 0.0
         Write Memblock Float 1, 68+216*(Num-1), -1.0
         Write Memblock Dword 1, 72+216*(Num-1), -1
         Write Memblock Float 1, 76+216*(Num-1), 0.0
         Write Memblock Float 1, 80+216*(Num-1), 0.0

         Write Memblock Float 1, 84+216*(Num-1), Poly1X3#
         Write Memblock Float 1, 88+216*(Num-1), Poly1Y3#
         Write Memblock Float 1, 92+216*(Num-1), Poly1Z3#
         Write Memblock Float 1, 96+216*(Num-1), 0.0
         Write Memblock Float 1, 100+216*(Num-1), 0.0
         Write Memblock Float 1, 104+216*(Num-1), -1.0
         Write Memblock Dword 1, 108+216*(Num-1), -1
         Write Memblock Float 1, 112+216*(Num-1), 1.0
         Write Memblock Float 1, 116+216*(Num-1), 1.0

         Write Memblock Float 1, 120+216*(Num-1), Poly2X1#
         Write Memblock Float 1, 124+216*(Num-1), Poly2Y1#
         Write Memblock Float 1, 128+216*(Num-1), Poly2Z1#
         Write Memblock Float 1, 132+216*(Num-1), 0.0
         Write Memblock Float 1, 136+216*(Num-1), 0.0
         Write Memblock Float 1, 140+216*(Num-1), -1.0
         Write Memblock Dword 1, 144+216*(Num-1), -1
         Write Memblock Float 1, 148+216*(Num-1), 1.0
         Write Memblock Float 1, 152+216*(Num-1), 0.0

         Write Memblock Float 1, 156+216*(Num-1), Poly2X2#
         Write Memblock Float 1, 160+216*(Num-1), Poly2Y2#
         Write Memblock Float 1, 164+216*(Num-1), Poly2Z2#
         Write Memblock Float 1, 168+216*(Num-1), 0.0
         Write Memblock Float 1, 172+216*(Num-1), 0.0
         Write Memblock Float 1, 176+216*(Num-1), -1.0
         Write Memblock Dword 1, 180+216*(Num-1), -1
         Write Memblock Float 1, 184+216*(Num-1), 1.0
         Write Memblock Float 1, 188+216*(Num-1), 1.0

         Write Memblock Float 1, 192+216*(Num-1), Poly2X3#
         Write Memblock Float 1, 196+216*(Num-1), Poly2Y3#
         Write Memblock Float 1, 200+216*(Num-1), Poly2Z3#
         Write Memblock Float 1, 204+216*(Num-1), 0.0
         Write Memblock Float 1, 208+216*(Num-1), 0.0
         Write Memblock Float 1, 212+216*(Num-1), -1.0
         Write Memblock Dword 1, 216+216*(Num-1), -1
         Write Memblock Float 1, 220+216*(Num-1), 0.0
         Write Memblock Float 1, 224+216*(Num-1), 1.0

      Next SZ

   Next SX

  Make Mesh From Memblock 1, 1
  If Object Exist(Object) Then Delete Object Object
  Make Object Object, 1, 22
  Hide Object Object
EndFunction

Function GetMatrixCloneIntersection(Obj)
   DestPoint=2
   ViewMatrix=101
   ProjectionMatrix=102
   r=Make Matrix4(ViewMatrix)
   r=Make Matrix4(ProjectionMatrix)
   r=Make Vector3(DestPoint)
   X#=2.0*(Mousex()*1.0)/(Screen Width()*1.0)-1.0
   Y#=1.0-2.0*(MouseY()*1.0)/(Screen Height()*1.0)
   Set Vector3 DestPoint, X#, Y#, 1.0
   View Matrix4 ViewMatrix
   Projection Matrix4 ProjectionMatrix
   r=Inverse Matrix4(ViewMatrix, ViewMatrix)
   r=Inverse Matrix4(ProjectionMatrix, ProjectionMatrix)
   Transform Coords Vector3 DestPoint, DestPoint, ProjectionMatrix
   Normalize Vector3 DestPoint, DestPoint
   Multiply Vector3 DestPoint, 1000.0
   Transform Coords Vector3 DestPoint, DestPoint, ViewMatrix

   Distance#=Intersect Object(Obj, Camera Position X(), Camera Position Y(), Camera Position Z(), X Vector3(DestPoint), Y Vector3(DestPoint), Z Vector3(DestPoint))

   Set Vector3 DestPoint, (X Vector3(DestPoint)-Camera Position X()), (Y Vector3(DestPoint)-Camera Position Y()), (Z Vector3(DestPoint)-Camera Position Z())
   Normalize Vector3 DestPoint, DestPoint
   Multiply Vector3 DestPoint, Distance#
   Set Vector3 DestPoint, X Vector3(DestPoint)+Camera Position X(), Y Vector3(DestPoint)+Camera Position Y(), Z Vector3(DestPoint)+Camera Position Z()

   If Distance#>0.0
      IntersectionPointX#=X Vector3(DestPoint)
      IntersectionPointY#=Y Vector3(DestPoint)
      IntersectionPointZ#=Z Vector3(DestPoint)
   Else
      IntersectionPointX#=0.0
      IntersectionPointY#=0.0
      IntersectionPointZ#=0.0
   EndIf
EndFunction

Function GetIntersectionPointX()
EndFunction IntersectionPointX#

Function GetIntersectionPointY()
EndFunction IntersectionPointY#

Function GetIntersectionPointZ()
EndFunction IntersectionPointZ#

Function GetDistanceToObject()
EndFunction Distance#

Delete Memblock 1

For I=1 to 50
   Make Object Cube I, RND(10)
   Position Object I, 100-RND(200), 100-RND(200), 100-RND(200)
Next I
Text$="Bla Bla"
Sync Rate 0
Sync On
Do
   Inc YRotAngle#, 0.2
   Position Camera NewXValue(0.0, YRotAngle#, 200.0), 10.0, NewZValue(0.0, YRotAngle#, 200.0)
   Point Camera 0.0, 0.0, 0.0
   For I=1 To 50
      X#=Object Screen X(I)
      Y#=Object Screen Y(I)
      Box X#, Y#-Text Height(Text$)*1.5, X#+Text Width(Text$)*1.5, Y#, RGB(255,0,255), RGB(255,0,255), RGB(255,0,255), RGB(255,0,255)
      Center Text X#+(Text Width(Text$)*1.5)/2.0, Y#-(Text Height(Text$)*1.5), Text$
   Next I
   Sync
Loop