rmath3d_plain 1.1.0 → 1.2.4

checksums.yaml

@@ -1,7 +1,7 @@
 ---
-SHA1:
-  metadata.gz: 5b2a5f2fa521f8b303ef746634b8574967b7b703
-  data.tar.gz: ecfe68cac549e516052ea77ddd885051fd36fc81
+SHA256:
+  metadata.gz: 569f5334bd7df6cf7e0f7a857ae589bc8258b0940913807ebaa2dcae60bf711e
+  data.tar.gz: 5f5cf50b25ad23c0a51574f6a4d83ac45cfc820e9c5a337591d66a7574c47145
 SHA512:
-  metadata.gz: 1bda4cabbc91bc183b5881f5c51c0716c3fb8b0f1b3562eb8843e47c3bd3f774c69de6111c4f7412103e81b2258b3476cd7f1b974c6aaa9abbc8b05c3d447750
-  data.tar.gz: a58f95113e0d9e3bef1488fdf26d79c2069f3fb0654eeeaec0bd6741ceb43941ec2af3e8cd7c016d3524d5adcda99933da38e1c4496cfb7f1c1335aea23ef245
+  metadata.gz: 13a0b47c640f875de55331962edc75ee9f48f43554e196b98073c297b7ebf6f02e47449c0280f4f9a377a281d0c2607a154a37838cdbaf2d90aa6536a7429c72
+  data.tar.gz: da55b4c5ed4322637162396c67ff45d2f2c870ec7242592309d9bc31dce68e9b3014449bc4f781cde6754cc02e1ce5f5aacbcb98f57cbe867d031dc37a07338d

data/ChangeLog

@@ -1,3 +1,29 @@
+2020-07-23  vaiorabbit  <http://twitter.com/vaiorabbit>
+
+	* rmath3d.c, rmath3d_plain.rb (RMtx4): Added lookAtLH, perspectiveLH, etc.
+
+2020-06-21  vaiorabbit  <http://twitter.com/vaiorabbit>
+
+	* rmath3d.c, rmath3d_plain.rb (RMtx4): Added argument 'ndc_homogeneous' for projection matrix APIs.
+	Set true for the environment with Z coordinate ranges from -1 to +1 (OpenGL), and false otherwise (Direct3D, Metal)
+	https://www.slideshare.net/Mark_Kilgard/opengl-32-and-more/26-Direct3D_vs_OpenGL_Coordinate_System
+	https://metashapes.com/blog/opengl-metal-projection-matrix-problem/
+	* rmath3d_plain.rb: Removed Fixnum and Bignum symbols (deprecated and unified into Integer class since Ruby 2.4)
+
+2017-07-22  vaiorabbit  <http://twitter.com/vaiorabbit>
+
+	* Added 'Integer' type for argument type branching (constant ::Fixnum is deprecated since ruby 2.4).
+
+2015-08-23  vaiorabbit  <http://twitter.com/vaiorabbit>
+
+	* rmath3d.c (def ==): Removed control path returning Qnil.
+	* rmath3d_plain.rb (def ==): Removed control path returning Qnil.
+
+2015-05-02  vaiorabbit  <http://twitter.com/vaiorabbit>
+
+	* Ruby 1.9.3 and prior versions are no longer supported. Ref.: https://www.ruby-lang.org/en/news/2015/02/23/support-for-ruby-1-9-3-has-ended/
+	* rmath3d.c: Data_Get_Struct -> TypedData_Get_Struct, etc.
+
 2015-04-12  vaiorabbit  <http://twitter.com/vaiorabbit>
 
 	* RVec2.c|h, RMtx2.c|h: Added.

data/LICENSE.txt

@@ -1,5 +1,5 @@
 rmath3d : Ruby math module for 3D Applications
-Copyright (c) 2008-2015 vaiorabbit  <http://twitter.com/vaiorabbit>
+Copyright (c) 2008-2020 vaiorabbit  <http://twitter.com/vaiorabbit>
 
 This software is provided 'as-is', without any express or implied
 warranty. In no event will the authors be held liable for any damages

data/README.md

@@ -4,9 +4,12 @@
 
 rmath3d is a math module for 3D game programming and computer graphics.
 
-*   Last Update: Apr 12, 2015
+*   Last Update: Jul 23, 2020
 *   Since: Jul 20, 2008
 
+* rmath3d (C Extension Library Implementation) [![Gem Version](https://badge.fury.io/rb/rmath3d.svg)](https://badge.fury.io/rb/rmath3d) [![Gem](https://img.shields.io/gem/dt/rmath3d.svg)](https://rubygems.org/gems/rmath3d)
+* rmath3d_plain (Pure Ruby Implementation) [![Gem Version](https://badge.fury.io/rb/rmath3d_plain.svg)](https://badge.fury.io/rb/rmath3d_plain) [![Gem](https://img.shields.io/gem/dt/rmath3d_plain.svg)](https://rubygems.org/gems/rmath3d_plain)
+
 ## Features ##
 
 ### Supports frequently-used vector and matrix classes ###
@@ -30,11 +33,8 @@ rmath3d is a math module for 3D game programming and computer graphics.
 Notice: This library provides native extension. You must setup develop environment (or DevKit) before installation.
 
 *   Ruby
-    *   ruby 2.0.0p247 (2013-06-27 revision 41674) [x86_64-darwin12.4.0]
-    *   ruby 2.0.0p247 (2013-06-27) [i386-mingw32] With Development Kit installed.
-        *   I used: DevKit-mingw64-32-4.7.2-20130224-1151-sfx.exe
-        *   Unpack the archive -> "> ruby dk.rb init" -> edit config.yml (just add your ruby foldier) -> "> ruby dk.rb install" 
-        *   Ref.: http://blog.mattwynne.net/2010/10/12/installing-ruby-gems-with-native-extensions-on-windows/
+    *   ruby 2.7.1p83 (2020-03-31 revision a0c7c23c9c) [x64-mingw32]
+    *   Ruby 2.3 and prior versions are no longer supported.
 
 ## Building rmath3d.{so|bundle} ##
 

data/lib/rmath3d/rmath3d_plain.rb

@@ -29,7 +29,7 @@ module RMath3D
                0.0, 0.0 ]
       when 1
         case a[0]
-        when Fixnum, Float
+        when Float, Integer
           @e = [ a[0], a[0],
                  a[0], a[0] ]
         when RMtx2
@@ -46,7 +46,7 @@ module RMath3D
           for col in 0...2 do
             index = 2*row + col
             case a[index]
-            when Fixnum, Float
+            when Float, Integer
               setElement( row, col, a[index] )
             else
               raise TypeError, "RMtx2#initialize : Unknown type #{a[0].class}."
@@ -88,7 +88,7 @@ module RMath3D
     #
     def coerce
       case arg
-      when Fixnum, Float, Bignum
+      when Float, Integer
         return [ self, arg ]
       else
         raise TypeError, "RMtx2#coerce : #{arg.self} can't be coerced into  #{self.class}."
@@ -411,7 +411,7 @@ module RMath3D
     #
     def *( arg )
       case arg
-      when Fixnum, Float, Bignum
+      when Float, Integer
         return RMtx2.new( arg*self.e00, arg*self.e01,
                           arg*self.e10, arg*self.e11 )
 
@@ -440,19 +440,18 @@ module RMath3D
     # mtx1 == mtx2 : evaluates equality.
     #
     def ==( other )
-      if ( other.class != RMtx2 )
-        raise TypeError, "RMtx2#==(other) : Unknown type #{other.class} given as RMtx2."
-        return nil
-      end
-
-      for row in 0...2 do
-        for col in 0...2 do
-          if ( (getElement(row,col) - other.getElement(row,col)).abs > TOLERANCE )
-            return false
+      if other.class == RMtx2
+        for row in 0...2 do
+          for col in 0...2 do
+            if ( (getElement(row,col) - other.getElement(row,col)).abs > TOLERANCE )
+              return false
+            end
           end
         end
+        return true
+      else
+        return false
       end
-      return true
     end
 
     #
@@ -504,7 +503,7 @@ module RMath3D
     #
     def mul!( other )
       case other
-      when Fixnum, Float, Bignum
+      when Float, Integer
         self.e00 = other*self.e00
         self.e01 = other*self.e01
         self.e10 = other*self.e10
@@ -561,7 +560,7 @@ module RMath3D
                0.0, 0.0, 0.0 ]
       when 1
         case a[0]
-        when Fixnum, Float
+        when Float, Integer
           @e = [ a[0], a[0], a[0],
                  a[0], a[0], a[0],
                  a[0], a[0], a[0] ]
@@ -580,7 +579,7 @@ module RMath3D
           for col in 0...3 do
             index = 3*row + col
             case a[index]
-            when Fixnum, Float
+            when Float, Integer
               setElement( row, col, a[index] )
             else
               raise TypeError, "RMtx3#initialize : Unknown type #{a[0].class}."
@@ -623,7 +622,7 @@ module RMath3D
     #
     def coerce
       case arg
-      when Fixnum, Float, Bignum
+      when Float, Integer
         return [ self, arg ]
       else
         raise TypeError, "RMtx3#coerce : #{arg.self} can't be coerced into  #{self.class}."
@@ -1099,7 +1098,7 @@ module RMath3D
     #
     def *( arg )
       case arg
-      when Fixnum, Float, Bignum
+      when Float, Integer
         return RMtx3.new( arg*self.e00, arg*self.e01, arg*self.e02,
                           arg*self.e10, arg*self.e11, arg*self.e12,
                           arg*self.e20, arg*self.e21, arg*self.e22 )
@@ -1129,19 +1128,18 @@ module RMath3D
     # mtx1 == mtx2 : evaluates equality.
     #
     def ==( other )
-      if ( other.class != RMtx3 )
-        raise TypeError, "RMtx3#==(other) : Unknown type #{other.class} given as RMtx3."
-        return nil
-      end
-
-      for row in 0...3 do
-        for col in 0...3 do
-          if ( (getElement(row,col) - other.getElement(row,col)).abs > TOLERANCE )
-            return false
+      if other.class == RMtx3
+        for row in 0...3 do
+          for col in 0...3 do
+            if ( (getElement(row,col) - other.getElement(row,col)).abs > TOLERANCE )
+              return false
+            end
           end
         end
+        return true
+      else
+        return false
       end
-      return true
     end
 
     #
@@ -1193,7 +1191,7 @@ module RMath3D
     #
     def mul!( other )
       case other
-      when Fixnum, Float, Bignum
+      when Float, Integer
         self.e00 = other*self.e00
         self.e01 = other*self.e01
         self.e02 = other*self.e02
@@ -1261,7 +1259,7 @@ module RMath3D
                0.0, 0.0, 0.0, 0.0 ]
       when 1
         case a[0]
-        when Fixnum, Float
+        when Float, Integer
           @e = [ a[0], a[0], a[0], a[0],
                  a[0], a[0], a[0], a[0],
                  a[0], a[0], a[0], a[0],
@@ -1282,7 +1280,7 @@ module RMath3D
           for col in 0...4 do
             index = 4*row + col
             case a[index]
-            when Fixnum, Float
+            when Float, Integer
               setElement( row, col, a[index] )
             else
               raise TypeError, "RMtx4#initialize : Unknown type #{a[0].class}."
@@ -1326,7 +1324,7 @@ module RMath3D
     #
     def coerce
       case arg
-      when Fixnum, Float, Bignum
+      when Float, Integer
         return [ self, arg ]
       else
         raise TypeError, "RMtx4#coerce : #{arg.self} can't be coerced into  #{self.class}."
@@ -1839,6 +1837,39 @@ module RMath3D
       return self
     end
 
+    #
+    # call-seq: lookAtLH(eye,at,up) -> self
+    #
+    # Builds a viewing matrix for a left-handed coordinate system from:
+    # * eye position (+eye+: RVec3)
+    # * a point looking at (+at+: RVec3)
+    # * up vector (+up+: RVec3)
+    #
+    def lookAtLH( eye, at, up )
+      setIdentity()
+
+      axis_z = (at - eye).normalize!
+      axis_x = RVec3.cross( up, axis_z ).normalize!
+      axis_y = RVec3.cross( axis_z, axis_x )
+
+      self.e00 = axis_x[0]
+      self.e01 = axis_x[1]
+      self.e02 = axis_x[2]
+      self.e03 = -RVec3.dot( axis_x, eye )
+
+      self.e10 = axis_y[0]
+      self.e11 = axis_y[1]
+      self.e12 = axis_y[2]
+      self.e13 = -RVec3.dot( axis_y, eye )
+
+      self.e20 = axis_z[0]
+      self.e21 = axis_z[1]
+      self.e22 = axis_z[2]
+      self.e23 = -RVec3.dot( axis_z, eye )
+
+      return self
+    end
+
     #
     # call-seq: lookAtRH(eye,at,up) -> self
     #
@@ -1873,43 +1904,146 @@ module RMath3D
     end
 
     #
-    # call-seq: perspectiveRH(width,height,znear,zfar) -> self
+    # call-seq: perspectiveLH(width,height,znear,zfar,ndc_convention) -> self
     #
     # Builds a perspective projection matrix for a right-handed coordinate system from:
     # * View volume width (+width+)
     # * View volume height (+height+)
     # * Near clip plane distance (+znear+)
     # * Far clip plane distance (+zfar+)
+    # * Set true for the environment with Z coordinate ranges from -1 to +1 (OpenGL), and false otherwise (Direct3D, Metal) (+ndc_homogeneous+)
     #
-    def perspectiveRH( width, height, znear, zfar )
-      perspectiveOffCenterRH(-width/2.0, width/2.0, -height/2.0, height/2.0, znear, zfar )
-      return self
+    def perspectiveLH( width, height, znear, zfar, ndc_homogeneous = true)
+      return perspectiveOffCenterLH(-width/2.0, width/2.0, -height/2.0, height/2.0, znear, zfar, ndc_homogeneous )
     end
 
     #
-    # call-seq: perspectiveFovRH(fovy,aspect,znear,zfar) -> self
+    # call-seq: perspectiveFovLH(fovy,aspect,znear,zfar,ndc_homogeneous) -> self
     #
     # Builds a perspective projection matrix for a right-handed coordinate system from:
     # * Field of view in y direction (+fovy+ radian)
     # * Aspect ratio (+aspect+)
     # * Near clip plane distance (+znear+)
     # * Far clip plane distance (+zfar+)
+    # * Set true for the environment with Z coordinate ranges from -1 to +1 (OpenGL), and false otherwise (Direct3D, Metal) (+ndc_homogeneous+)
+    #
+    def perspectiveFovLH( fovy_radian, aspect, znear, zfar, ndc_homogeneous = true)
+      # Ref.: https://www.scratchapixel.com/lessons/3d-basic-rendering/perspective-and-orthographic-projection-matrix/opengl-perspective-projection-matrix
+      top = Math::tan(fovy_radian / 2.0) * znear
+      bottom = -top
+      right = top * aspect
+      left = -right
+      return perspectiveOffCenterLH(left, right, bottom, top, znear, zfar, ndc_homogeneous)
+    end
+
+    #
+    # call-seq: perspectiveOffCenterLH(left,right,bottom,top,znear,zfar) -> self
+    #
+    # Builds a perspective projection matrix for a right-handed coordinate system from:
+    # * Minimum value of the view volume width (+left+)
+    # * Maximum value of the view volume width (+right+)
+    # * Minimum value of the view volume height (+bottom+)
+    # * Maximum value of the view volume height (+top+)
+    # * Near clip plane distance (+znear+)
+    # * Far clip plane distance (+zfar+)
+    # * Set true for the environment with Z coordinate ranges from -1 to +1 (OpenGL), and false otherwise (Direct3D, Metal) (+ndc_homogeneous+)
+    #
+    def perspectiveOffCenterLH( left, right, bottom, top, znear, zfar, ndc_homogeneous = true)
+      a = (right+left) / (right-left)
+      b = (top+bottom) / (top-bottom)
+      c = ndc_homogeneous ? -(zfar+znear) / (zfar-znear) : -zfar / (zfar-znear)
+      d = ndc_homogeneous ? -(2*znear*zfar) / (zfar-znear) : -(znear*zfar) / (zfar-znear)
+
+      setZero()
+
+      setElement( 0, 0, 2*znear/(right-left) )
+      setElement( 0, 2, -a )
+      setElement( 1, 1, 2*znear/(top-bottom) )
+      setElement( 1, 2, -b )
+      setElement( 2, 2, -c )
+      setElement( 2, 3, d )
+      setElement( 3, 2, 1.0 )
+
+      return self
+    end
+
+    #
+    # call-seq: orthoLH(width,height,znear,zfar) -> self
+    #
+    # Builds a orthogonal projection matrix for a right-handed coordinate system from:
+    # * View volume width (+width+)
+    # * View volume height (+height+)
+    # * Near clip plane distance (+znear+)
+    # * Far clip plane distance (+zfar+)
+    # * Set true for the environment with Z coordinate ranges from -1 to +1 (OpenGL), and false otherwise (Direct3D, Metal) (+ndc_homogeneous+)
     #
-    def perspectiveFovRH( fovy_radian, aspect, znear, zfar )
-      f = Math::tan( fovy_radian / 2.0 )
-      f = 1.0 / f
+    def orthoLH( width, height, znear, zfar, ndc_homogeneous = true)
+      orthoOffCenterLH( -width/2.0, width/2.0, -height/2.0, height/2.0, znear, zfar, ndc_homogeneous )
+      return self
+    end
+
+    #
+    # call-seq: orthoOffCenterLH(left,right,bottom,top,znear,zfar) -> self
+    #
+    # Builds a orthogonal projection matrix for a right-handed coordinate system from:
+    # * Minimum value of the view volume width (+left+)
+    # * Maximum value of the view volume width (+right+)
+    # * Minimum value of the view volume height (+bottom+)
+    # * Maximum value of the view volume height (+top+)
+    # * Near clip plane distance (+znear+)
+    # * Far clip plane distance (+zfar+)
+    # * Set true for the environment with Z coordinate ranges from -1 to +1 (OpenGL), and false otherwise (Direct3D, Metal) (+ndc_homogeneous+)
+    #
+    def orthoOffCenterLH( left, right, bottom, top, znear, zfar, ndc_homogeneous = true)
+      tx = -(right+left) / (right-left)
+      ty = -(top+bottom) / (top-bottom)
+      tz = ndc_homogeneous ? -(zfar+znear) / (zfar-znear) : -znear / (zfar-znear)
 
       setIdentity()
-      setElement( 0, 0, f / aspect )
-      setElement( 1, 1, f )
-      setElement( 2, 2, (zfar+znear)/(znear-zfar) )
-      setElement( 2, 3, 2*zfar*znear/(znear-zfar) )
-      setElement( 3, 2, -1.0 )
-      setElement( 3, 3, 0.0 )
+
+      setElement( 0, 0, 2.0/(right-left) )
+      setElement( 0, 3, tx )
+      setElement( 1, 1, 2.0/(top-bottom) )
+      setElement( 1, 3, ty )
+      setElement( 2, 2, (ndc_homogeneous ? 2.0 : 1.0)/(zfar-znear) )
+      setElement( 2, 3, tz )
 
       return self
     end
 
+    #
+    # call-seq: perspectiveRH(width,height,znear,zfar,ndc_convention) -> self
+    #
+    # Builds a perspective projection matrix for a right-handed coordinate system from:
+    # * View volume width (+width+)
+    # * View volume height (+height+)
+    # * Near clip plane distance (+znear+)
+    # * Far clip plane distance (+zfar+)
+    # * Set true for the environment with Z coordinate ranges from -1 to +1 (OpenGL), and false otherwise (Direct3D, Metal) (+ndc_homogeneous+)
+    #
+    def perspectiveRH( width, height, znear, zfar, ndc_homogeneous = true)
+      return perspectiveOffCenterRH(-width/2.0, width/2.0, -height/2.0, height/2.0, znear, zfar, ndc_homogeneous )
+    end
+
+    #
+    # call-seq: perspectiveFovRH(fovy,aspect,znear,zfar,ndc_homogeneous) -> self
+    #
+    # Builds a perspective projection matrix for a right-handed coordinate system from:
+    # * Field of view in y direction (+fovy+ radian)
+    # * Aspect ratio (+aspect+)
+    # * Near clip plane distance (+znear+)
+    # * Far clip plane distance (+zfar+)
+    # * Set true for the environment with Z coordinate ranges from -1 to +1 (OpenGL), and false otherwise (Direct3D, Metal) (+ndc_homogeneous+)
+    #
+    def perspectiveFovRH( fovy_radian, aspect, znear, zfar, ndc_homogeneous = true)
+      # Ref.: https://www.scratchapixel.com/lessons/3d-basic-rendering/perspective-and-orthographic-projection-matrix/opengl-perspective-projection-matrix
+      top = Math::tan(fovy_radian / 2.0) * znear
+      bottom = -top
+      right = top * aspect
+      left = -right
+      return perspectiveOffCenterRH(left, right, bottom, top, znear, zfar, ndc_homogeneous)
+    end
+
     #
     # call-seq: perspectiveOffCenterRH(left,right,bottom,top,znear,zfar) -> self
     #
@@ -1920,14 +2054,15 @@ module RMath3D
     # * Maximum value of the view volume height (+top+)
     # * Near clip plane distance (+znear+)
     # * Far clip plane distance (+zfar+)
+    # * Set true for the environment with Z coordinate ranges from -1 to +1 (OpenGL), and false otherwise (Direct3D, Metal) (+ndc_homogeneous+)
     #
-    def perspectiveOffCenterRH( left, right, bottom, top, znear, zfar )
+    def perspectiveOffCenterRH( left, right, bottom, top, znear, zfar, ndc_homogeneous = true)
       a = (right+left) / (right-left)
       b = (top+bottom) / (top-bottom)
-      c = -(zfar+znear) / (zfar-znear)
-      d = -(2*znear*zfar) / (zfar-znear)
+      c = ndc_homogeneous ? -(zfar+znear) / (zfar-znear) : -zfar / (zfar-znear)
+      d = ndc_homogeneous ? -(2*znear*zfar) / (zfar-znear) : -(znear*zfar) / (zfar-znear)
 
-      setIdentity()
+      setZero()
 
       setElement( 0, 0, 2*znear/(right-left) )
       setElement( 0, 2, a )
@@ -1936,7 +2071,6 @@ module RMath3D
       setElement( 2, 2, c )
       setElement( 2, 3, d )
       setElement( 3, 2, -1.0 )
-      setElement( 3, 3,  0.0 )
 
       return self
     end
@@ -1949,9 +2083,10 @@ module RMath3D
     # * View volume height (+height+)
     # * Near clip plane distance (+znear+)
     # * Far clip plane distance (+zfar+)
+    # * Set true for the environment with Z coordinate ranges from -1 to +1 (OpenGL), and false otherwise (Direct3D, Metal) (+ndc_homogeneous+)
     #
-    def orthoRH( width, height, znear, zfar )
-      orthoOffCenterRH( -width/2.0, width/2.0, -height/2.0, height/2.0, znear, zfar )
+    def orthoRH( width, height, znear, zfar, ndc_homogeneous = true)
+      orthoOffCenterRH( -width/2.0, width/2.0, -height/2.0, height/2.0, znear, zfar, ndc_homogeneous )
       return self
     end
 
@@ -1965,11 +2100,12 @@ module RMath3D
     # * Maximum value of the view volume height (+top+)
     # * Near clip plane distance (+znear+)
     # * Far clip plane distance (+zfar+)
+    # * Set true for the environment with Z coordinate ranges from -1 to +1 (OpenGL), and false otherwise (Direct3D, Metal) (+ndc_homogeneous+)
     #
-    def orthoOffCenterRH( left, right, bottom, top, znear, zfar )
-      tx = (right+left) / (right-left)
-      ty = (top+bottom) / (top-bottom)
-      tz = (zfar+znear) / (zfar-znear)
+    def orthoOffCenterRH( left, right, bottom, top, znear, zfar, ndc_homogeneous = true)
+      tx = -(right+left) / (right-left)
+      ty = -(top+bottom) / (top-bottom)
+      tz = ndc_homogeneous ? -(zfar+znear) / (zfar-znear) : -znear / (zfar-znear)
 
       setIdentity()
 
@@ -1977,7 +2113,7 @@ module RMath3D
       setElement( 0, 3, tx )
       setElement( 1, 1, 2.0/(top-bottom) )
       setElement( 1, 3, ty )
-      setElement( 2, 2, -2.0/(zfar-znear) )
+      setElement( 2, 2, -(ndc_homogeneous ? 2.0 : 1.0)/(zfar-znear) )
       setElement( 2, 3, tz )
 
       return self
@@ -2050,7 +2186,7 @@ module RMath3D
     #
     def *( arg )
       case arg
-      when Fixnum, Float, Bignum
+      when Float, Integer
         return RMtx4.new( arg*self.e00, arg*self.e01, arg*self.e02, arg*self.e03,
                           arg*self.e10, arg*self.e11, arg*self.e12, arg*self.e13,
                           arg*self.e20, arg*self.e21, arg*self.e22, arg*self.e23,
@@ -2081,19 +2217,18 @@ module RMath3D
     # mtx1 == mtx2 : evaluates equality.
     #
     def ==( other )
-      if ( other.class != RMtx4 )
-        raise TypeError, "RMtx4#==(other) : Unknown type #{other.class} given as RMtx4."
-        return nil
-      end
-
-      for row in 0...4 do
-        for col in 0...4 do
-          if ( (getElement(row,col) - other.getElement(row,col)).abs > TOLERANCE )
-            return false
+      if other.class == RMtx4
+        for row in 0...4 do
+          for col in 0...4 do
+            if ( (getElement(row,col) - other.getElement(row,col)).abs > TOLERANCE )
+              return false
+            end
           end
         end
+        return true
+      else
+        return false
       end
-      return true
     end
 
     #
@@ -2145,7 +2280,7 @@ module RMath3D
     #
     def mul!( other )
       case other
-      when Fixnum, Float, Bignum
+      when Float, Integer
         self.e00 = other*self.e00
         self.e01 = other*self.e01
         self.e02 = other*self.e02
@@ -2227,7 +2362,7 @@ module RMath3D
         @e = [0.0, 0.0, 0.0, 0.0]
       when 1
         case a[0]
-        when Fixnum, Float
+        when Float, Integer
           @e = [ a[0], a[0], a[0], a[0] ]
         when RQuat
           @e = [ a[0].x, a[0].y, a[0].z, a[0].w ]
@@ -2238,7 +2373,7 @@ module RMath3D
       when 4
         a.each_with_index do |elem, index|
           case elem
-          when Fixnum, Float
+          when Float, Integer
             @e[index] = elem
           else
             raise TypeError, "RQuat#initialize : Unknown type #{elem.class}."
@@ -2258,7 +2393,7 @@ module RMath3D
     # Returns human-readable string.
     #
     def to_s
-      return "( #{@e[0]}, #{@e[1]}, #{@e[2]}, #{@e[3]} )\n"
+      return "( #{@e[0]}, #{@e[1]}, #{@e[2]}, #{@e[3]} )"
     end
 
     #
@@ -2277,7 +2412,7 @@ module RMath3D
     #
     def coerce( arg )
       case arg
-      when Fixnum, Float, Bignum
+      when Float, Integer
         return [ self, arg ]
       else
         raise TypeError, "RQuat#coerce : #{arg.self} can't be coerced into  #{self.class}."
@@ -2616,7 +2751,7 @@ module RMath3D
         z = q1w*q2z + q1x*q2y - q1y*q2x + q1z*q2w
         w = q1w*q2w - q1x*q2x - q1y*q2y - q1z*q2z
         return RQuat.new( x, y, z, w )
-      when Fixnum, Float
+      when Float, Integer
         return RQuat.new( @e[0]*arg, @e[1]*arg, @e[2]*arg, @e[3]*arg )
       else
         raise TypeError, "RQuat#* : Unknown type #{arg}."
@@ -2630,16 +2765,15 @@ module RMath3D
     # quat1 == quat2 : evaluates equality.
     #
     def ==( other )
-      if other.class != RQuat
-        raise TypeError, "RQuat#== : Unknown type #{other.class}."
-        return nil
-      end
-
-      if  (x-other.x).abs<=Float::EPSILON &&
-          (y-other.y).abs<=Float::EPSILON &&
-          (z-other.z).abs<=Float::EPSILON &&
-          (w-other.w).abs<=Float::EPSILON
-        return true
+      if other.class == RQuat
+        if  (x-other.x).abs<=Float::EPSILON &&
+            (y-other.y).abs<=Float::EPSILON &&
+            (z-other.z).abs<=Float::EPSILON &&
+            (w-other.w).abs<=Float::EPSILON
+          return true
+        else
+          return false
+        end
       else
         return false
       end
@@ -2712,7 +2846,7 @@ module RMath3D
 
         return self
 
-      when Fixnum, Float
+      when Float, Integer
         self.x *= other
         self.y *= other
         self.z *= other
@@ -2828,7 +2962,7 @@ module RMath3D
         @e = [0.0, 0.0]
       when 1
         case a[0]
-        when Fixnum, Float
+        when Float, Integer
           @e = [ a[0], a[0] ]
         when RVec2
           @e = [ a[0].x, a[0].y ]
@@ -2839,7 +2973,7 @@ module RMath3D
       when 2
         a.each_with_index do |elem, index|
           case elem
-          when Fixnum, Float
+          when Float, Integer
             @e[index] = elem
           else
             raise TypeError, "RVec2#initialize : Unknown type #{elem.class}."
@@ -2859,7 +2993,7 @@ module RMath3D
     # Returns human-readable string.
     #
     def to_s
-      return "( #{@e[0]}, #{@e[1]} )\n"
+      return "( #{@e[0]}, #{@e[1]} )"
     end
 
     #
@@ -2878,7 +3012,7 @@ module RMath3D
     #
     def coerce( arg )
       case arg
-      when Fixnum, Float, Bignum
+      when Float, Integer
         return [ self, arg ]
       else
         raise TypeError, "RVec2#coerce : #{arg.self} can't be coerced into  #{self.class}."
@@ -3067,7 +3201,7 @@ module RMath3D
     #
     def *( arg )
       case arg
-      when Fixnum, Float
+      when Float, Integer
         return RVec2.new( @e[0]*arg, @e[1]*arg )
       else
         raise TypeError, "RVec2#* : Unknown type #{arg}."
@@ -3081,14 +3215,13 @@ module RMath3D
     # vec1 == vec2 : evaluates equality.
     #
     def ==( other )
-      if other.class != RVec2
-        raise TypeError, "RVec2#== : Unknown type #{other.class}."
-        return nil
-      end
-
-      if  (x-other.x).abs<=Float::EPSILON &&
-          (y-other.y).abs<=Float::EPSILON
-        return true
+      if other.class == RVec2
+        if  (x-other.x).abs<=Float::EPSILON &&
+            (y-other.y).abs<=Float::EPSILON
+          return true
+        else
+          return false
+        end
       else
         return false
       end
@@ -3134,7 +3267,7 @@ module RMath3D
     # vec1 *= vec2
     #
     def mul!( arg )
-      if arg.class != Fixnum && arg.class != Float
+      if !(arg.class == Float || arg.class == Integer)
         raise TypeError, "RVec2#mul! : Unknown type #{arg.class}."
         return nil
       end
@@ -3168,7 +3301,7 @@ module RMath3D
         @e = [0.0, 0.0, 0.0]
       when 1
         case a[0]
-        when Fixnum, Float
+        when Float, Integer
           @e = [ a[0], a[0], a[0] ]
         when RVec3
           @e = [ a[0].x, a[0].y, a[0].z ]
@@ -3179,7 +3312,7 @@ module RMath3D
       when 3
         a.each_with_index do |elem, index|
           case elem
-          when Fixnum, Float
+          when Float, Integer
             @e[index] = elem
           else
             raise TypeError, "RVec3#initialize : Unknown type #{elem.class}."
@@ -3199,7 +3332,7 @@ module RMath3D
     # Returns human-readable string.
     #
     def to_s
-      return "( #{@e[0]}, #{@e[1]}, #{@e[2]} )\n"
+      return "( #{@e[0]}, #{@e[1]}, #{@e[2]} )"
     end
 
     #
@@ -3218,7 +3351,7 @@ module RMath3D
     #
     def coerce( arg )
       case arg
-      when Fixnum, Float, Bignum
+      when Float, Integer
         return [ self, arg ]
       else
         raise TypeError, "RVec3#coerce : #{arg.self} can't be coerced into  #{self.class}."
@@ -3522,9 +3655,9 @@ module RMath3D
       t_z = q.w*self[2] + q.x*self[1] - q.y*self[0]
       t_w =             - q.x*self[0] - q.y*self[1] - q.z*self[2]
 
-      result.x = -t_w*q.x + t_x*q.w - t_y*q.z + t_z*q.y;
-      result.y = -t_w*q.y + t_x*q.z + t_y*q.w - t_z*q.x;
-      result.z = -t_w*q.z - t_x*q.y + t_y*q.x + t_z*q.w;
+      result.x = -t_w*q.x + t_x*q.w - t_y*q.z + t_z*q.y
+      result.y = -t_w*q.y + t_x*q.z + t_y*q.w - t_z*q.x
+      result.z = -t_w*q.z - t_x*q.y + t_y*q.x + t_z*q.w
 
       return result
     end
@@ -3538,9 +3671,9 @@ module RMath3D
       t_z = q.w*self[2] + q.x*self[1] - q.y*self[0]
       t_w =             - q.x*self[0] - q.y*self[1] - q.z*self[2]
 
-      self.x = -t_w*q.x + t_x*q.w - t_y*q.z + t_z*q.y;
-      self.y = -t_w*q.y + t_x*q.z + t_y*q.w - t_z*q.x;
-      self.z = -t_w*q.z - t_x*q.y + t_y*q.x + t_z*q.w;
+      self.x = -t_w*q.x + t_x*q.w - t_y*q.z + t_z*q.y
+      self.y = -t_w*q.y + t_x*q.z + t_y*q.w - t_z*q.x
+      self.z = -t_w*q.z - t_x*q.y + t_y*q.x + t_z*q.w
 
       return self
     end
@@ -3621,7 +3754,7 @@ module RMath3D
     #
     def *( arg )
       case arg
-      when Fixnum, Float
+      when Float, Integer
         return RVec3.new( @e[0]*arg, @e[1]*arg, @e[2]*arg )
       else
         raise TypeError, "RVec3#* : Unknown type #{arg}."
@@ -3635,15 +3768,14 @@ module RMath3D
     # vec1 == vec2 : evaluates equality.
     #
     def ==( other )
-      if other.class != RVec3
-        raise TypeError, "RVec3#== : Unknown type #{other.class}."
-        return nil
-      end
-
-      if  (x-other.x).abs<=Float::EPSILON &&
-          (y-other.y).abs<=Float::EPSILON &&
-          (z-other.z).abs<=Float::EPSILON
-        return true
+      if other.class == RVec3
+        if  (x-other.x).abs<=Float::EPSILON &&
+            (y-other.y).abs<=Float::EPSILON &&
+            (z-other.z).abs<=Float::EPSILON
+          return true
+        else
+          return false
+        end
       else
         return false
       end
@@ -3691,7 +3823,7 @@ module RMath3D
     # vec1 *= vec2
     #
     def mul!( arg )
-      if arg.class != Fixnum && arg.class != Float
+      if !(arg.class == Float || arg.class == Integer)
         raise TypeError, "RVec3#mul! : Unknown type #{arg.class}."
         return nil
       end
@@ -3726,7 +3858,7 @@ module RMath3D
         @e = [0.0, 0.0, 0.0, 0.0]
       when 1
         case a[0]
-        when Fixnum, Float
+        when Float, Integer
           @e = [ a[0], a[0], a[0], a[0] ]
         when RVec3
           @e = [ a[0].x, a[0].y, a[0].z, 0.0 ]
@@ -3739,7 +3871,7 @@ module RMath3D
       when 4
         a.each_with_index do |elem, index|
           case elem
-          when Fixnum, Float
+          when Float, Integer
             @e[index] = elem
           else
             raise TypeError, "RVec4#initialize : Unknown type #{elem.class}."
@@ -3759,7 +3891,7 @@ module RMath3D
     # Returns human-readable string.
     #
     def to_s
-      return "( #{@e[0]}, #{@e[1]}, #{@e[2]}, #{@e[3]} )\n"
+      return "( #{@e[0]}, #{@e[1]}, #{@e[2]}, #{@e[3]} )"
     end
 
     #
@@ -3778,7 +3910,7 @@ module RMath3D
     #
     def coerce( arg )
       case arg
-      when Fixnum, Float, Bignum
+      when Float, Integer
         return [ self, arg ]
       else
         raise TypeError, "RVec4#coerce : #{arg.self} can't be coerced into  #{self.class}."
@@ -4068,7 +4200,7 @@ module RMath3D
     #
     def *( arg )
       case arg
-      when Fixnum, Float
+      when Float, Integer
         return RVec4.new( @e[0]*arg, @e[1]*arg, @e[2]*arg, @e[3]*arg )
       else
         raise TypeError, "RVec4#* : Unknown type #{arg}."
@@ -4082,16 +4214,15 @@ module RMath3D
     # vec1 == vec2 : evaluates equality.
     #
     def ==( other )
-      if other.class != RVec4
-        raise TypeError, "RVec4#== : Unknown type #{other.class}."
-        return nil
-      end
-
-      if  (x-other.x).abs<=Float::EPSILON &&
-          (y-other.y).abs<=Float::EPSILON &&
-          (z-other.z).abs<=Float::EPSILON &&
-          (w-other.w).abs<=Float::EPSILON
-        return true
+      if other.class == RVec4
+        if  (x-other.x).abs<=Float::EPSILON &&
+            (y-other.y).abs<=Float::EPSILON &&
+            (z-other.z).abs<=Float::EPSILON &&
+            (w-other.w).abs<=Float::EPSILON
+          return true
+        else
+          return false
+        end
       else
         return false
       end
@@ -4141,7 +4272,7 @@ module RMath3D
     # vec1 *= vec2
     #
     def mul!( other )
-      if other.class != Fixnum && other.class != Float
+      if !(other.class == Float || other.class == Integer)
         raise TypeError, "RVec4#mul! : Unknown type #{other.class}."
         return nil
       end
@@ -4159,7 +4290,7 @@ end
 
 =begin
 RMath : Ruby math module for 3D Applications
-Copyright (c) 2008- vaiorabbit  <http://twitter.com/vaiorabbit>
+Copyright (c) 2008-2020 vaiorabbit  <http://twitter.com/vaiorabbit>
 
 This software is provided 'as-is', without any express or implied
 warranty. In no event will the authors be held liable for any damages

data/sample/opengl-bindings/load_matrix.rb

@@ -104,7 +104,7 @@ class App
   def size_callback( window_handle, w, h )
     glViewport( 0, 0, w, h )
     glMatrixMode( GL_PROJECTION )
-    @mtxProj.perspectiveFovRH( 30.0*Math::PI/180.0, w.to_f/h.to_f, 0.1, 1000.0 )
+    @mtxProj.perspectiveFovRH( 30.0*Math::PI/180.0, w.to_f/h.to_f, 0.1, 1000.0, true )
     glLoadMatrixf( @mtxProj.to_a.pack('F16') )
 
     @window_width  = w
@@ -119,7 +119,7 @@ class App
     @at  = RVec3.new(0.0,  0.0,  0.0)
     @up  = RVec3.new(0.0,  1.0,  0.0)
     @mtxLookAt = RMtx4.new.lookAtRH( @eye, @at, @up )
-    @mtxProj = RMtx4.new.perspectiveFovRH( 30.0*Math::PI/180.0, @window_width.to_f/@window_height.to_f, 0.1, 1000.0 )
+    @mtxProj = RMtx4.new.perspectiveFovRH( 30.0*Math::PI/180.0, @window_width.to_f/@window_height.to_f, 0.1, 1000.0, true )
 
     @light_pos = [2.5,0,5,1]
     @light_diffuse  = [1,1,1,1]

data/sample/opengl2/load_matrix.rb

@@ -42,7 +42,7 @@ class App
   def reshape( width, height )
     glViewport( 0, 0, width, height )
     glMatrixMode( GL_PROJECTION )
-    @mtxProj.perspectiveFovRH( 30.0*Math::PI/180.0, width.to_f/height.to_f, 0.1, 1000.0 )
+    @mtxProj.perspectiveFovRH( 30.0*Math::PI/180.0, width.to_f/height.to_f, 0.1, 1000.0, true )
     glLoadMatrix( @mtxProj )
 
     @window_width  = width
@@ -74,7 +74,7 @@ class App
     @at  = RVec3.new(0.0,  0.0,  0.0)
     @up  = RVec3.new(0.0,  1.0,  0.0)
     @mtxLookAt = RMtx4.new.lookAtRH( @eye, @at, @up )
-    @mtxProj = RMtx4.new.perspectiveFovRH( 30.0*Math::PI/180.0, @window_width.to_f/@window_height.to_f, 0.1, 1000.0 )
+    @mtxProj = RMtx4.new.perspectiveFovRH( 30.0*Math::PI/180.0, @window_width.to_f/@window_height.to_f, 0.1, 1000.0, true )
 
     @light_pos = [2.5,0,5,1]
     @light_diffuse  = [1,1,1,1]

data/test/test.rb

@@ -1,4 +1,4 @@
-begin 
+begin
   require 'rmath3d/rmath3d'
 rescue LoadError
   require 'rmath3d/rmath3d_plain'

data/test/test_RMtx4.rb

@@ -485,6 +485,29 @@ class TC_RMtx4 < Minitest::Test
     end
   end
 
+  def test_lookAtLH
+    pEye = RVec3.new( 10, 10, 10 )
+    vDir = ( RVec3.new(0,0,0) - pEye ).normalize! # staring at (0,0,0)
+    vUp = RVec3.new( 0, 1, 0 )
+    vRight = RVec3.cross( vUp, vDir ).normalize!
+    vUp = RVec3.cross( vDir, vRight ).normalize!
+
+    m0 = RMtx4.new( vRight.x, vRight.y, vRight.z, -RVec3.dot(pEye,vRight),
+                    vUp.x,    vUp.y,    vUp.z,    -RVec3.dot(pEye,vUp),
+                    vDir.x,   vDir.y,   vDir.z,   -RVec3.dot(pEye,vDir),
+                    0.0,      0.0,      0.0,       1.0    )
+
+    m1 = RMtx4.new.lookAtLH( RVec3.new(10,10,10), # posistion
+                             RVec3.new(0,0,0),    # at
+                             RVec3.new(0,1,0) )   # up
+
+    for r in 0...4 do
+      for c in 0...4 do
+        assert_in_delta( m0.getElement(r,c), m1.getElement(r,c), @tolerance )
+      end
+    end
+  end
+
   def test_lookAtRH
     pEye = RVec3.new( 10, 10, 10 )
     vDir = ( pEye - RVec3.new(0,0,0) ).normalize! # staring at (0,0,0)
@@ -526,7 +549,7 @@ class TC_RMtx4 < Minitest::Test
                     0.0,         2*z_n/height, 0.0,                  0.0,
                     0.0,         0.0,         -(z_f+z_n)/(z_f-z_n), -2.0*z_f*z_n / (z_f-z_n),
                     0.0,         0.0,         -1.0,                  0.0 )
-    m1 = RMtx4.new.perspectiveRH( width, height, z_n, z_f )
+    m1 = RMtx4.new.perspectiveRH( width, height, z_n, z_f, true )
 
     for r in 0...4 do
       for c in 0...4 do
@@ -540,15 +563,15 @@ class TC_RMtx4 < Minitest::Test
     fovy = 2.0 * Math::atan( (height/2.0) / z_n )
     f = 1.0/Math::tan( fovy/2.0 )
 
-    m2 = RMtx4.new( f/aspect, 0.0,  0.0,                 0.0,
-                    0.0,      f,    0.0,                 0.0,
-                    0.0,      0.0,  (z_f+z_n)/(z_n-z_f), 2*z_f*z_n/(z_n-z_f),
-                    0.0,      0.0, -1.0,                 0.0 )
-    m3 = RMtx4.new.perspectiveFovRH( fovy, aspect, z_n, z_f );
+    m2 = RMtx4.new( f/aspect, 0.0,  0.0,                   0.0,
+                    0.0,      f,    0.0,                   0.0,
+                    0.0,      0.0,  -(z_f+z_n)/(z_f-z_n), -2*z_f*z_n/(z_f-z_n),
+                    0.0,      0.0,  -1.0,                  0.0 )
+    m3 = RMtx4.new.perspectiveFovRH( fovy, aspect, z_n, z_f, true );
 
     for r in 0...4 do
       for c in 0...4 do
-        assert_in_delta( m2.getElement(r,c), m2.getElement(r,c), @tolerance )
+        assert_in_delta( m2.getElement(r,c), m3.getElement(r,c), @tolerance )
       end
     end
 
@@ -563,7 +586,68 @@ class TC_RMtx4 < Minitest::Test
                     0.0,                0.0,                c,   d,
                     0.0,                0.0,               -1.0, 0.0 )
 
-    m5 = RMtx4.new.perspectiveOffCenterRH( left, right, bottom, top, z_n, z_f )
+    m5 = RMtx4.new.perspectiveOffCenterRH( left, right, bottom, top, z_n, z_f, true )
+
+    for r in 0...4 do
+      for c in 0...4 do
+        assert_in_delta( m4.getElement(r,c), m5.getElement(r,c), @tolerance )
+      end
+    end
+  end
+
+  def test_perspectiveLH
+    left  = -640.0
+    right =  640.0
+    bottom = -360.0
+    top    =  360.0
+    z_n = 1.0
+    z_f = 1000.0
+    width = right - left
+    height = top - bottom
+    aspect = width/height
+
+    # RMtx4#perspectiveLH
+    m0 = RMtx4.new( 2*z_n/width, 0.0,          0.0,                  0.0,
+                    0.0,         2*z_n/height, 0.0,                  0.0,
+                    0.0,         0.0,          (z_f+z_n)/(z_f-z_n), -2.0*z_f*z_n / (z_f-z_n),
+                    0.0,         0.0,          1.0,                  0.0 )
+    m1 = RMtx4.new.perspectiveLH( width, height, z_n, z_f, true )
+
+    for r in 0...4 do
+      for c in 0...4 do
+        assert_in_delta( m0.getElement(r,c), m1.getElement(r,c), @tolerance )
+      end
+    end
+
+    # RMtx4#perspectiveFovLH
+
+    # [NOTE] tan(fovy/2) == (height/2)/z_n
+    fovy = 2.0 * Math::atan( (height/2.0) / z_n )
+    f = 1.0/Math::tan( fovy/2.0 )
+
+    m2 = RMtx4.new( f/aspect, 0.0,  0.0,                 0.0,
+                    0.0,      f,    0.0,                 0.0,
+                    0.0,      0.0,  (z_f+z_n)/(z_f-z_n), -2.0*z_f*z_n/(z_f-z_n),
+                    0.0,      0.0,  1.0,                 0.0 )
+    m3 = RMtx4.new.perspectiveFovLH( fovy, aspect, z_n, z_f, true );
+    for r in 0...4 do
+      for c in 0...4 do
+        assert_in_delta( m2.getElement(r,c), m3.getElement(r,c), @tolerance )
+      end
+    end
+
+    # RMtx4#perspectiveOffCenterLH
+
+    a = (right+left)/(right-left)
+    b = (top+bottom)/(top-bottom)
+    c = -(z_f+z_n)/(z_f-z_n)
+    d = -2.0*z_f*z_n/(z_f-z_n)
+    m4 = RMtx4.new( 2*z_n/(right-left), 0.0,                -a,   0.0,
+                    0.0,                2*z_n/(top-bottom), -b,   0.0,
+                    0.0,                0.0,                -c,   d,
+                    0.0,                0.0,                 1.0, 0.0 )
+
+    m5 = RMtx4.new.perspectiveOffCenterLH( left, right, bottom, top, z_n, z_f, true )
 
     for r in 0...4 do
       for c in 0...4 do
@@ -584,14 +668,14 @@ class TC_RMtx4 < Minitest::Test
     height = top - bottom
 
     # RMtx4#orthoRH
-    tx = (right+left)/width
-    ty = (top+bottom)/height
-    tz = (z_f+z_n)/(z_f-z_n)
+    tx = -(right+left)/width
+    ty = -(top+bottom)/height
+    tz = -(z_f+z_n)/(z_f-z_n)
     m0 = RMtx4.new( 2.0/width, 0.0,         0.0,           tx,
                     0.0,       2.0/height,  0.0,           ty,
                     0.0,       0.0,        -2.0/(z_f-z_n), tz,
                     0.0,       0.0,         0.0,           1.0 )
-    m1 = RMtx4.new.orthoRH( width, height, z_n, z_f )
+    m1 = RMtx4.new.orthoRH( width, height, z_n, z_f, true )
 
     for r in 0...4 do
       for c in 0...4 do
@@ -600,14 +684,57 @@ class TC_RMtx4 < Minitest::Test
     end
 
     # RMtx4#orthoOffCenterRH
-    tx = (right+left)/(right-left)
-    ty = (top+bottom)/(top-bottom)
-    tz = (z_f+z_n)/(z_f-z_n)
+    tx = -(right+left)/(right-left)
+    ty = -(top+bottom)/(top-bottom)
+    tz = -(z_f+z_n)/(z_f-z_n)
     m2 = RMtx4.new( 2.0/(right-left), 0.0,               0.0,           tx,
                     0.0,              2.0/(top-bottom),  0.0,           ty,
                     0.0,              0.0,              -2.0/(z_f-z_n), tz,
                     0.0,              0.0,               0.0,           1.0 )
-    m3 = RMtx4.new.orthoOffCenterRH( left, right, bottom, top, z_n, z_f )
+    m3 = RMtx4.new.orthoOffCenterRH( left, right, bottom, top, z_n, z_f, true )
+
+    for r in 0...4 do
+      for c in 0...4 do
+        assert_in_delta( m2.getElement(r,c), m3.getElement(r,c), @tolerance )
+      end
+    end
+  end
+
+  def test_orthoLH
+    left  = -640.0
+    right =  640.0
+    bottom = -360.0
+    top    =  360.0
+    z_n = 1.0
+    z_f = 1000.0
+    width = right - left
+    height = top - bottom
+
+    # RMtx4#orthoLH
+    tx = -(right+left)/width
+    ty = -(top+bottom)/height
+    tz = -(z_f+z_n)/(z_f-z_n)
+    m0 = RMtx4.new( 2.0/width, 0.0,         0.0,           tx,
+                    0.0,       2.0/height,  0.0,           ty,
+                    0.0,       0.0,         2.0/(z_f-z_n), tz,
+                    0.0,       0.0,         0.0,           1.0 )
+    m1 = RMtx4.new.orthoLH( width, height, z_n, z_f, true )
+
+    for r in 0...4 do
+      for c in 0...4 do
+        assert_in_delta( m0.getElement(r,c), m1.getElement(r,c), @tolerance )
+      end
+    end
+
+    # RMtx4#orthoOffCenterLH
+    tx = -(right+left)/(right-left)
+    ty = -(top+bottom)/(top-bottom)
+    tz = -(z_f+z_n)/(z_f-z_n)
+    m2 = RMtx4.new( 2.0/(right-left), 0.0,               0.0,           tx,
+                    0.0,              2.0/(top-bottom),  0.0,           ty,
+                    0.0,              0.0,               2.0/(z_f-z_n), tz,
+                    0.0,              0.0,               0.0,           1.0 )
+    m3 = RMtx4.new.orthoOffCenterLH( left, right, bottom, top, z_n, z_f, true )
 
     for r in 0...4 do
       for c in 0...4 do

metadata

@@ -1,17 +1,19 @@
 --- !ruby/object:Gem::Specification
 name: rmath3d_plain
 version: !ruby/object:Gem::Version
-  version: 1.1.0
+  version: 1.2.4
 platform: ruby
 authors:
 - vaiorabbit
-autorequire: 
+autorequire:
 bindir: bin
 cert_chain: []
-date: 2015-04-12 00:00:00.000000000 Z
+date: 2020-07-23 00:00:00.000000000 Z
 dependencies: []
-description: |
-  Provides vector2/3/4, matrix2x2/3x3/4x4 and quaternion in plain Ruby form.
+description: 'Provides vector2/3/4, matrix2x2/3x3/4x4 and quaternion in plain Ruby
+  form.
+
+  '
 email:
 - vaiorabbit@gmail.com
 executables: []
@@ -35,9 +37,9 @@ files:
 - test/test_RVec4.rb
 homepage: https://github.com/vaiorabbit/rmath3d
 licenses:
-- zlib/libpng
+- Zlib
 metadata: {}
-post_install_message: 
+post_install_message:
 rdoc_options: []
 require_paths:
 - lib
@@ -45,16 +47,15 @@ required_ruby_version: !ruby/object:Gem::Requirement
   requirements:
   - - ">="
     - !ruby/object:Gem::Version
-      version: '0'
+      version: 2.4.0
 required_rubygems_version: !ruby/object:Gem::Requirement
   requirements:
   - - ">="
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubyforge_project: 
-rubygems_version: 2.4.6
-signing_key: 
+rubygems_version: 3.1.2
+signing_key:
 specification_version: 4
 summary: Ruby Math Module for 3D Applications
 test_files: []