Merge pull request #412 from recp/sse_only

separate SSE and SSE2

.github/workflows/ci.yml

@@ -0,0 +1,236 @@
+name: CI
+
+on:
+  push:
+    branches: [ "master" ]
+  pull_request:
+    branches: [ "master" ]
+
+jobs:
+  build_autotools:
+    name: Autotools / ${{ matrix.os }}
+    runs-on: ${{ matrix.os }}
+    strategy:
+      fail-fast: false
+      matrix:
+        os: [macos-12, macos-14, ubuntu-22.04]
+
+    steps:
+    - uses: actions/checkout@v4
+
+    - name: Install Autotools
+      if: runner.os == 'macOS'
+      run: brew upgrade && brew install autoconf automake libtool
+
+    - name: Generate Autotools
+      run: ./autogen.sh
+
+    - name: Configure Autotools
+      run: ./configure
+
+    - name: Build
+      run: make
+
+    - name: Test
+      run: make check
+
+  build_cmake_ios:
+    name: CMake / iOS
+    runs-on: macos-14
+
+    steps:
+    - uses: actions/checkout@v4
+
+    - name: Configure CMake
+      run: |
+        cmake \
+          -B build \
+          -GXcode \
+          -DCMAKE_SYSTEM_NAME=iOS \
+          -DCMAKE_BUILD_TYPE=Release \
+          -DCMAKE_XCODE_ATTRIBUTE_CODE_SIGNING_ALLOWED=NO \
+          -DCGLM_STATIC=ON \
+          -DCGLM_USE_TEST=ON
+
+    - name: Build
+      run: cmake --build build
+
+  build_cmake_macos:
+    name: CMake / ${{ matrix.os }}
+    runs-on: ${{ matrix.os }}
+    strategy:
+      fail-fast: false
+      matrix:
+        os: [macos-12, macos-14]
+
+    steps:
+    - uses: actions/checkout@v4
+
+    - name: Install Ninja
+      if: runner.os == 'macOS'
+      run: brew upgrade && brew install ninja
+
+    - name: Configure CMake
+      run: |
+        cmake \
+          -B build \
+          -GNinja \
+          -DCMAKE_BUILD_TYPE=Release \
+          -DCGLM_STATIC=ON \
+          -DCGLM_USE_TEST=ON
+
+    - name: Build
+      run: cmake --build build
+
+    - name: Test
+      working-directory: build
+      run: ./tests
+
+  build_cmake_ubuntu:
+    name: CMake / ${{ matrix.target.os }} / ${{ matrix.target.cc }}
+    runs-on: ${{ matrix.target.os }}
+    strategy:
+      fail-fast: false
+      matrix:
+        target:
+          - { os: ubuntu-20.04, cc: gcc-11 }
+          - { os: ubuntu-22.04, cc: gcc-12 }
+          - { os: ubuntu-22.04, cc: gcc-13 }
+          - { os: ubuntu-20.04, cc: clang-12 }
+          - { os: ubuntu-22.04, cc: clang-15 }
+    steps:
+    - uses: actions/checkout@v4
+
+    - name: Install Compiler and Ninja
+      run: |
+        sudo apt-get update -y
+        sudo apt-get install -y ${{ matrix.target.cc }} ninja-build
+
+    - name: Configure CMake
+      run: |
+        cmake \
+          -B build \
+          -GNinja \
+          -DCMAKE_C_COMPILER=${{ matrix.target.cc }} \
+          -DCMAKE_BUILD_TYPE=Release \
+          -DCGLM_STATIC=ON \
+          -DCGLM_USE_TEST=ON
+
+    - name: Build
+      run: cmake --build build
+
+    - name: Test
+      working-directory: build
+      run: ./tests
+
+  build_cmake_windows:
+    name: CMake / ${{ matrix.platform.name }}
+    runs-on: windows-2022
+    strategy:
+      fail-fast: false
+      matrix:
+        platform:
+        - { name: Windows (x64),          flags: -A x64 }
+        - { name: Windows (x86),          flags: -A Win32 }
+        - { name: Windows (clang-cl x64), flags: -T ClangCL -A x64 }
+        - { name: Windows (clang-cl x86), flags: -T ClangCL -A Win32 }
+        - { name: Windows (ARM),          flags: -A ARM, skip_tests: true, skip_build: true } # This fails to build.
+        - { name: Windows (ARM64),        flags: -A ARM64, skip_tests: true }
+        - { name: UWP (ARM64),            flags: -A ARM64, -DCMAKE_SYSTEM_NAME=WindowsStore -DCMAKE_SYSTEM_VERSION="10.0", skip_tests: true }
+        - { name: UWP (x64),              flags: -A x64 -DCMAKE_SYSTEM_NAME=WindowsStore -DCMAKE_SYSTEM_VERSION="10.0", skip_tests: true }
+
+    steps:
+    - uses: actions/checkout@v4
+
+    - name: Configure CMake
+      run: cmake -B build `
+          -DCGLM_STATIC=ON `
+          -DCGLM_USE_TEST=ON `
+          ${{ matrix.platform.flags }}
+
+    - name: Build
+      if: ${{ !matrix.platform.skip_build }}
+      run: cmake --build build --config Release --parallel
+
+    - name: Test
+      if: ${{ !matrix.platform.skip_tests }}
+      working-directory: build
+      run: .\Release\tests.exe
+
+  build_documentation:
+    name: Documentation
+    runs-on: ubuntu-22.04
+
+    steps:
+    - uses: actions/checkout@v4
+
+    - uses: actions/setup-python@v5
+      with:
+        python-version: '3.12'
+
+    - name: Install Dependencies
+      working-directory: docs
+      run: python3 -m pip install -r requirements.txt
+
+    - name: Build
+      working-directory: docs
+      run: sphinx-build -W --keep-going source build
+
+  build_meson:
+    name: Meson / ${{ matrix.os }}
+    runs-on: ${{ matrix.os }}
+    strategy:
+      fail-fast: false
+      matrix:
+        os: [macos-14, ubuntu-22.04]
+
+    steps:
+    - uses: actions/checkout@v4
+
+    - uses: actions/setup-python@v5
+      with:
+        python-version: '3.12'
+        cache: 'pip'
+
+    - name: Install meson
+      run: python3 -m pip install meson ninja
+
+    - name: Build
+      run: meson setup build -Dbuildtype=release --default-library=static -Dbuild_tests=true
+
+    - name: Test
+      run: meson test -C build
+
+  build_msbuild:
+    name: MSBuild / Windows
+    runs-on: windows-2022
+
+    # This has no test yet.
+    # It could also try building for ARM, ARM64, ARM64EC, but those fail currently.
+    steps:
+    - uses: actions/checkout@v4
+
+    - uses: microsoft/setup-msbuild@v2
+
+    - name: Build (x86)
+      working-directory: win
+      run: msbuild cglm.vcxproj /p:Configuration=Release /p:Platform=x86 /p:BuildInParallel=true
+
+    - name: Build (x64)
+      working-directory: win
+      run: msbuild cglm.vcxproj /p:Configuration=Release /p:Platform=x64 /p:BuildInParallel=true
+
+  build_swift:
+    name: Swift ${{ matrix.swift }} / ${{ matrix.os }}
+    runs-on: ${{ matrix.os }}
+    strategy:
+      fail-fast: false
+      matrix:
+        os: [macos-12, macos-14, ubuntu-22.04]
+
+    # This has no test yet.
+    steps:
+      - uses: actions/checkout@v4
+
+      - name: Build
+        run: swift build

.github/workflows/cmake-wasm.yml

@@ -0,0 +1,97 @@
+name: CMake WebAssembly
+
+on:
+  push:
+    branches: [ "master" ]
+  pull_request:
+    branches: [ "master" ]
+
+env:
+  wasmtime_version: v7.0.0
+  wasmer_version: v3.1.1
+
+jobs:
+  build_wasi_sdk:
+    strategy:
+      matrix:
+        BUILD_TYPE: [Release, Debug, RelWithDebInfo, MinSizeRel]
+        C_FLAGS: ['', '-msimd128']
+        wasi_sdk_version: [19, 20]
+    # The CMake configure and build commands are platform agnostic and should work equally well on Windows or Mac.
+    # You can convert this to a matrix build if you need cross-platform coverage.
+    # See: https://docs.github.com/en/free-pro-team@latest/actions/learn-github-actions/managing-complex-workflows#using-a-build-matrix
+    runs-on: ubuntu-latest
+
+    steps:
+    - uses: actions/checkout@v4
+
+    - name: Downloading wasi-sdk
+      run: |
+        cd ${{github.workspace}}
+        wget --no-verbose https://github.com/WebAssembly/wasi-sdk/releases/download/wasi-sdk-${{matrix.wasi_sdk_version}}/wasi-sdk-${{matrix.wasi_sdk_version}}.0-linux.tar.gz
+        tar xf wasi-sdk-${{matrix.wasi_sdk_version}}.0-linux.tar.gz
+
+    - name: Configure CMake
+      # Configure CMake in a 'build' subdirectory. `CMAKE_BUILD_TYPE` is only required if you are using a single-configuration generator such as make.
+      # See https://cmake.org/cmake/help/latest/variable/CMAKE_BUILD_TYPE.html?highlight=cmake_build_type
+      run: cmake -B ${{github.workspace}}/build -DCMAKE_BUILD_TYPE=${{matrix.BUILD_TYPE}} -DCMAKE_C_FLAGS="${{matrix.C_FLAGS}}" -DCMAKE_TOOLCHAIN_FILE=${{github.workspace}}/wasi-sdk-${{matrix.wasi_sdk_version}}.0/share/cmake/wasi-sdk.cmake -DWASI_SDK_PREFIX=${{github.workspace}}/wasi-sdk-${{matrix.wasi_sdk_version}}.0 -DCGLM_USE_TEST=ON
+
+    - name: Build
+      # Build your program with the given configuration
+      run: cmake --build ${{github.workspace}}/build --config ${{matrix.BUILD_TYPE}}
+
+    - name: Test with wasmtime
+      run: |
+        cd ${{github.workspace}}
+        ls -lh ${{github.workspace}}/build/
+        wget --no-verbose https://github.com/bytecodealliance/wasmtime/releases/download/${{env.wasmtime_version}}/wasmtime-${{env.wasmtime_version}}-x86_64-linux.tar.xz
+        tar xf wasmtime-${{env.wasmtime_version}}-x86_64-linux.tar.xz
+        ./wasmtime-${{env.wasmtime_version}}-x86_64-linux/wasmtime run --wasm-features simd ${{github.workspace}}/build/tests
+
+    - name: Test with wasmer
+      run: |
+        cd ${{github.workspace}}
+        mkdir wasmer
+        cd wasmer
+        wget --no-verbose https://github.com/wasmerio/wasmer/releases/download/${{env.wasmer_version}}/wasmer-linux-amd64.tar.gz
+        tar xf wasmer-linux-amd64.tar.gz
+        ./bin/wasmer run --enable-simd ${{github.workspace}}/build/tests
+
+  build_emsdk:
+    strategy:
+      matrix:
+        BUILD_TYPE: [Release, Debug, RelWithDebInfo, MinSizeRel]
+        C_FLAGS: ['', '-msimd128', '-msse -msse2 -msimd128', '-msse -msse2 -msse3 -msse4 -msimd128']
+    runs-on: ubuntu-latest
+    steps:
+    - uses: actions/checkout@v4
+    - name: Setup emsdk
+      uses: mymindstorm/setup-emsdk@v12
+
+    - name: Verify emsdk
+      run: emcc -v
+
+    - name: Configure CMake
+      # Configure CMake in a 'build' subdirectory. `CMAKE_BUILD_TYPE` is only required if you are using a single-configuration generator such as make.
+      # See https://cmake.org/cmake/help/latest/variable/CMAKE_BUILD_TYPE.html?highlight=cmake_build_type
+      run: emcmake cmake -B ${{github.workspace}}/build -DCMAKE_BUILD_TYPE=${{matrix.BUILD_TYPE}} -DCMAKE_C_FLAGS="${{matrix.C_FLAGS}}" -DCMAKE_EXE_LINKER_FLAGS="-s STANDALONE_WASM" -DCGLM_STATIC=ON -DCGLM_USE_TEST=ON
+
+    - name: Build
+      # Build your program with the given configuration
+      run: cmake --build ${{github.workspace}}/build --config ${{matrix.BUILD_TYPE}}
+
+    - name: Test with wasmtime
+      run: |
+        cd ${{github.workspace}}
+        ls -lh ${{github.workspace}}/build/
+        wget --no-verbose https://github.com/bytecodealliance/wasmtime/releases/download/${{env.wasmtime_version}}/wasmtime-${{env.wasmtime_version}}-x86_64-linux.tar.xz
+        tar xf wasmtime-${{env.wasmtime_version}}-x86_64-linux.tar.xz
+        ./wasmtime-${{env.wasmtime_version}}-x86_64-linux/wasmtime run --wasm-features simd ${{github.workspace}}/build/tests.wasm
+    - name: Test with wasmer
+      run: |
+        cd ${{github.workspace}}
+        mkdir wasmer
+        cd wasmer
+        wget --no-verbose https://github.com/wasmerio/wasmer/releases/download/${{env.wasmer_version}}/wasmer-linux-amd64.tar.gz
+        tar xf wasmer-linux-amd64.tar.gz
+        ./bin/wasmer run --enable-simd ${{github.workspace}}/build/tests.wasm

.github/workflows/meson-wasm.yml

@@ -0,0 +1,79 @@
+name: Meson WebAssembly
+
+on:
+  push:
+    branches: [ "master" ]
+  pull_request:
+    branches: [ "master" ]
+
+env:
+  wasmtime_version: v7.0.0
+  wasmer_version: v3.1.1
+
+jobs:
+  build_emsdk:
+    strategy:
+      matrix:
+        BUILD_TYPE: [debug, debugoptimized, release, minsize]
+        C_FLAGS: ['', '-msimd128', '-msse -msse2 -msimd128', '-msse -msse2 -msse3 -msse4 -msimd128']
+    runs-on: ubuntu-latest
+    steps:
+    - uses: actions/checkout@v4
+    - name: Setup emsdk
+      uses: mymindstorm/setup-emsdk@v13
+
+    - name: Verify emsdk
+      run: emcc -v
+
+    - name: Creating cross file
+      run: |
+        cat << EOF > ${{github.workspace}}/meson_cross_emsdk.txt
+        [binaries]
+        c = '`which emcc`'
+        cpp = '`which em++`'
+        ar = '`which emar`'
+        
+        [built-in options]
+        
+        c_args = ['-Wno-unused-parameter']
+        c_link_args = ['-s', 'STANDALONE_WASM']
+        cpp_args = ['-Wno-unused-parameter']
+        cpp_link_args = ['-s', 'STANDALONE_WASM']
+        
+        [host_machine]
+        
+        system = 'emscripten'
+        cpu_family = 'wasm32'
+        cpu = 'wasm32'
+        endian = 'little'
+
+        EOF
+        cat ${{github.workspace}}/meson_cross_emsdk.txt
+
+    - uses: actions/setup-python@v4
+    
+    - name: Install meson
+      run: |
+        sudo python3 -m pip install meson ninja
+
+    - name: Build with meson
+      run: |
+        meson setup build -Dbuildtype=${{matrix.BUILD_TYPE}} --cross-file ${{github.workspace}}/meson_cross_emsdk.txt --default-library=static -Dbuild_tests=true
+        meson test -C build
+
+    - name: Test with wasmtime
+      run: |
+        cd ${{github.workspace}}
+        ls -lh ${{github.workspace}}/build/
+        wget --no-verbose https://github.com/bytecodealliance/wasmtime/releases/download/${{env.wasmtime_version}}/wasmtime-${{env.wasmtime_version}}-x86_64-linux.tar.xz
+        tar xf wasmtime-${{env.wasmtime_version}}-x86_64-linux.tar.xz
+        ./wasmtime-${{env.wasmtime_version}}-x86_64-linux/wasmtime run --wasm-features simd ${{github.workspace}}/build/tests.wasm
+
+    - name: Test with wasmer
+      run: |
+        cd ${{github.workspace}}
+        mkdir wasmer
+        cd wasmer
+        wget --no-verbose https://github.com/wasmerio/wasmer/releases/download/${{env.wasmer_version}}/wasmer-linux-amd64.tar.gz
+        tar xf wasmer-linux-amd64.tar.gz
+        ./bin/wasmer run --enable-simd ${{github.workspace}}/build/tests.wasm

.readthedocs.yaml

@@ -0,0 +1,39 @@
+# Read the Docs configuration file for Sphinx projects
+# See https://docs.readthedocs.io/en/stable/config-file/v2.html for details
+
+# Required
+version: 2
+
+# Set the OS, Python version and other tools you might need
+build:
+  os: ubuntu-22.04
+  tools:
+    python: "3.12"
+    # You can also specify other tool versions:
+    # nodejs: "20"
+    # rust: "1.70"
+    # golang: "1.20"
+
+# Build documentation in the "docs/" directory with Sphinx
+sphinx:
+  configuration: docs/source/conf.py
+  # You can configure Sphinx to use a different builder, for instance use the dirhtml builder for simpler URLs
+  # builder: "dirhtml"
+  # Fail on all warnings to avoid broken references
+  # fail_on_warning: true
+
+# Optionally build your docs in additional formats such as PDF and ePub
+# formats:
+#   - pdf
+#   - epub
+
+# Optional but recommended, declare the Python requirements required
+# to build your documentation
+# See https://docs.readthedocs.io/en/stable/guides/reproducible-builds.html
+# python:
+#   install:
+#     - requirements: docs/requirements.txt
+
+python:
+  install:
+    - requirements: docs/requirements.txt

.vscode/c_cpp_properties.json

@@ -0,0 +1,20 @@
+{
+    "configurations": [
+        {
+            "name": "Mac",
+            "includePath": [
+                "${workspaceFolder}/**"
+            ],
+            "defines": [],
+            "macFrameworkPath": [
+                "/Applications/Xcode.app/Contents/Developer/Platforms/MacOSX.platform/Developer/SDKs/MacOSX.sdk/System/Library/Frameworks"
+            ],
+            "compilerPath": "/usr/bin/clang",
+            "cStandard": "c23",
+            "cppStandard": "c++23",
+            "intelliSenseMode": "macos-clang-arm64",
+            "configurationProvider": "vector-of-bool.cmake-tools"
+        }
+    ],
+    "version": 4
+}

.vscode/settings.json

@@ -1,4 +1,31 @@
 {
   "C_Cpp.default.configurationProvider": "vector-of-bool.cmake-tools",
-  "restructuredtext.confPath": "${workspaceFolder}/docs/source"
+  "restructuredtext.confPath": "${workspaceFolder}/docs/source",
+  "workbench.colorTheme": "Default Light+ Experimental",
+  "editor.fontSize": 14,
+  "workbench.colorCustomizations": {
+	"[Default Light*]": {
+  	    "editor.background": "#fefefe",
+  	    "sideBar.background": "#fefefe",
+  	    "sideBar.foreground": "#343436",
+  	    "sideBarTitle.foreground": "#343436",
+  	    "sideBar.border": "#e2e2e4",
+  	    "statusBar.background": "#fefefe",
+  	    "titleBar.activeBackground": "#fefefe",
+  	    "tab.activeBackground": "#f4fff4aa",
+  	    "tab.inactiveBackground": "#fefefe",
+  	    "activityBar.background": "#fefefe",
+  	    "editorGroupHeader.tabsBackground": "#fefefe"
+	},
+	"[Default Dark*]": {
+		"editor.background": "#1D1D25",
+		"sideBar.background": "#1D1D25",
+		"statusBar.background": "#1D1D25",
+		"titleBar.activeBackground": "#1D1D25",
+		"tab.activeBackground": "#2C2C3A",
+		"tab.inactiveBackground": "#1D1D25",
+		"activityBar.background": "#1D1D25",
+		"editorGroupHeader.tabsBackground": "#1D1D25"
+    }
+  },
 }

CMakeLists.txt

@@ -1,6 +1,6 @@
 cmake_minimum_required(VERSION 3.8.2)
 project(cglm
-  VERSION 0.8.8
+  VERSION 0.9.4
   HOMEPAGE_URL https://github.com/recp/cglm
   DESCRIPTION "OpenGL Mathematics (glm) for C"
   LANGUAGES C
@@ -16,6 +16,11 @@ option(CGLM_STATIC "Static build" OFF)
 option(CGLM_USE_C99 "" OFF)
 option(CGLM_USE_TEST "Enable Tests" OFF)
 
+if(CMAKE_SYSTEM_NAME STREQUAL WASI)
+  set(CGLM_STATIC ON CACHE BOOL "Static option" FORCE)
+  set(CGLM_SHARED OFF CACHE BOOL "Shared option" FORCE)
+endif()
+
 if(NOT CGLM_STATIC AND CGLM_SHARED)
   set(CGLM_BUILD SHARED)
 else(CGLM_STATIC)
@@ -27,21 +32,28 @@ if(CGLM_USE_C99)
 endif()
 
 if(MSVC)
-  add_definitions(-DNDEBUG -D_WINDOWS -D_USRDLL)
-  add_compile_options(/W3 /Ox /Gy /Oi /TC)
-  
-  # Ref: https://skia.googlesource.com/third_party/sdl/+/refs/heads/master/CMakeLists.txt#225
-  # Make sure /RTC1 is disabled, otherwise it will use functions from the CRT
-  foreach(flag_var
-      CMAKE_C_FLAGS CMAKE_C_FLAGS_DEBUG CMAKE_C_FLAGS_RELEASE
-      CMAKE_C_FLAGS_MINSIZEREL CMAKE_C_FLAGS_RELWITHDEBINFO)
-    string(REGEX REPLACE "/RTC(su|[1su])" "" ${flag_var} "${${flag_var}}")
-  endforeach(flag_var)
+  add_definitions(-D_WINDOWS -D_USRDLL)
+
+  if(NOT CMAKE_BUILD_TYPE MATCHES Debug)
+    add_definitions(-DNDEBUG)
+    add_compile_options(/W3 /Ox /Gy /Oi /TC)
+    foreach(flag_var
+        CMAKE_C_FLAGS CMAKE_C_FLAGS_DEBUG CMAKE_C_FLAGS_RELEASE
+        CMAKE_C_FLAGS_MINSIZEREL CMAKE_C_FLAGS_RELWITHDEBINFO)
+      string(REGEX REPLACE "/RTC(su|[1su])" "" ${flag_var} "${${flag_var}}")
+    endforeach(flag_var)
+  endif()
 else()
-  add_compile_options(-Wall -Werror -O3)
+  add_compile_options(-Wall)
+
+  if(NOT CMAKE_BUILD_TYPE MATCHES Debug)
+    add_compile_options(-O3)
+  endif()
 endif()
 
-if(NOT CMAKE_BUILD_TYPE AND NOT CMAKE_CONFIGURATION_TYPES)
+get_directory_property(hasParent PARENT_DIRECTORY)
+
+if(NOT hasParent AND NOT CMAKE_BUILD_TYPE AND NOT CMAKE_CONFIGURATION_TYPES)
   message(STATUS "Setting build type to '${DEFAULT_BUILD_TYPE}' as none was specified.")
   set(CMAKE_BUILD_TYPE "${DEFAULT_BUILD_TYPE}" CACHE STRING "Choose the type of build." FORCE)
   # Set the possible values of build type for cmake-gui
@@ -72,8 +84,14 @@ add_library(${PROJECT_NAME}
   src/vec4.c
   src/ivec4.c
   src/mat2.c
+  src/mat2x3.c
+  src/mat2x4.c
   src/mat3.c
+  src/mat3x2.c
+  src/mat3x4.c
   src/mat4.c
+  src/mat4x2.c
+  src/mat4x3.c
   src/plane.c
   src/frustum.c
   src/box.c
@@ -157,6 +175,20 @@ install(EXPORT      ${PROJECT_NAME}
         NAMESPACE   ${PROJECT_NAME}::
         DESTINATION ${CMAKE_INSTALL_LIBDIR}/cmake/${PROJECT_NAME})
 
+set(PACKAGE_NAME ${PROJECT_NAME})
+set(prefix ${CMAKE_INSTALL_PREFIX})
+set(exec_prefix ${CMAKE_INSTALL_PREFIX})
+if (IS_ABSOLUTE "${CMAKE_INSTALL_INCLUDEDIR}")
+  set(includedir "${CMAKE_INSTALL_INCLUDEDIR}")
+else()
+  set(includedir "\${prefix}/${CMAKE_INSTALL_INCLUDEDIR}")
+endif()
+if (IS_ABSOLUTE "${CMAKE_INSTALL_LIBDIR}")
+  set(libdir "${CMAKE_INSTALL_LIBDIR}")
+else()
+  set(libdir "\${exec_prefix}/${CMAKE_INSTALL_LIBDIR}")
+endif()
+set(PACKAGE_VERSION "${PROJECT_VERSION}")
 configure_file(cglm.pc.in cglm.pc @ONLY)
 
 install(FILES ${CMAKE_CURRENT_BINARY_DIR}/cglm.pc

CREDITS

@@ -82,3 +82,17 @@ http://github.com/microsoft/DirectXMath
 17. Pick Matrix
 
 glu project -> project.c
+
+18. Ray sphere intersection
+
+RAY TRACING GEMS
+HIGH-QUALITY AND REAL-TIME RENDERING WITH DXR AND OTHER APIS
+
+CHAPTER 7
+Precision Improvements for Ray/Sphere Intersection
+Eric Haines (1), Johannes Günther (2), and Tomas Akenine-Möller (1)
+  (1) NVIDIA
+  (2) Intel
+
+Wyman, C., and Haines, E. Getting Started with RTX Ray Tracing.
+https://github.com/NVIDIAGameWorks/GettingStartedWithRTXRayTracing

Makefile.am

@@ -11,8 +11,7 @@ AM_CFLAGS = -Wall \
             -std=gnu11 \
             -O3 \
             -Wstrict-aliasing=2 \
-            -fstrict-aliasing \
-            -Werror=strict-prototypes
+            -fstrict-aliasing
 
 lib_LTLIBRARIES = libcglm.la
 libcglm_la_LDFLAGS = -no-undefined -version-info 0:1:0
@@ -43,8 +42,14 @@ cglm_HEADERS = include/cglm/version.h \
                include/cglm/cam.h \
                include/cglm/io.h \
                include/cglm/mat4.h \
+               include/cglm/mat4x2.h \
+               include/cglm/mat4x3.h \
                include/cglm/mat3.h \
+               include/cglm/mat3x2.h \
+               include/cglm/mat3x4.h \
                include/cglm/mat2.h \
+               include/cglm/mat2x3.h \
+               include/cglm/mat2x4.h \
                include/cglm/affine-pre.h \
                include/cglm/affine-post.h \
                include/cglm/affine.h \
@@ -64,6 +69,7 @@ cglm_HEADERS = include/cglm/version.h \
                include/cglm/plane.h \
                include/cglm/frustum.h \
                include/cglm/box.h \
+               include/cglm/aabb2d.h \
                include/cglm/color.h \
                include/cglm/project.h \
                include/cglm/sphere.h \
@@ -95,8 +101,14 @@ cglm_clipspace_HEADERS = include/cglm/clipspace/persp.h \
 
 cglm_calldir=$(includedir)/cglm/call
 cglm_call_HEADERS = include/cglm/call/mat4.h \
+                    include/cglm/call/mat4x2.h \
+                    include/cglm/call/mat4x3.h \
                     include/cglm/call/mat3.h \
+                    include/cglm/call/mat3x2.h \
+                    include/cglm/call/mat3x4.h \
                     include/cglm/call/mat2.h \
+                    include/cglm/call/mat2x3.h \
+                    include/cglm/call/mat2x4.h \
                     include/cglm/call/vec2.h \
                     include/cglm/call/vec3.h \
                     include/cglm/call/vec4.h \
@@ -138,7 +150,8 @@ cglm_call_clipspace_HEADERS = include/cglm/call/clipspace/persp_lh_no.h \
 cglm_simddir=$(includedir)/cglm/simd
 cglm_simd_HEADERS = include/cglm/simd/intrin.h \
                     include/cglm/simd/x86.h \
-                    include/cglm/simd/arm.h
+                    include/cglm/simd/arm.h \
+                    include/cglm/simd/wasm.h
 
 cglm_simd_sse2dir=$(includedir)/cglm/simd/sse2
 cglm_simd_sse2_HEADERS = include/cglm/simd/sse2/affine.h \
@@ -157,20 +170,41 @@ cglm_simd_neon_HEADERS = include/cglm/simd/neon/affine.h \
                          include/cglm/simd/neon/mat4.h \
                          include/cglm/simd/neon/quat.h
 
+cglm_simd_wasmdir=$(includedir)/cglm/simd/wasm
+cglm_simd_wasm_HEADERS = include/cglm/simd/wasm/affine.h \
+                         include/cglm/simd/wasm/mat2.h \
+                         include/cglm/simd/wasm/mat3.h \
+                         include/cglm/simd/wasm/mat4.h \
+                         include/cglm/simd/wasm/quat.h
+
+cglm_handeddir=$(includedir)/cglm/handed
+cglm_handed_HEADERS = include/cglm/handed/euler_to_quat_lh.h \
+                      include/cglm/handed/euler_to_quat_rh.h
+
 cglm_structdir=$(includedir)/cglm/struct
 cglm_struct_HEADERS = include/cglm/struct/mat4.h \
+                      include/cglm/struct/mat4x2.h \
+                      include/cglm/struct/mat4x3.h \
                       include/cglm/struct/mat3.h \
+                      include/cglm/struct/mat3x2.h \
+                      include/cglm/struct/mat3x4.h \
                       include/cglm/struct/mat2.h \
+                      include/cglm/struct/mat2x3.h \
+                      include/cglm/struct/mat2x4.h \
                       include/cglm/struct/affine-pre.h \
                       include/cglm/struct/affine-post.h \
+                      include/cglm/struct/affine-mat.h \
                       include/cglm/struct/affine.h \
                       include/cglm/struct/affine2d.h \
                       include/cglm/struct/vec2.h \
                       include/cglm/struct/vec2-ext.h \
+                      include/cglm/struct/ivec2.h \
                       include/cglm/struct/vec3.h \
                       include/cglm/struct/vec3-ext.h \
+                      include/cglm/struct/ivec3.h \
                       include/cglm/struct/vec4.h \
                       include/cglm/struct/vec4-ext.h \
+                      include/cglm/struct/ivec4.h \
                       include/cglm/struct/io.h \
                       include/cglm/struct/cam.h \
                       include/cglm/struct/quat.h \
@@ -178,10 +212,12 @@ cglm_struct_HEADERS = include/cglm/struct/mat4.h \
                       include/cglm/struct/plane.h \
                       include/cglm/struct/frustum.h \
                       include/cglm/struct/box.h \
+                      include/cglm/struct/aabb2d.h \
                       include/cglm/struct/project.h \
                       include/cglm/struct/sphere.h \
                       include/cglm/struct/color.h \
-                      include/cglm/struct/curve.h
+                      include/cglm/struct/curve.h \
+                      include/cglm/struct/ray.h
 
 cglm_struct_clipspacedir=$(includedir)/cglm/struct/clipspace
 cglm_struct_clipspace_HEADERS = include/cglm/struct/clipspace/persp_lh_no.h \
@@ -199,6 +235,10 @@ cglm_struct_clipspace_HEADERS = include/cglm/struct/clipspace/persp_lh_no.h \
                                 include/cglm/struct/clipspace/project_no.h \
                                 include/cglm/struct/clipspace/project_zo.h
 
+cglm_struct_handeddir=$(includedir)/cglm/struct/handed
+cglm_struct_handed_HEADERS = include/cglm/struct/handed/euler_to_quat_lh.h \
+                             include/cglm/struct/handed/euler_to_quat_rh.h
+
 libcglm_la_SOURCES=\
     src/euler.c \
     src/affine.c \
@@ -212,8 +252,14 @@ libcglm_la_SOURCES=\
     src/vec4.c \
     src/ivec4.c \
     src/mat2.c \
+    src/mat2x3.c \
+    src/mat2x4.c \
     src/mat3.c \
+    src/mat3x2.c \
+    src/mat3x4.c \
     src/mat4.c \
+    src/mat4x2.c \
+    src/mat4x3.c \
     src/plane.c \
     src/frustum.c \
     src/box.c \
@@ -243,11 +289,6 @@ test_tests_SOURCES=\
     test/runner.c \
     test/src/test_common.c \
     test/src/tests.c \
-    test/src/test_cam.c \
-    test/src/test_cam_lh_zo.c \
-    test/src/test_cam_rh_zo.c \
-    test/src/test_cam_lh_no.c \
-    test/src/test_cam_rh_no.c \
     test/src/test_clamp.c \
     test/src/test_euler.c \
     test/src/test_bezier.c \

README.md

@@ -5,14 +5,10 @@
 </p>
 <br>
 <p align="center">
-    <a href="https://travis-ci.com/recp/cglm">
-        <img src="https://travis-ci.com/recp/cglm.svg?branch=master"
+    <a href="https://github.com/recp/cglm/actions/workflows/ci.yml">
+        <img src="https://github.com/recp/cglm/actions/workflows/ci.yml/badge.svg"
              alt="Build Status">
     </a>
-    <a href="https://ci.appveyor.com/project/recp/cglm/branch/master">
-        <img src="https://ci.appveyor.com/api/projects/status/av7l3gc0yhfex8y4/branch/master?svg=true"
-             alt="Windows Build Status">
-    </a>
     <a href="http://cglm.readthedocs.io/en/latest/?badge=latest">
         <img src="https://readthedocs.org/projects/cglm/badge/?version=latest"
              alt="Documentation Status">
@@ -66,7 +62,7 @@ you have the latest version
 - **[major change]** by starting v0.5.0, vec3 functions use **glm_vec3_** namespace, it was **glm_vec_** until v0.5.0
 - **[major change]** by starting v0.5.1, built-in alignment is removed from **vec3** and **mat3** types
 - **[major change]** by starting v0.7.3, inline print functions are disabled in release/production mode to eliminate print costs (see options in documentation). Print output also improved. You can disable colors if you need  (see documentation)
-- **[major change]** by starting v0.8.3, **cglm** supports alternative clipspace configuations e.g. Left Handed, Zero-to-One (_zo)... `CGLM_FORCE_DEPTH_ZERO_TO_ONE` and `CGLM_FORCE_LEFT_HANDED` is provided to control clipspace. You should be able to use **cglm** with Vulkan, DirectX and Metal now... see https://cglm.readthedocs.io/en/latest/opt.html#clipspace-option-s
+- **[major change]** by starting v0.8.3, **cglm** supports alternative clipspace configurations e.g. Left Handed, Zero-to-One (_zo)... `CGLM_FORCE_DEPTH_ZERO_TO_ONE` and `CGLM_FORCE_LEFT_HANDED` is provided to control clipspace. You should be able to use **cglm** with Vulkan, DirectX and Metal now... see https://cglm.readthedocs.io/en/latest/opt.html#clipspace-option-s
 
 #### 📌 Note for C++ developers:
 If you are not aware of the original GLM library yet, you may also want to look at:
@@ -108,7 +104,7 @@ https://github.com/g-truc/glm
 
 ## 🚀 Features
 - **scalar** and **simd** (sse, avx, neon...) optimizations
-- option to use different clipspaces e.g. Left Handed, Zero-to-One... (currrently right handed negative-one is default)
+- option to use different clipspaces e.g. Left Handed, Zero-to-One... (currently right handed negative-one is default)
 - array api and struct api, you can use arrays or structs.
 - general purpose matrix operations (mat4, mat3)
 - chain matrix multiplication (square only)
@@ -148,7 +144,7 @@ To call pre-compiled versions, just use `glmc_` (c stands for 'call') instead of
   glm_mul(trans, rot, rt);  /* inline */
   glmc_mul(trans, rot, rt); /* call from library */
 ```
-Most of math functions are optimized manualy with SSE2 if available, if not? Dont worry there are non-sse versions of all operations
+Most of math functions are optimized manually with SSE2 if available, if not? Dont worry there are non-sse versions of all operations
 
 You can pass matrices and vectors as array to functions rather than get address.
 
@@ -242,6 +238,38 @@ add_subdirectory(external/cglm/)
 # or you can use find_package to configure cglm
 ```
 
+#### Use CMake to build for WebAssembly
+
+Since math functions like `sinf` is used, this can not be targeted at `wasm32-unknown-unknown`, one of [wasi-sdk](https://github.com/WebAssembly/wasi-sdk) or [emscripten](https://github.com/emscripten-core/emsdk) should be used.
+
+Should note that shared build is not yet supported for WebAssembly.
+
+For [simd128](https://github.com/WebAssembly/simd) support, add `-msimd128` to `CMAKE_C_FLAGS`, in command line `-DCMAKE_C_FLAGS="-msimd128"`.
+
+For tests, the cmake option `CGLM_USE_TEST` would still work, you'll need a wasi runtime for running tests, see our [ci config file](.github/workflows/cmake-wasm.yml) for a detailed example.
+
+##### Use CMake and WASI SDK to build for WebAssembly
+
+```bash
+$ cmake .. \
+  -DCMAKE_TOOLCHAIN_FILE=/path/to/wasi-sdk-19.0/share/cmake/wasi-sdk.cmake \
+  -DWASI_SDK_PREFIX=/path/to/wasi-sdk-19.0
+```
+
+Where `/path/to/wasi-sdk-19.0/` is the path to extracted [wasi sdk](https://github.com/WebAssembly/wasi-sdk).
+
+In this case it would by default make a static build.
+
+##### Use CMake and Emscripten SDK to build for WebAssembly
+
+```bash
+$ emcmake cmake .. \
+  -DCMAKE_EXE_LINKER_FLAGS="-s STANDALONE_WASM" \
+  -DCGLM_STATIC=ON
+```
+
+The `emcmake` here is the cmake wrapper for Emscripten from installed [emsdk](https://github.com/emscripten-core/emsdk).
+
 ### Meson (All platforms)
 
 ```bash
@@ -257,7 +285,7 @@ $ sudo ninja install # [Optional]
 c_std=c11
 buildtype=release
 default_library=shared
-enable_tests=false # to run tests: ninja test
+build_tests=true # to run tests: ninja test
 ```
 #### Use with your Meson project
 * Example:

appveyor.yml

@@ -1,7 +0,0 @@
-image: Visual Studio 2017
-
-build_script:
-- ps: >-
-    cd win
-
-    .\build.bat

cglm.pc.in

@@ -1,11 +1,11 @@
-prefix=@CMAKE_INSTALL_PREFIX@
-exec_prefix="${prefix}"
-libdir="${exec_prefix}/lib"
-includedir="${prefix}/include"
+prefix=@prefix@
+exec_prefix=@exec_prefix@
+libdir=@libdir@
+includedir=@includedir@
 
-Name: @PROJECT_NAME@
-Description: @CMAKE_PROJECT_DESCRIPTION@
-URL: @CMAKE_PROJECT_HOMEPAGE_URL@
-Version: @PROJECT_VERSION@
+Name: @PACKAGE_NAME@
+Description: OpenGL Mathematics (glm) for C
+URL: https://github.com/recp/cglm
+Version: @PACKAGE_VERSION@
 Cflags: -I${includedir}
 Libs: -L${libdir} -lcglm @LIBS@

cglm.podspec

@@ -2,7 +2,7 @@ Pod::Spec.new do |s|
 
   # Description
   s.name         = "cglm"
-  s.version      = "0.8.7"
+  s.version      = "0.9.3"
   s.summary      = "📽 Highly Optimized Graphics Math (glm) for C"
   s.description  = <<-DESC
 cglm is math library for graphics programming for C. See the documentation or README for all features.

configure.ac

@@ -7,8 +7,8 @@
 #*****************************************************************************
 
 AC_PREREQ([2.69])
-AC_INIT([cglm], [0.8.8], [info@recp.me])
-AM_INIT_AUTOMAKE([-Wall -Werror foreign subdir-objects serial-tests])
+AC_INIT([cglm], [0.9.4], [info@recp.me])
+AM_INIT_AUTOMAKE([-Wall foreign subdir-objects serial-tests])
 
 # Don't use the default cflags (-O2 -g), we set ours manually in Makefile.am.
 : ${CFLAGS=""}

docs/requirements.txt

@@ -0,0 +1,4 @@
+# Defining the exact version will make sure things don't break
+sphinx==7.2.6
+sphinx_rtd_theme==2.0.0
+readthedocs-sphinx-search==0.3.2

docs/source/aabb2d.rst

@@ -0,0 +1,198 @@
+.. default-domain:: C
+
+2d axis aligned bounding box (AABB)
+================================================================================
+
+Header: cglm/aabb2d.h
+
+Some convenient functions provided for AABB.
+
+**Definition of aabb:**
+
+cglm defines an aabb as a two dimensional array of vec2's.
+The first element is the **min** point and the second one is the **max** point.
+If you have another type e.g. struct or even another representation then you must
+convert it before and after calling a cglm aabb2d function.
+
+Table of contents (click to go):
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+Macros:
+
+1. :c:func:`glm_aabb2d_size`
+
+Functions:
+
+1. :c:func:`glm_aabb2d_copy`
+#. :c:func:`glm_aabb2d_zero`
+#. :c:func:`glm_aabb2d_transform`
+#. :c:func:`glm_aabb2d_merge`
+#. :c:func:`glm_aabb2d_crop`
+#. :c:func:`glm_aabb2d_crop_until`
+#. :c:func:`glm_aabb2d_invalidate`
+#. :c:func:`glm_aabb2d_isvalid`
+#. :c:func:`glm_aabb2d_diag`
+#. :c:func:`glm_aabb2d_sizev`
+#. :c:func:`glm_aabb2d_radius`
+#. :c:func:`glm_aabb2d_center`
+#. :c:func:`glm_aabb2d_aabb`
+#. :c:func:`glm_aabb2d_circle`
+#. :c:func:`glm_aabb2d_point`
+#. :c:func:`glm_aabb2d_contains`
+
+Functions documentation
+~~~~~~~~~~~~~~~~~~~~~~~
+
+.. c:function:: void  glm_aabb2d_copy(vec2 aabb[2], vec2 dest[2])
+
+    | copy all members of [aabb] to [dest]
+
+    Parameters:
+      | *[in]*  **aabb**  bounding box
+      | *[out]* **dest**  destination
+
+.. c:function:: void  glm_aabb2d_zero(vec2 aabb[2])
+
+    | makes all members of [aabb] 0.0f (zero)
+
+    Parameters:
+      | *[in, out]*  **aabb**     bounding box
+
+.. c:function:: void  glm_aabb2d_transform(vec2 aabb[2], mat3 m, vec2 dest[2])
+
+    | apply transform to Axis-Aligned Bounding Box
+
+    Parameters:
+      | *[in]*  **aabb**   bounding box
+      | *[in]*  **m**     transform matrix
+      | *[out]* **dest**  transformed bounding box
+
+.. c:function:: void  glm_aabb2d_merge(vec2 aabb1[2], vec2 aabb2[2], vec2 dest[2])
+
+    | merges two AABB bounding box and creates new one
+
+    two aabb must be in the same space
+
+    Parameters:
+      | *[in]*  **aabb1** bounding box 1
+      | *[in]*  **aabb2** bounding box 2
+      | *[out]* **dest** merged bounding box
+
+.. c:function:: void  glm_aabb2d_crop(vec2 aabb[2], vec2 cropAabb[2], vec2 dest[2])
+
+    | crops a bounding box with another one.
+
+    this could be useful for gettng a bbox which fits with view frustum and
+    object bounding boxes. In this case you crop view frustum box with objects
+    box
+
+    Parameters:
+      | *[in]*  **aabb**      bounding box 1
+      | *[in]*  **cropAabb**  crop box
+      | *[out]* **dest**     cropped bounding box
+
+.. c:function:: void  glm_aabb2d_crop_until(vec2 aabb[2], vec2 cropAabb[2], vec2 clampAabb[2], vec2 dest[2])
+
+    | crops a bounding box with another one.
+
+    this could be useful for gettng a bbox which fits with view frustum and
+    object bounding boxes. In this case you crop view frustum box with objects
+    box
+
+    Parameters:
+      | *[in]*  **aabb**      bounding box
+      | *[in]*  **cropAabb**  crop box
+      | *[in]*  **clampAabb** minimum box
+      | *[out]* **dest**     cropped bounding box
+
+.. c:function:: void  glm_aabb2d_invalidate(vec2 aabb[2])
+
+    | invalidate AABB min and max values
+
+    | It fills *max* values with -FLT_MAX and *min* values with +FLT_MAX
+
+    Parameters:
+      | *[in, out]*   **aabb**     bounding box
+
+.. c:function:: bool  glm_aabb2d_isvalid(vec2 aabb[2])
+
+    | check if AABB is valid or not
+
+    Parameters:
+      | *[in]*   **aabb**     bounding box
+
+    Returns:
+      returns true if aabb is valid otherwise false
+
+.. c:function:: float  glm_aabb2d_diag(vec2 aabb[2])
+
+    | distance between min and max
+
+    Parameters:
+      | *[in]*   **aabb**     bounding box
+
+    Returns:
+      distance between min - max
+
+
+.. c:function:: void  glm_aabb2d_sizev(vec2 aabb[2], vec2 dest)
+
+    | size vector of aabb
+
+    Parameters:
+      | *[in]*   **aabb**     bounding box
+      | *[out]*  **dest**     size vector
+
+    Returns:
+      size vector of aabb max - min
+
+.. c:function:: float  glm_aabb2d_radius(vec2 aabb[2])
+
+    | radius of sphere which surrounds AABB
+
+    Parameters:
+      | *[in]*   **aabb**     bounding box
+
+.. c:function:: void  glm_aabb2d_center(vec2 aabb[2], vec2 dest)
+
+    | computes center point of AABB
+
+    Parameters:
+      | *[in]*    **aabb**      bounding box
+      | *[out]*   **dest**     center of bounding box
+
+.. c:function:: bool  glm_aabb2d_aabb(vec2 aabb[2], vec2 other[2])
+
+    | check if two AABB intersects
+
+    Parameters:
+      | *[in]*    **aabb**     bounding box
+      | *[out]*   **other**   other bounding box
+
+.. c:function:: bool  glm_aabb2d_circle(vec2 aabb[2], vec3 c)
+
+    | check if AABB intersects with sphere
+
+    | https://github.com/erich666/GraphicsGems/blob/master/gems/BoxSphere.c
+    | Solid Box - Solid Sphere test.
+
+    Parameters:
+      | *[in]*    **aabb**     solid bounding box
+      | *[out]*   **c**        solid circle
+
+.. c:function:: bool  glm_aabb2d_point(vec2 aabb[2], vec2 point)
+
+    | check if point is inside of AABB
+
+    Parameters:
+      | *[in]*    **aabb**     bounding box
+      | *[out]*   **point**   point
+
+.. c:function:: bool  glm_aabb2d_contains(vec2 aabb[2], vec2 other[2])
+
+    | check if AABB contains other AABB
+
+    Parameters:
+      | *[in]*    **aabb**     bounding box
+      | *[out]*   **other**   other bounding box
+

docs/source/affine-common.rst

@@ -3,7 +3,7 @@
 3D Affine Transforms (common)
 ================================================================================
 
-Common transfrom functions.
+Common transform functions.
 
 Table of contents (click to go):
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
@@ -30,7 +30,7 @@ Functions documentation
     creates NEW translate transform matrix by *v* vector.
 
     Parameters:
-      | *[in, out]* **m**  affine transfrom
+      | *[in, out]* **m**  affine transform
       | *[in]*      **v**  translate vector [x, y, z]
 
 .. c:function:: void  glm_scale_to(mat4 m, vec3 v, mat4 dest)
@@ -38,7 +38,7 @@ Functions documentation
     scale existing transform matrix by *v* vector and store result in dest
 
     Parameters:
-      | *[in]*  **m**    affine transfrom
+      | *[in]*  **m**    affine transform
       | *[in]*  **v**    scale vector [x, y, z]
       | *[out]* **dest** scaled matrix
 
@@ -47,7 +47,7 @@ Functions documentation
     creates NEW scale matrix by v vector
 
     Parameters:
-      | *[out]* **m** affine transfrom
+      | *[out]* **m** affine transform
       | *[in]*  **v** scale vector [x, y, z]
 
 .. c:function:: void  glm_scale(mat4 m, vec3 v)
@@ -56,7 +56,7 @@ Functions documentation
     and stores result in same matrix
 
     Parameters:
-      | *[in, out]* **m** affine transfrom
+      | *[in, out]* **m** affine transform
       | *[in]*      **v** scale vector [x, y, z]
 
 .. c:function:: void  glm_scale_uni(mat4 m, float s)
@@ -65,7 +65,7 @@ Functions documentation
     and stores result in same matrix
 
     Parameters:
-      | *[in, out]* **m** affine transfrom
+      | *[in, out]* **m** affine transform
       | *[in]*      **v** scale factor
 
 .. c:function:: void  glm_rotate_make(mat4 m, float angle, vec3 axis)
@@ -74,7 +74,7 @@ Functions documentation
     axis will be normalized so you don't need to normalize it
 
     Parameters:
-      | *[out]* **m**    affine transfrom
+      | *[out]* **m**    affine transform
       | *[in]*  **axis** angle (radians)
       | *[in]*  **axis** axis
 
@@ -86,7 +86,7 @@ Functions documentation
     | this should work faster than glm_rotate_at because it reduces one glm_translate.
 
     Parameters:
-      | *[in, out]* **m**     affine transfrom
+      | *[in, out]* **m**     affine transform
       | *[in]*      **pivot** pivot, anchor point, rotation center
       | *[in]*      **angle** angle (radians)
       | *[in]*      **axis**  axis
@@ -123,7 +123,7 @@ Functions documentation
     DON'T pass projected matrix here
 
     Parameters:
-      | *[in]*  **m** affine transfrom
+      | *[in]*  **m** affine transform
       | *[out]* **t** translation vector
       | *[out]* **r** rotation matrix (mat4)
       | *[out]* **s** scaling vector [X, Y, Z]

docs/source/affine-mat.rst

@@ -25,7 +25,7 @@ You cannot use :c:func:`glm_mul` anymore.
 Same is also true for :c:func:`glm_inv_tr` if you only have rotation and
 translation then it will work as expected, otherwise you cannot use that.
 
-In the future it may accept scale factors too but currectly it does not.
+In the future it may accept scale factors too but currently it does not.
 
 Table of contents (click func go):
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

docs/source/affine-post.rst

@@ -3,10 +3,8 @@
 3D Affine Transforms (post)
 ================================================================================
 
-Post transfrom functions are similar to pre transform functions except order of application is reversed.  
-Post transform functions are applied after the object is transformed with given (model matrix) transfrom. 
-
-Ther are named af
+Post transform functions are similar to pre transform functions except order of application is reversed.
+Post transform functions are applied after the object is transformed with given (model matrix) transform.
 
 Table of contents (click to go):
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
@@ -33,7 +31,7 @@ Functions documentation
     translate existing transform matrix by *v* vector and store result in dest
 
     Parameters:
-      | *[in]*  **m**    affine transfrom
+      | *[in]*  **m**    affine transform
       | *[in]*  **v**    translate vector [x, y, z]
       | *[out]* **dest** translated matrix
 
@@ -43,7 +41,7 @@ Functions documentation
     and stores result in same matrix
 
     Parameters:
-      | *[in, out]* **m**  affine transfrom
+      | *[in, out]* **m**  affine transform
       | *[in]*      **v**  translate vector [x, y, z]
 
 .. c:function:: void  glm_translated_x(mat4 m, float x)
@@ -51,7 +49,7 @@ Functions documentation
     translate existing transform matrix by x factor
 
     Parameters:
-      | *[in, out]* **m**  affine transfrom
+      | *[in, out]* **m**  affine transform
       | *[in]*      **v**  x factor
 
 .. c:function:: void  glm_translated_y(mat4 m, float y)
@@ -59,7 +57,7 @@ Functions documentation
     translate existing transform matrix by *y* factor
 
     Parameters:
-      | *[in, out]* **m**  affine transfrom
+      | *[in, out]* **m**  affine transform
       | *[in]*      **v**  y factor
 
 .. c:function:: void  glm_translated_z(mat4 m, float z)
@@ -67,7 +65,7 @@ Functions documentation
     translate existing transform matrix by *z* factor
 
     Parameters:
-      | *[in, out]* **m**  affine transfrom
+      | *[in, out]* **m**  affine transform
       | *[in]*      **v**  z factor
 
 .. c:function:: void  glm_rotated_x(mat4 m, float angle, mat4 dest)
@@ -76,7 +74,7 @@ Functions documentation
     and store result in dest
 
     Parameters:
-      | *[in]*  **m**     affine transfrom
+      | *[in]*  **m**     affine transform
       | *[in]*  **angle** angle (radians)
       | *[out]* **dest**  rotated matrix
 
@@ -86,7 +84,7 @@ Functions documentation
     and store result in dest
 
     Parameters:
-      | *[in]*  **m**     affine transfrom
+      | *[in]*  **m**     affine transform
       | *[in]*  **angle** angle (radians)
       | *[out]* **dest**  rotated matrix
 
@@ -96,7 +94,7 @@ Functions documentation
     and store result in dest
 
     Parameters:
-      | *[in]*  **m**     affine transfrom
+      | *[in]*  **m**     affine transform
       | *[in]*  **angle** angle (radians)
       | *[out]* **dest**  rotated matrix
 
@@ -105,7 +103,7 @@ Functions documentation
     rotate existing transform matrix around Z axis by angle and axis
 
     Parameters:
-      | *[in, out]* **m**     affine transfrom
+      | *[in, out]* **m**     affine transform
       | *[in]*      **angle** angle (radians)
       | *[in]*      **axis**  axis
 
@@ -114,7 +112,7 @@ Functions documentation
     rotate existing transform around given axis by angle at given pivot point (rotation center)
 
     Parameters:
-      | *[in, out]* **m**     affine transfrom
+      | *[in, out]* **m**     affine transform
       | *[in]*      **pivot** pivot, anchor point, rotation center
       | *[in]*      **angle** angle (radians)
       | *[in]*      **axis**  axis
@@ -124,6 +122,6 @@ Functions documentation
     | rotate existing transform matrix around given axis by angle around self (doesn't affected by position)
 
     Parameters:
-      | *[in, out]* **m**     affine transfrom
+      | *[in, out]* **m**     affine transform
       | *[in]*      **angle** angle (radians)
       | *[in]*      **axis**  axis

docs/source/affine-pre.rst

@@ -3,7 +3,7 @@
 3D Affine Transforms (pre)
 ================================================================================
 
-Pre transfrom functions which are regular transfrom functions.
+Pre transform functions which are regular transform functions.
 
 Table of contents (click to go):
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
@@ -41,7 +41,7 @@ Functions documentation
     translate existing transform matrix by *v* vector and store result in dest
 
     Parameters:
-      | *[in]*  **m**    affine transfrom
+      | *[in]*  **m**    affine transform
       | *[in]*  **v**    translate vector [x, y, z]
       | *[out]* **dest** translated matrix
 
@@ -51,7 +51,7 @@ Functions documentation
     and stores result in same matrix
 
     Parameters:
-      | *[in, out]* **m**  affine transfrom
+      | *[in, out]* **m**  affine transform
       | *[in]*      **v**  translate vector [x, y, z]
 
 .. c:function:: void  glm_translate_x(mat4 m, float x)
@@ -59,7 +59,7 @@ Functions documentation
     translate existing transform matrix by x factor
 
     Parameters:
-      | *[in, out]* **m**  affine transfrom
+      | *[in, out]* **m**  affine transform
       | *[in]*      **v**  x factor
 
 .. c:function:: void  glm_translate_y(mat4 m, float y)
@@ -67,7 +67,7 @@ Functions documentation
     translate existing transform matrix by *y* factor
 
     Parameters:
-      | *[in, out]* **m**  affine transfrom
+      | *[in, out]* **m**  affine transform
       | *[in]*      **v**  y factor
 
 .. c:function:: void  glm_translate_z(mat4 m, float z)
@@ -75,59 +75,16 @@ Functions documentation
     translate existing transform matrix by *z* factor
 
     Parameters:
-      | *[in, out]* **m**  affine transfrom
+      | *[in, out]* **m**  affine transform
       | *[in]*      **v**  z factor
 
-.. c:function:: void  glm_translate_make(mat4 m, vec3 v)
-
-    creates NEW translate transform matrix by *v* vector.
-
-    Parameters:
-      | *[in, out]* **m**  affine transfrom
-      | *[in]*      **v**  translate vector [x, y, z]
-
-.. c:function:: void  glm_scale_to(mat4 m, vec3 v, mat4 dest)
-
-    scale existing transform matrix by *v* vector and store result in dest
-
-    Parameters:
-      | *[in]*  **m**    affine transfrom
-      | *[in]*  **v**    scale vector [x, y, z]
-      | *[out]* **dest** scaled matrix
-
-.. c:function:: void  glm_scale_make(mat4 m, vec3 v)
-
-    creates NEW scale matrix by v vector
-
-    Parameters:
-      | *[out]* **m** affine transfrom
-      | *[in]*  **v** scale vector [x, y, z]
-
-.. c:function:: void  glm_scale(mat4 m, vec3 v)
-
-    scales existing transform matrix by v vector
-    and stores result in same matrix
-
-    Parameters:
-      | *[in, out]* **m** affine transfrom
-      | *[in]*      **v** scale vector [x, y, z]
-
-.. c:function:: void  glm_scale_uni(mat4 m, float s)
-
-    applies uniform scale to existing transform matrix v = [s, s, s]
-    and stores result in same matrix
-
-    Parameters:
-      | *[in, out]* **m** affine transfrom
-      | *[in]*      **v** scale factor
-
 .. c:function:: void  glm_rotate_x(mat4 m, float angle, mat4 dest)
 
     rotate existing transform matrix around X axis by angle
     and store result in dest
 
     Parameters:
-      | *[in]*  **m**     affine transfrom
+      | *[in]*  **m**     affine transform
       | *[in]*  **angle** angle (radians)
       | *[out]* **dest**  rotated matrix
 
@@ -137,7 +94,7 @@ Functions documentation
     and store result in dest
 
     Parameters:
-      | *[in]*  **m**     affine transfrom
+      | *[in]*  **m**     affine transform
       | *[in]*  **angle** angle (radians)
       | *[out]* **dest**  rotated matrix
 
@@ -147,26 +104,33 @@ Functions documentation
     and store result in dest
 
     Parameters:
-      | *[in]*  **m**     affine transfrom
+      | *[in]*  **m**     affine transform
       | *[in]*  **angle** angle (radians)
       | *[out]* **dest**  rotated matrix
 
-.. c:function:: void  glm_rotate_make(mat4 m, float angle, vec3 axis)
-
-    creates NEW rotation matrix by angle and axis,
-    axis will be normalized so you don't need to normalize it
-
-    Parameters:
-      | *[out]* **m**    affine transfrom
-      | *[in]*  **axis** angle (radians)
-      | *[in]*  **axis** axis
-
 .. c:function:: void  glm_rotate(mat4 m, float angle, vec3 axis)
 
-    rotate existing transform matrix around Z axis by angle and axis
+    rotate existing transform matrix around given axis by angle at ORIGIN (0,0,0)
+
+    **❗️IMPORTANT ❗️**
+    
+    If you need to rotate object around itself e.g. center of object or at
+    some point [of object] then `glm_rotate_at()` would be better choice to do so.
+    
+    Even if object's model transform is identiy, rotation may not be around
+    center of object if object does not lay out at ORIGIN perfectly.
+    
+    Using `glm_rotate_at()` with center of bounding shape ( AABB, Sphere ... )
+    would be an easy option to rotate around object if object is not at origin.
+    
+    One another option to rotate around itself at any point is `glm_spin()`
+    which is perfect if only rotating around model position is desired e.g. not
+    specific point on model for instance center of geometry or center of mass,
+    again if geometry is not perfectly centered at origin at identity transform,
+    rotation may not be around geometry.
 
     Parameters:
-      | *[in, out]* **m**     affine transfrom
+      | *[in, out]* **m**     affine transform
       | *[in]*      **angle** angle (radians)
       | *[in]*      **axis**  axis
 
@@ -175,66 +139,16 @@ Functions documentation
     rotate existing transform around given axis by angle at given pivot point (rotation center)
 
     Parameters:
-      | *[in, out]* **m**     affine transfrom
+      | *[in, out]* **m**     affine transform
       | *[in]*      **pivot** pivot, anchor point, rotation center
       | *[in]*      **angle** angle (radians)
       | *[in]*      **axis**  axis
 
-.. c:function:: void  glm_rotate_atm(mat4 m, vec3 pivot, float angle, vec3 axis)
-
-    | creates NEW rotation matrix by angle and axis at given point
-    | this creates rotation matrix, it assumes you don't have a matrix
-
-    | this should work faster than glm_rotate_at because it reduces one glm_translate.
-
-    Parameters:
-      | *[in, out]* **m**     affine transfrom
-      | *[in]*      **pivot** pivot, anchor point, rotation center
-      | *[in]*      **angle** angle (radians)
-      | *[in]*      **axis**  axis
-
-.. c:function:: void  glm_decompose_scalev(mat4 m, vec3 s)
-
-    decompose scale vector
-
-    Parameters:
-      | *[in]*  **m**  affine transform
-      | *[out]* **s**  scale vector (Sx, Sy, Sz)
-
-.. c:function:: bool  glm_uniscaled(mat4 m)
-
-    returns true if matrix is uniform scaled.
-    This is helpful for creating normal matrix.
-
-    Parameters:
-      | *[in]*  **m**   matrix
-
-.. c:function:: void  glm_decompose_rs(mat4 m, mat4 r, vec3 s)
-
-    decompose rotation matrix (mat4) and scale vector [Sx, Sy, Sz]
-    DON'T pass projected matrix here
-
-    Parameters:
-      | *[in]*  **m** affine transform
-      | *[out]* **r** rotation matrix
-      | *[out]* **s** scale matrix
-
-.. c:function:: void  glm_decompose(mat4 m, vec4 t, mat4 r, vec3 s)
-
-    decompose affine transform, TODO: extract shear factors.
-    DON'T pass projected matrix here
-
-    Parameters:
-      | *[in]*  **m** affine transfrom
-      | *[out]* **t** translation vector
-      | *[out]* **r** rotation matrix (mat4)
-      | *[out]* **s** scaling vector [X, Y, Z]
-
 .. c:function:: void  glm_spin(mat4 m, float angle, vec3 axis)
 
     | rotate existing transform matrix around given axis by angle around self (doesn't affected by position)
 
     Parameters:
-      | *[in, out]* **m**     affine transfrom
+      | *[in, out]* **m**     affine transform
       | *[in]*      **angle** angle (radians)
       | *[in]*      **axis**  axis

docs/source/affine.rst

@@ -13,7 +13,7 @@ Post functions (`T' = T * Tnew`) are like `glm_translated`, `glm_rotated` which
 `glm_translate`, `glm_rotate` are pre functions and are similar to C++ **glm** which you are familiar with. 
 
 In new versions of **cglm** we added `glm_translated`, `glm_rotated`... which are post functions, 
-they are useful in some cases, e.g. append transform to existing transform (apply/append transform as last transfrom T' = T * Tnew).
+they are useful in some cases, e.g. append transform to existing transform (apply/append transform as last transform T' = T * Tnew).
 
 Post functions are named after pre functions with `ed` suffix, e.g. `glm_translate` -> `glm_translated`. So don't mix them up.
 
@@ -24,7 +24,7 @@ a matrix for you. You don't need to pass identity matrix.
 
 But other functions expect you have a matrix and you want to transform them. If
 you didn't have any existing matrix you have to initialize matrix to identity
-before sending to transfrom functions.
+before sending to transform functions.
 
 There are also functions to decompose transform matrix. These functions can't
 decompose matrix after projected.
@@ -35,7 +35,7 @@ Rotation Center
 Rotating functions uses origin as rotation center (pivot/anchor point),
 since scale factors are stored in rotation matrix, same may also true for scalling.
 cglm provides some functions for rotating around at given point e.g.
-**glm_rotate_at**, **glm_quat_rotate_at**. Use them or follow next section for algorihm ("Rotate or Scale around specific Point (Pivot Point / Anchor Point)").
+**glm_rotate_at**, **glm_quat_rotate_at**. Use them or follow next section for algorithm ("Rotate or Scale around specific Point (Pivot Point / Anchor Point)").
 
 Also **cglm** provides :c:func:`glm_spin` and :c:func:`glm_spinned` functions to rotate around itself. No need to give pivot.
 These functions are useful for rotating around center of object. 
@@ -43,7 +43,7 @@ These functions are useful for rotating around center of object.
 Rotate or Scale around specific Point (Anchor Point)
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
-If you want to rotate model around arbibtrary point follow these steps:
+If you want to rotate model around arbitrary point follow these steps:
 
 1. Move model from pivot point to origin: **translate(-pivot.x, -pivot.y, -pivot.z)**
 2. Apply rotation (or scaling maybe)
@@ -82,11 +82,11 @@ helpers functions works like this (cglm provides reverse order as `ed` suffix e.
 .. code-block:: c
   :linenos:
 
-  TransformMatrix = TransformMatrix * TraslateMatrix; // glm_translate()
+  TransformMatrix = TransformMatrix * TranslateMatrix; // glm_translate()
   TransformMatrix = TransformMatrix * RotateMatrix;   // glm_rotate(), glm_quat_rotate()
   TransformMatrix = TransformMatrix * ScaleMatrix;    // glm_scale()
 
-As you can see it is multipled as right matrix. For instance what will happen if you call `glm_translate` twice?
+As you can see it is multiplied as right matrix. For instance what will happen if you call `glm_translate` twice?
 
 .. code-block:: c
   :linenos:

docs/source/affine2d.rst

@@ -14,7 +14,7 @@ a matrix for you. You don't need to pass identity matrix.
 
 But other functions expect you have a matrix and you want to transform them. If
 you didn't have any existing matrix you have to initialize matrix to identity
-before sending to transfrom functions.
+before sending to transform functions.
 
 Transforms Order
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
@@ -45,7 +45,7 @@ Functions:
     translate existing 2d transform matrix by *v* vector and stores result in same matrix
 
     Parameters:
-      | *[in, out]*  **m**    2d affine transfrom
+      | *[in, out]*  **m**    2d affine transform
       | *[in]*       **v**    translate vector [x, y]
 
 .. c:function:: void glm_translate2d_to(mat3 m, vec2 v, mat3 dest)
@@ -53,7 +53,7 @@ Functions:
     translate existing 2d transform matrix by *v* vector and store result in dest
 
     Parameters:
-      | *[in]*  **m**    2d affine transfrom
+      | *[in]*  **m**    2d affine transform
       | *[in]*  **v**    translate vector [x, y]
       | *[out]* **dest** translated matrix
 
@@ -62,7 +62,7 @@ Functions:
     translate existing 2d transform matrix by x factor
 
     Parameters:
-      | *[in, out]*  **m**    2d affine transfrom
+      | *[in, out]*  **m**    2d affine transform
       | *[in]*       **x**    x factor
 
 .. c:function:: void glm_translate2d_y(mat3 m, float y)
@@ -70,7 +70,7 @@ Functions:
     translate existing 2d transform matrix by y factor
 
     Parameters:
-      | *[in, out]*  **m**    2d affine transfrom
+      | *[in, out]*  **m**    2d affine transform
       | *[in]*       **y**    y factor
 
 .. c:function:: void glm_translate2d_make(mat3 m, vec2 v)
@@ -78,7 +78,7 @@ Functions:
     creates NEW translate 2d transform matrix by *v* vector
 
     Parameters:
-      | *[in, out]*  **m**    affine transfrom
+      | *[in, out]*  **m**    affine transform
       | *[in]*       **v**    translate vector [x, y]
 
 .. c:function:: void glm_scale2d_to(mat3 m, vec2 v, mat3 dest)
@@ -86,7 +86,7 @@ Functions:
     scale existing 2d transform matrix by *v* vector and store result in dest
 
     Parameters:
-      | *[in]*  **m**    affine transfrom
+      | *[in]*  **m**    affine transform
       | *[in]*  **v**    scale vector [x, y]
       | *[out]* **dest** scaled matrix
 
@@ -95,7 +95,7 @@ Functions:
     creates NEW 2d scale matrix by *v* vector
 
     Parameters:
-      | *[in, out]*  **m**    affine transfrom
+      | *[in, out]*  **m**    affine transform
       | *[in]*       **v**    scale vector [x, y]
 
 .. c:function:: void glm_scale2d(mat3 m, vec2 v)
@@ -103,7 +103,7 @@ Functions:
     scales existing 2d transform matrix by *v* vector and stores result in same matrix
 
     Parameters:
-      | *[in, out]*  **m**    affine transfrom
+      | *[in, out]*  **m**    affine transform
       | *[in]*       **v**    translate vector [x, y]
 
 .. c:function:: void glm_scale2d_uni(mat3 m, float s)
@@ -111,7 +111,7 @@ Functions:
     applies uniform scale to existing 2d transform matrix v = [s, s] and stores result in same matrix
 
     Parameters:
-      | *[in, out]*  **m**  affine transfrom
+      | *[in, out]*  **m**  affine transform
       | *[in]*       **s**  scale factor
 
 .. c:function:: void glm_rotate2d_make(mat3 m, float angle)
@@ -119,7 +119,7 @@ Functions:
     creates NEW rotation matrix by angle around *Z* axis
 
     Parameters:
-      | *[in, out]*  **m**      affine transfrom
+      | *[in, out]*  **m**      affine transform
       | *[in]*       **angle**  angle (radians)
 
 .. c:function:: void glm_rotate2d(mat3 m, float angle)
@@ -127,7 +127,7 @@ Functions:
     rotate existing 2d transform matrix around *Z* axis by angle and store result in same matrix
 
     Parameters:
-      | *[in, out]*  **m**      affine transfrom
+      | *[in, out]*  **m**      affine transform
       | *[in]*       **angle**  angle (radians)
 
 .. c:function:: void glm_rotate2d_to(mat3 m, float angle, mat3 dest)
@@ -135,6 +135,6 @@ Functions:
     rotate existing 2d transform matrix around *Z* axis by angle and store result in dest
 
     Parameters:
-      | *[in]*  **m**      affine transfrom
+      | *[in]*  **m**      affine transform
       | *[in]*  **angle**  angle (radians)
       | *[out]* **dest**   rotated matrix

docs/source/api.rst

@@ -1,59 +1,28 @@
-API documentation
+📚 API documentation
 ================================
 
-Some functions may exist twice,
-once for their namespace and once for global namespace
-to make easier to write very common functions
+**cglm** provides a few APIs for similar functions.
 
-For instance, in general we use :code:`glm_vec3_dot` to get dot product
-of two **vec3**. Now we can also do this with :code:`glm_dot`,
-same for *_cross* and so on...
+* 📦 **Inline API**: All functions are inline. You can include **cglm/cglm.h** header
+  to use this API. This is the default API. `glm_` is namespace/prefix for this API. 
+* 📦 **Call API**: All functions are not inline. You need to build *cglm* and link it
+  to your project. You can include **cglm/call.h** header to use this API. `glmc_` is namespace/prefix for this API. 
 
-The original function stays where it is, the function in global namespace
-of same name is just an alias, so there is no call version of those functions.
-e.g there is no func like :code:`glmc_dot` because *glm_dot* is just alias for
-:code:`glm_vec3_dot`
+And also there are also sub categories:
 
-By including **cglm/cglm.h** header you will include all inline version
-of functions. Since functions in this header[s] are inline you don't need to
-build or link *cglm* against your project.
+* 📦 **Array API**: Types are raw arrays and functions takes array as argument. You can include **cglm/cglm.h** header
+  to use this API. This is the default API. `glm_` is namespace/prefix for this API. 
+* 📦 **Struct API**: Types are union/struct and functions takes struct as argument and return structs if needed. You can include **cglm/struct.h** header
+  to use this API. This also includes **cglm/cglm.h** header.`glms_` is namespace/prefix for this API but your can omit or change it, see struct api docs.
+* 📦 **SIMD API**: SIMD functions and helpers. `glmm_` is namespace/prefix for this API. 
 
-But by including **cglm/call.h** header you will include all *non-inline*
-version of functions. You need to build *cglm* and link it.
-Follow the :doc:`build` documentation for this
+📌 Since struct api and call api are built top of inline array api, follow inline array api docs for individual functions. 
 
 .. toctree::
    :maxdepth: 1
-   :caption: API categories:
+   :caption: API documentations:
 
-   affine
-   affine-mat
-   affine2d
-   cam
-   frustum
-   box
-   quat
-   euler
-   mat2
-   mat3
-   mat4
-   vec2
-   vec2-ext
-   vec3
-   vec3-ext
-   vec4
-   vec4-ext
-   ivec2
-   ivec3
-   ivec4
-   color
-   plane
-   project
-   util
-   io
-   call
-   sphere
-   curve
-   bezier
-   version
-   ray
+   api_inline_array
+   api_struct
+   api_call
+   api_simd

docs/source/api_call.rst

@@ -0,0 +1,11 @@
+Call API
+================================
+
+Call API is pre-built API for making calls from library. It is built on top of the array api. **glmc_** is the namespace for the call api. 
+**c** stands for call. 
+
+You need to built cglm to use call api. See build instructions (:doc:`build`) for more details.
+
+The functions except namespace **glmc_** are same as inline api. See ( :doc:`api_inline_array` ) for more details.
+
+📌 In the future we can define option to forward inline functions or struct api to call api.

docs/source/api_inline_array.rst

@@ -0,0 +1,77 @@
+Array API - Inline (Default)
+========================================
+
+This is the default API of *cglm*. All functions are forced to be inlined 
+and struct api, call api uses this inline api to share implementation.
+
+📌 Call api is also array api but it is not inlined. 
+In the future there may be option to forward struct api to call api instead of inline api to reduce binary size if needed.
+
+📌 **USE this API docs for similar functions in struct and call api** 
+
+📌 In struct api you can omit namespace e.g :code:`glms_vec3_dot` can be called as :code:`vec3_dot` in struct api, see :doc:`api_struct` to configure struct api for more details.
+📌 In struct api functions can return struct/union
+📌 In struct api you can access items like **.x**, **.y**, **.z**, **.w**, **.r**, **.g**, **.b**, **.a**, **.m00**, **m01**... 
+
+Some functions may exist twice, once for their namespace and once for global namespace
+to make easier to write very common functions
+
+For instance, in general we use :code:`glm_vec3_dot` to get dot product
+of two **vec3**. Now we can also do this with :code:`glm_dot`,
+same for *_cross* and so on...
+
+The original function stays where it is, the function in global namespace
+of same name is just an alias, so there is no call version of those functions.
+e.g there is no func like :code:`glmc_dot` because *glm_dot* is just alias for
+:code:`glm_vec3_dot`
+
+By including **cglm/cglm.h** header you will include all inline version
+of functions. Since functions in this header[s] are inline you don't need to
+build or link *cglm* against your project.
+
+But by including **cglm/call.h** header you will include all *non-inline*
+version of functions. You need to build *cglm* and link it.
+Follow the :doc:`build` documentation for this
+
+.. toctree::
+   :maxdepth: 1
+   :caption: API categories:
+
+   affine
+   affine-mat
+   affine2d
+   cam
+   frustum
+   box
+   aabb2d
+   quat
+   euler
+   mat2
+   mat2x3
+   mat2x4
+   mat3
+   mat3x2
+   mat3x4
+   mat4
+   mat4x2
+   mat4x3
+   vec2
+   vec2-ext
+   vec3
+   vec3-ext
+   vec4
+   vec4-ext
+   ivec2
+   ivec3
+   ivec4
+   color
+   plane
+   project
+   util
+   io
+   call
+   sphere
+   curve
+   bezier
+   version
+   ray

docs/source/api_simd.rst

@@ -0,0 +1,12 @@
+SIMD API
+================================
+
+SIMD api is special api for SIMD operations. **glmm_** prefix is used for SIMD operations in cglm. It is used in many places in cglm.
+You can use it for your own SIMD operations too. In the future the api may be extended by time.
+
+Supported SIMD architectures ( may vary by time )
+
+* SSE / SSE2
+* AVX
+* NEON
+* WASM 128

docs/source/api_struct.rst

@@ -0,0 +1,101 @@
+Struct API 
+================================
+
+Struct API is alternative API to array api to use **cglm** with improved type safety and easy to use. 
+Since struct api is built top of array api, every struct API is not documented here. 
+See array api documentation for more information for individual functions.
+
+By default struct api adds `s` suffix to every type name e.g. vec3s, mat4s, versors etc. 
+Also struct api `s` suffix to namespace e.g. `glms_vec3_add`, `glms_mat4_mul` etc.
+
+By starting v0.9.0, struct api namespace is configurable. We can omit **glms_** namespace or 
+even change it with custom name to move existing api integrations to **cglm** more easliy...
+We can also add **s** to function names if we want e.g. `glms_vec3_add()` -> `vec3_add()` or `vec3s_add()`.
+
+By including **cglm/struct.h** header you will include all struct api. It will also include **cglm/cglm.h** too. 
+Since struct apis are inline you don't need to build or link *cglm* against 
+your project unless if you want to use pre-built call-api too.
+
+Struct API is built top of array api. So you can mix them. 
+Use **.raw** union member to access raw array data to use it with array api. 
+
+Unlike array api ([0], [1], [0][0] ...), it is also possible to use struct api 
+with **.x**, **.y**, **.z**, **.w**, **.r**, **.g**, **.b**, **.a**, **.m00**, **m01**... 
+accessors to access individual elements/properties of vectors and matrices.
+
+Struct API usage:
+-----------------
+
+.. code-block:: c
+
+   #include <cglm/struct.h>
+   
+   mat4s m1  = glms_mat4_identity(); /* init...  */
+   mat4s m2  = glms_mat4_identity(); /* init...  */
+   mat4s m3  = glms_mat4_mul(glms_mat4_mul(m1, m2), glms_mat4_mul(m3, m4));
+   
+   vec3s v1  = { 1.0f, 0.0f, 0.0f };
+   vec3s v2  = { 0.0f, 1.0f, 0.0f };
+   vec4s v4  = { 0.0f, 1.0f, 0.0f, 0.0f };
+   vec4  v5a = { 0.0f, 1.0f, 0.0f, 0.0f };
+   
+   mat4s m4 = glms_rotate(m3, M_PI_2, 
+                          glms_vec3_cross(glms_vec3_add(v1, v6) 
+                                          glms_vec3_add(v1, v7)));
+
+   v1.x = 1.0f; v1.y = 0.0f; v1.z = 0.0f;
+   // or
+   v1.raw[0] = 1.0f; v1.raw[1] = 0.0f; v1.raw[2] = 0.0f;
+
+   /* use struct api with array api (mix them). */
+   /* use .raw to access array and use it with array api */
+
+   glm_vec4_add(m4.col[0].raw, v5a, m4.col[0].raw);
+   glm_mat4_mulv(m4.raw, v4.raw, v5a);
+
+or omit `glms_` namespace completely (see options below):
+
+.. code-block:: c
+
+   #define CGLM_OMIT_NS_FROM_STRUCT_API
+
+   #include <cglm/struct.h>
+
+   mat4s m1  = mat4_identity(); /* init...  */
+   mat4s m2  = mat4_identity(); /* init...  */
+   mat4s m3  = mat4_mul(mat4_mul(m1, m2), mat4_mul(m3, m4));
+   
+   vec3s v1  = { 1.0f, 0.0f, 0.0f };
+   vec3s v2  = { 0.0f, 1.0f, 0.0f };
+   vec4s v4  = { 0.0f, 1.0f, 0.0f, 0.0f };
+   vec4  v5a = { 0.0f, 1.0f, 0.0f, 0.0f };
+   
+   mat4s m4 = glms_rotate(m3, M_PI_2, 
+                          vec3_cross(vec3_add(v1, v6) 
+                                     vec3_add(v1, v7)));
+   
+   v1.x = 1.0f; v1.y = 0.0f; v1.z = 0.0f;
+   // or
+   v1.raw[0] = 1.0f; v1.raw[1] = 0.0f; v1.raw[2] = 0.0f;
+
+   /* use struct api with array api (mix them) */
+   glm_vec4_add(m4.col[0].raw, v5a, m4.col[0].raw);
+   glm_mat4_mulv(m4.raw, v4.raw, v5a);
+   
+Configuring the Struct API:
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+To configure the Struct API namespace, you can define the following macros before including the cglm/struct.h header:
+
+- **CGLM_OMIT_NS_FROM_STRUCT_API**: omits CGLM_STRUCT_API_NS (`glms_`) namespace completely if there is sub namespace e.g `mat4_`, `vec4_` ... DEFAULT is not defined
+- **CGLM_STRUCT_API_NS**:           define name space for struct api, DEFAULT is **glms**
+- **CGLM_STRUCT_API_NAME_SUFFIX**:  define name suffix, DEFAULT is **empty** e.g defining it as #define CGLM_STRUCT_API_NAME_SUFFIX  s will add s suffix to mat4_mul -> mat4s_mul
+
+❗️ IMPORTANT ❗️
+
+It's a good idea to set up your config macros in build settings like CMake, Xcode, or Visual Studio. This is especially important if you're using features like Modules in Xcode, where adding macros directly before the **cglm** headers might not work.
+
+Detailed documentation for Struct API:
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+Since struct api if built top of array api, see array api functions for more information about individual functions.

docs/source/box.rst

@@ -82,7 +82,7 @@ Functions documentation
     Parameters:
       | *[in]*  **box**      bounding box
       | *[in]*  **cropBox**  crop box
-      | *[in]*  **clampBox** miniumum box
+      | *[in]*  **clampBox** minimum box
       | *[out]* **dest**     cropped bounding box
 
 .. c:function:: bool  glm_aabb_frustum(vec3 box[2], vec4 planes[6])

docs/source/build.rst

@@ -3,9 +3,9 @@ Build cglm
 
 | **cglm** does not have any external dependencies.
 
-**NOTE:**
-If you only need to inline versions, you don't need to build **cglm**, you don't need to link it to your program.
-Just import cglm to your project as dependency / external lib by copy-paste then use it as usual
+.. note::
+   If you only need to inline versions, you don't need to build **cglm**, you don't need to link it to your program.
+   Just import cglm to your project as dependency / external lib by copy-paste then use it as usual
 
 CMake (All platforms):
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

docs/source/call.rst

@@ -3,17 +3,17 @@
 precompiled functions (call)
 ================================================================================
 
-All funcitons in **glm_** namespace are forced to **inline**.
+All functions in **glm_** namespace are forced to **inline**.
 Most functions also have pre-compiled version.
 
 Precompiled versions are in **glmc_** namespace. *c* in the namespace stands for
 "call".
 
-Since precompiled functions are just wrapper for inline verisons,
+Since precompiled functions are just wrapper for inline versions,
 these functions are not documented individually.
 It would be duplicate documentation also it
-would be hard to sync documentation between inline and call verison for me.
+would be hard to sync documentation between inline and call version for me.
 
-By including **clgm/cglm.h** you include all inline verisons. To get precompiled
+By including **clgm/cglm.h** you include all inline versions. To get precompiled
 versions you need to include **cglm/call.h** header it also includes all
-call versions plus *clgm/cglm.h* (inline verisons)
+call versions plus *clgm/cglm.h* (inline versions)

docs/source/cam.rst

@@ -18,10 +18,10 @@ fast if you don't care specific projection values.
 *_decomp* means decompose; these function can help to decompose projection
 matrices.
 
- **NOTE**: Be careful when working with high range (very small near, very large
- far) projection matrices. You may not get exact value you gave.
- **float** type cannot store very high precision so you will lose precision.
- Also your projection matrix will be inaccurate due to losing precision
+.. note:: Be careful when working with high range (very small near, very large
+   far) projection matrices. You may not get exact value you gave.
+   **float** type cannot store very high precision so you will lose precision.
+   Also your projection matrix will be inaccurate due to losing precision
 
 Table of contents (click to go):
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
@@ -101,7 +101,7 @@ Functions documentation
     Parameters:
       | *[in]*  **box**      AABB
       | *[in]*  **padding**  padding
-      | *[out]* **d**        result matrix
+      | *[out]* **dest**        result matrix
 
 .. c:function:: void  glm_ortho_aabb_pz(vec3 box[2], float padding, mat4 dest)
 
@@ -113,7 +113,7 @@ Functions documentation
     Parameters:
       | *[in]*  **box**      AABB
       | *[in]*  **padding**  padding for near and far
-      | *[out]* **d**        result matrix
+      | *[out]* **dest**        result matrix
 
     Returns:
       square of norm / magnitude
@@ -178,7 +178,7 @@ Functions documentation
 
     | set up view matrix
 
-    **NOTE:** The UP vector must not be parallel to the line of sight from the eye point to the reference point.
+    .. note:: The UP vector must not be parallel to the line of sight from the eye point to the reference point.
 
     Parameters:
       | *[in]*  **eye**     eye vector
@@ -194,7 +194,7 @@ Functions documentation
     target self then this might be useful. Because you need to get target
     from direction.
 
-    **NOTE:** The UP vector must not be parallel to the line of sight from the eye point to the reference point.
+    .. note:: The UP vector must not be parallel to the line of sight from the eye point to the reference point.
 
     Parameters:
       | *[in]*  **eye**     eye vector
@@ -250,7 +250,7 @@ Functions documentation
 .. c:function:: void  glm_persp_decomp_y(mat4 proj, float *top, float *bottom)
 
     | decomposes top and bottom values of perspective projection.
-    | y stands for y axis (top / botom axis)
+    | y stands for y axis (top / bottom axis)
 
     Parameters:
       | *[in]*   **proj**    perspective projection matrix

docs/source/conf.py

@@ -62,16 +62,16 @@ author = u'Recep Aslantas'
 # built documents.
 #
 # The short X.Y version.
-version = u'0.8.8'
+version = u'0.9.4'
 # The full version, including alpha/beta/rc tags.
-release = u'0.8.8'
+release = u'0.9.4'
 
 # The language for content autogenerated by Sphinx. Refer to documentation
 # for a list of supported languages.
 #
 # This is also used if you do content translation via gettext catalogs.
 # Usually you set "language" from the command line for these cases.
-language = None
+language = 'en'
 
 # List of patterns, relative to source directory, that match files and
 # directories to ignore when looking for source files.
@@ -91,6 +91,7 @@ todo_include_todos = False
 # a list of builtin themes.
 #
 html_theme = 'sphinx_rtd_theme'
+pygments_style = 'monokai'
 
 # Theme options are theme-specific and customize the look and feel of a theme
 # further.  For a list of options available for each theme, see the
@@ -111,8 +112,11 @@ html_theme_options = {
 # Add any paths that contain custom static files (such as style sheets) here,
 # relative to this directory. They are copied after the builtin static files,
 # so a file named "default.css" will overwrite the builtin "default.css".
-html_static_path = ['_static']
+html_static_path = ['sphinx-static']
 
+# Add customm CSS and JS files
+html_css_files = ['theme_overrides.css']
+html_js_files = []
 
 # -- Options for HTMLHelp output ------------------------------------------
 

docs/source/features.rst

@@ -1,8 +1,8 @@
 Features
 ================================================================================
 
-* **scalar** and **simd** (sse, avx, neon...) optimizations
-* option to use different clipspaces e.g. Left Handed, Zero-to-One... (currrently right handed negative-one is default)
+* **scalar** and **simd** (sse, avx, neon, wasm...) optimizations
+* option to use different clipspaces e.g. Left Handed, Zero-to-One... (currently right handed negative-one is default)
 * array api and struct api, you can use arrays or structs.
 * general purpose matrix operations (mat4, mat3)
 * chain matrix multiplication (square only)
@@ -18,6 +18,7 @@ Features
 * inline or pre-compiled function call
 * frustum (extract view frustum planes, corners...)
 * bounding box (AABB in Frustum (culling), crop, merge...)
+* 2d bounding box (crop, merge...)
 * bounding sphere
 * project, unproject
 * easing functions

docs/source/frustum.rst

@@ -11,9 +11,9 @@ not **vec3**. If you want to store them to save space you msut convert them
 yourself.
 
 **vec4** is used to speed up functions need to corners. This is why frustum
-fucntions use *vec4* instead of *vec3*
+functions use *vec4* instead of *vec3*
 
-Currenty related-functions use [-1, 1] clip space configuration to extract
+Currently related-functions use [-1, 1] clip space configuration to extract
 corners but you can override it by prodiving **GLM_CUSTOM_CLIPSPACE** macro.
 If you provide it then you have to all bottom macros as *vec4*
 

docs/source/getting_started.rst

@@ -22,7 +22,7 @@ Types:
   typedef CGLM_ALIGN_IF(16) vec4  mat4[4];
   #endif
 
-As you can see types don't store extra informations in favor of space.
+As you can see types don't store extra information in favor of space.
 You can send these values e.g. matrix to OpenGL directly without casting or calling a function like *value_ptr*
 
 Alignment Is Required:
@@ -35,16 +35,16 @@ Alignment Is Required:
 
   | Check :doc:`opt` page for more details
 
-  Also alignment is disabled for older msvc verisons as default. Now alignment is only required in Visual Studio 2017 version 15.6+ if CGLM_ALL_UNALIGNED macro is not defined.
+  Also alignment is disabled for older msvc versions as default. Now alignment is only required in Visual Studio 2017 version 15.6+ if CGLM_ALL_UNALIGNED macro is not defined.
 
 Allocations:
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 *cglm* doesn't alloc any memory on heap. So it doesn't provide any allocator.
 You must allocate memory yourself. You should alloc memory for out parameters too if you pass pointer of memory location. When allocating memory, don't forget that **vec4** and **mat4** require alignment.
 
-**NOTE:** Unaligned **vec4** and unaligned **mat4** operations will be supported in the future. Check todo list.
-Because you may want to multiply a CGLM matrix with external matrix.
-There is no guarantee that non-CGLM matrix is aligned. Unaligned types will have *u* prefix e.g. **umat4**
+.. note:: Unaligned **vec4** and unaligned **mat4** operations will be supported in the future. Check todo list.
+   Because you may want to multiply a CGLM matrix with external matrix.
+   There is no guarantee that non-CGLM matrix is aligned. Unaligned types will have *u* prefix e.g. **umat4**
 
 Array vs Struct:
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

docs/source/index.rst

@@ -28,7 +28,7 @@ considered to be supported as optional.
    opengl
 
 .. toctree::
-   :maxdepth: 2
+   :maxdepth: 3
    :caption: API:
 
    api

docs/source/io.rst

@@ -24,12 +24,12 @@ Example to print mat4 matrix:
    /* ... */
    glm_mat4_print(transform, stderr);
 
-**NOTE:** print functions use **%0.4f** precision if you need more
-(you probably will in some cases), you can change it temporary.
-cglm may provide precision parameter in the future
+.. note:: print functions use **%0.4f** precision if you need more
+   (you probably will in some cases), you can change it temporary.
+   cglm may provide precision parameter in the future.
 
 Changes since **v0.7.3**:
-* Now mis-alignment of columns are fixed: larger numbers are printed via %g and others are printed via %f. Column withs are calculated before print.
+* Now mis-alignment of columns are fixed: larger numbers are printed via %g and others are printed via %f. Column widths are calculated before print.
 * Now values are colorful ;)
 * Some print improvements
 * New options with default values:
@@ -53,8 +53,11 @@ Functions:
 1. :c:func:`glm_mat4_print`
 #. :c:func:`glm_mat3_print`
 #. :c:func:`glm_vec4_print`
+#. :c:func:`glm_ivec4_print`
 #. :c:func:`glm_vec3_print`
 #. :c:func:`glm_ivec3_print`
+#. :c:func:`glm_vec2_print`
+#. :c:func:`glm_ivec2_print`
 #. :c:func:`glm_versor_print`
 #. :c:func:`glm_aabb_print`
 
@@ -63,7 +66,7 @@ Functions documentation
 
 .. c:function:: void  glm_mat4_print(mat4 matrix, FILE * __restrict ostream)
 
-    | print mat4 to given stream
+    | print matrix to given stream
 
     Parameters:
       | *[in]*  **matrix**   matrix
@@ -71,7 +74,7 @@ Functions documentation
 
 .. c:function:: void  glm_mat3_print(mat3 matrix, FILE * __restrict ostream)
 
-    | print mat3 to given stream
+    | print matrix to given stream
 
     Parameters:
       | *[in]*  **matrix**   matrix
@@ -79,7 +82,15 @@ Functions documentation
 
 .. c:function:: void  glm_vec4_print(vec4 vec, FILE * __restrict ostream)
 
-    | print vec4 to given stream
+    | print vector to given stream
+
+    Parameters:
+      | *[in]*  **vec**      vector
+      | *[in]*  **ostream**  FILE to write
+
+.. c:function:: void  glm_ivec4_print(ivec4 vec, FILE * __restrict ostream)
+
+    | print vector to given stream
 
     Parameters:
       | *[in]*  **vec**      vector
@@ -87,7 +98,7 @@ Functions documentation
 
 .. c:function:: void  glm_vec3_print(vec3 vec, FILE * __restrict ostream)
 
-    | print vec3 to given stream
+    | print vector to given stream
 
     Parameters:
       | *[in]*  **vec**      vector
@@ -95,12 +106,29 @@ Functions documentation
 
 .. c:function:: void  glm_ivec3_print(ivec3 vec, FILE * __restrict ostream)
 
-    | print ivec3 to given stream
+    | print vector to given stream
 
     Parameters:
       | *[in]*  **vec**      vector
       | *[in]*  **ostream**  FILE to write
 
+.. c:function:: void  glm_vec2_print(vec2 vec, FILE * __restrict ostream)
+
+    | print vector to given stream
+
+    Parameters:
+      | *[in]*  **vec**      vector
+      | *[in]*  **ostream**  FILE to write
+
+.. c:function:: void  glm_ivec2_print(ivec2 vec, FILE * __restrict ostream)
+
+    | print vector to given stream
+
+    Parameters:
+      | *[in]*  **vec**      vector
+      | *[in]*  **ostream**  FILE to write
+
+
 .. c:function:: void  glm_versor_print(versor vec, FILE * __restrict ostream)
 
     | print quaternion to given stream
@@ -115,5 +143,5 @@ Functions documentation
 
     Parameters:
       | *[in]*  **vec**      aabb (axis-aligned bounding box)
-      | *[in]*  **tag**      tag to find it more easly in logs
+      | *[in]*  **tag**      tag to find it more easily in logs
       | *[in]*  **ostream**  FILE to write

docs/source/ivec2.rst

@@ -8,18 +8,30 @@ Header: cglm/ivec2.h
 Table of contents (click to go):
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
+Macros:
+
+1. GLM_IVEC2_ONE_INIT
+#. GLM_IVEC2_ZERO_INIT
+#. GLM_IVEC2_ONE
+#. GLM_IVEC2_ZERO
+
 Functions:
 
 1. :c:func:`glm_ivec2`
 #. :c:func:`glm_ivec2_copy`
 #. :c:func:`glm_ivec2_zero`
 #. :c:func:`glm_ivec2_one`
+#. :c:func:`glm_ivec2_dot`
+#. :c:func:`glm_ivec2_cross`
 #. :c:func:`glm_ivec2_add`
 #. :c:func:`glm_ivec2_adds`
 #. :c:func:`glm_ivec2_sub`
 #. :c:func:`glm_ivec2_subs`
 #. :c:func:`glm_ivec2_mul`
 #. :c:func:`glm_ivec2_scale`
+#. :c:func:`glm_ivec2_div`
+#. :c:func:`glm_ivec2_divs`
+#. :c:func:`glm_ivec2_mod`
 #. :c:func:`glm_ivec2_distance2`
 #. :c:func:`glm_ivec2_distance`
 #. :c:func:`glm_ivec2_maxv`
@@ -60,6 +72,30 @@ Functions documentation
     Parameters:
       | *[out]* **v** vector
 
+.. c:function:: int glm_ivec2_dot(ivec2 a, ivec2 b)
+
+    dot product of ivec2
+
+    Parameters:
+      | *[in]*  **a**  vector1
+      | *[in]*  **b**  vector2
+
+    Returns:
+      dot product
+
+.. c:function:: int glm_ivec2_cross(ivec2 a, ivec2 b)
+
+    cross product of two vector (RH)
+
+    | ref: http://allenchou.net/2013/07/cross-product-of-2d-vectors/
+
+    Parameters:
+      | *[in]*  **a**     vector 1
+      | *[in]*  **b**     vector 2
+
+    Returns:
+      Z component of cross product
+
 .. c:function:: void glm_ivec2_add(ivec2 a, ivec2 b, ivec2 dest)
 
     add vector [a] to vector [b] and store result in [dest]
@@ -114,6 +150,33 @@ Functions documentation
       | *[in]*  **s**    scalar
       | *[out]* **dest** destination
 
+.. c:function:: void glm_ivec2_div(ivec2 a, ivec2 b, ivec2 dest)
+
+    div vector with another component-wise division: d = a / b
+
+    Parameters:
+      | *[in]*  **a**     vector 1
+      | *[in]*  **b**     vector 2
+      | *[out]* **dest**  result = (a[0] / b[0], a[1] / b[1], a[2] / b[2])
+
+.. c:function:: void glm_ivec2_divs(ivec2 v, int s, ivec2 dest)
+
+    div vector with scalar: d = v / s
+
+    Parameters:
+      | *[in]*  **v**     vector
+      | *[in]*  **s**     scalar
+      | *[out]* **dest**  result = (a[0] / s, a[1] / s, a[2] / s)
+
+.. c:function:: void glm_ivec2_mod(ivec2 a, ivec2 b, ivec2 dest)
+
+    mod vector with another component-wise modulo: d = a % b
+
+    Parameters:
+      | *[in]*  **a**     vector
+      | *[in]*  **b**     scalar
+      | *[out]* **dest**  result = (a[0] % b[0], a[1] % b[1])
+
 .. c:function:: int glm_ivec2_distance2(ivec2 a, ivec2 b)
 
     squared distance between two vectors
@@ -136,6 +199,31 @@ Functions documentation
     Returns:
         distance
 
+.. c:function:: void  glm_ivec2_fill(ivec2 v, int val)
+
+    fill a vector with specified value
+
+    Parameters:
+      | *[out]*  **v**    vector
+      | *[in]*   **val**  value
+
+.. c:function:: bool  glm_ivec2_eq(ivec2 v, int val)
+
+    check if vector is equal to value
+
+    Parameters:
+      | *[in]*  **v**    vector
+      | *[in]*  **val**  value
+
+.. c:function:: bool  glm_ivec2_eqv(ivec2 v1, ivec2 v2)
+
+    check if vector is equal to another vector
+
+    Parameters:
+      | *[in]*  **vec**   vector 1
+      | *[in]*  **vec**   vector 2
+
+
 .. c:function:: void glm_ivec2_maxv(ivec2 a, ivec2 b, ivec2 dest)
 
     set each member of dest to greater of vector a and b

docs/source/ivec3.rst

@@ -8,20 +8,36 @@ Header: cglm/ivec3.h
 Table of contents (click to go):
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
+Macros:
+
+1. GLM_IVEC3_ONE_INIT
+#. GLM_IVEC3_ZERO_INIT
+#. GLM_IVEC3_ONE
+#. GLM_IVEC3_ZERO
+
 Functions:
 
 1. :c:func:`glm_ivec3`
 #. :c:func:`glm_ivec3_copy`
 #. :c:func:`glm_ivec3_zero`
 #. :c:func:`glm_ivec3_one`
+#. :c:func:`glm_ivec3_dot`
+#. :c:func:`glm_ivec3_norm2`
+#. :c:func:`glm_ivec3_norm`
 #. :c:func:`glm_ivec3_add`
 #. :c:func:`glm_ivec3_adds`
 #. :c:func:`glm_ivec3_sub`
 #. :c:func:`glm_ivec3_subs`
 #. :c:func:`glm_ivec3_mul`
 #. :c:func:`glm_ivec3_scale`
+#. :c:func:`glm_ivec3_div`
+#. :c:func:`glm_ivec3_divs`
+#. :c:func:`glm_ivec3_mod`
 #. :c:func:`glm_ivec3_distance2`
 #. :c:func:`glm_ivec3_distance`
+#. :c:func:`glm_ivec3_fill`
+#. :c:func:`glm_ivec3_eq`
+#. :c:func:`glm_ivec3_eqv`
 #. :c:func:`glm_ivec3_maxv`
 #. :c:func:`glm_ivec3_minv`
 #. :c:func:`glm_ivec3_clamp`
@@ -60,6 +76,39 @@ Functions documentation
     Parameters:
       | *[out]* **v** vector
 
+.. c:function:: int glm_ivec3_dot(ivec3 a, ivec3 b)
+
+    dot product of ivec3
+
+    Parameters:
+      | *[in]*  **a**  vector1
+      | *[in]*  **b**  vector2
+
+    Returns:
+      dot product
+
+.. c:function:: int glm_ivec3_norm2(ivec3 v)
+
+    norm * norm (magnitude) of vector
+
+    we can use this func instead of calling norm * norm, because it would call
+    sqrtf function twice but with this func we can avoid func call, maybe this is
+    not good name for this func
+
+    Parameters:
+      | *[in]*  **v**   vector
+
+    Returns:
+      square of norm / magnitude, cast to an integer
+
+.. c:function:: int glm_ivec3_norm(ivec3 vec)
+
+    | euclidean norm (magnitude), also called L2 norm
+    | this will give magnitude of vector in euclidean space
+
+    Parameters:
+      | *[in]*  **vec**   vector
+
 .. c:function:: void glm_ivec3_add(ivec3 a, ivec3 b, ivec3 dest)
 
     add vector [a] to vector [b] and store result in [dest]
@@ -114,6 +163,33 @@ Functions documentation
       | *[in]*  **s**    scalar
       | *[out]* **dest** destination
 
+.. c:function:: void glm_ivec3_div(ivec3 a, ivec3 b, ivec3 dest)
+
+    div vector with another component-wise division: d = a / b
+
+    Parameters:
+      | *[in]*  **a**     vector 1
+      | *[in]*  **b**     vector 2
+      | *[out]* **dest**  result = (a[0] / b[0], a[1] / b[1], a[2] / b[2])
+
+.. c:function:: void glm_ivec3_divs(ivec3 v, int s, ivec3 dest)
+
+    div vector with scalar: d = v / s
+
+    Parameters:
+      | *[in]*  **v**     vector
+      | *[in]*  **s**     scalar
+      | *[out]* **dest**  result = (a[0] / s, a[1] / s, a[2] / s)
+
+.. c:function:: void glm_ivec3_mod(ivec3 a, ivec3 b, ivec3 dest)
+
+    mod vector with another component-wise modulo: d = a % b
+
+    Parameters:
+      | *[in]*  **a**     vector
+      | *[in]*  **b**     scalar
+      | *[out]* **dest**  result = (a[0] % b[0], a[1] % b[1], a[2] % b[2])
+
 .. c:function:: int glm_ivec3_distance2(ivec3 a, ivec3 b)
 
     squared distance between two vectors
@@ -136,6 +212,30 @@ Functions documentation
     Returns:
         distance
 
+.. c:function:: void  glm_ivec3_fill(ivec3 v, int val)
+
+    fill a vector with specified value
+
+    Parameters:
+      | *[out]*  **v**    vector
+      | *[in]*   **val**  value
+
+.. c:function:: bool  glm_ivec3_eq(ivec3 v, int val)
+
+    check if vector is equal to value
+
+    Parameters:
+      | *[in]*  **v**    vector
+      | *[in]*  **val**  value
+
+.. c:function:: bool  glm_ivec3_eqv(ivec3 v1, ivec3 v2)
+
+    check if vector is equal to another vector
+
+    Parameters:
+      | *[in]*  **vec**   vector 1
+      | *[in]*  **vec**   vector 2
+
 .. c:function:: void glm_ivec3_maxv(ivec3 a, ivec3 b, ivec3 dest)
 
     set each member of dest to greater of vector a and b

docs/source/ivec4.rst

@@ -8,6 +8,13 @@ Header: cglm/ivec4.h
 Table of contents (click to go):
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
+Macros:
+
+1. GLM_IVEC4_ONE_INIT
+#. GLM_IVEC4_ZERO_INIT
+#. GLM_IVEC4_ONE
+#. GLM_IVEC4_ZERO
+
 Functions:
 
 1. :c:func:`glm_ivec4`

docs/source/mat2.rst

@@ -32,6 +32,7 @@ Functions:
 #. :c:func:`glm_mat2_swap_col`
 #. :c:func:`glm_mat2_swap_row`
 #. :c:func:`glm_mat2_rmc`
+#. :c:func:`glm_mat2_make`
 
 Functions documentation
 ~~~~~~~~~~~~~~~~~~~~~~~
@@ -46,7 +47,7 @@ Functions documentation
 
 .. c:function:: void glm_mat2_identity(mat2 mat)
 
-    copy identity mat2 to mat, or makes mat to identiy
+    copy identity mat2 to mat, or makes mat to identity
 
     Parameters:
       | *[out]* **mat**  matrix
@@ -64,11 +65,12 @@ Functions documentation
     make given matrix zero
 
     Parameters:
-      | *[in,out]* **mat**  matrix to
+      | *[in,out]* **mat**  matrix
 
 .. c:function:: void glm_mat2_mul(mat2 m1, mat2 m2, mat2 dest)
 
     multiply m1 and m2 to dest
+
     m1, m2 and dest matrices can be same matrix, it is possible to write this:
 
     .. code-block:: c
@@ -92,7 +94,7 @@ Functions documentation
 
 .. c:function:: void glm_mat2_transpose(mat2 m)
 
-    tranpose mat2 and store result in same matrix
+    transpose mat2 and store result in same matrix
 
     Parameters:
       | *[in]*  **mat**   source
@@ -100,7 +102,7 @@ Functions documentation
 
 .. c:function:: void glm_mat2_mulv(mat2 m, vec2 v, vec2 dest)
 
-    multiply mat4 with vec4 (column vector) and store in dest vector
+    multiply mat2 with vec2 (column vector) and store in dest vector
 
     Parameters:
       | *[in]*  **mat**   mat2 (left)
@@ -112,8 +114,8 @@ Functions documentation
     multiply matrix with scalar
 
     Parameters:
-      | *[in, out]* **mat**   matrix
-      | *[in]*      **dest**  scalar
+      | *[in, out]* **m**  matrix
+      | *[in]*      **s**  scalar
 
 .. c:function:: float  glm_mat2_det(mat2 mat)
 
@@ -177,3 +179,13 @@ Functions documentation
 
     Returns:
         scalar value e.g. Matrix1x1
+
+.. c:function:: void glm_mat2_make(const float * __restrict src, mat2 dest)
+
+    Create mat2 matrix from pointer
+
+    .. note:: **@src** must contain at least 4 elements.
+
+    Parameters:
+      | *[in]*  **src**  pointer to an array of floats
+      | *[out]* **dest** destination matrix2x2

docs/source/mat2x3.rst

@@ -0,0 +1,140 @@
+.. default-domain:: C
+
+mat2x3
+======
+
+Header: cglm/mat2x3.h
+
+Table of contents (click to go):
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+Macros:
+
+1. GLM_MAT2X3_ZERO_INIT
+#. GLM_MAT2X3_ZERO
+
+Functions:
+
+1. :c:func:`glm_mat2x3_copy`
+#. :c:func:`glm_mat2x3_zero`
+#. :c:func:`glm_mat2x3_make`
+#. :c:func:`glm_mat2x3_mul`
+#. :c:func:`glm_mat2x3_mulv`
+#. :c:func:`glm_mat2x3_transpose`
+#. :c:func:`glm_mat2x3_scale`
+
+Represented
+~~~~~~~~~~~
+
+.. csv-table:: mat2x3
+   :header: "", "column 1", "column 2"
+
+   "row 1", "m00", "m10"
+   "row 2", "m01", "m11"
+   "row 3", "m02", "m12"
+
+Functions documentation
+~~~~~~~~~~~~~~~~~~~~~~~
+
+.. c:function:: void glm_mat2x3_copy(mat2x3 mat, mat2x3 dest)
+
+    copy mat2x3 to another one (dest).
+
+    Parameters:
+      | *[in]*  **mat**   source
+      | *[out]* **dest**  destination
+
+.. c:function:: void glm_mat2x3_zero(mat2x3 mat)
+
+    make given matrix zero
+
+    Parameters:
+      | *[in,out]* **mat**  matrix
+
+.. c:function:: void glm_mat2x3_make(const float * __restrict src, mat2x3 dest)
+
+    Create mat2x3 matrix from pointer
+
+    .. note:: **@src** must contain at least 6 elements.
+
+    Parameters:
+      | *[in]*  **src**  pointer to an array of floats
+      | *[out]* **dest** destination matrix2x3
+
+.. c:function:: void glm_mat2x3_mul(mat2x3 m1, mat3x2 m2, mat3 dest)
+
+    multiply m1 and m2 to dest
+
+    .. code-block:: c
+
+       glm_mat2x3_mul(mat2x3, mat3x2, mat3);
+
+    Parameters:
+      | *[in]*  **m1**    left matrix (mat2x3)
+      | *[in]*  **m2**    right matrix (mat3x2)
+      | *[out]* **dest**  destination matrix (mat3)
+
+    .. csv-table:: mat2x3
+        :header: "", "column 1", "column 2"
+
+        "row 1", "a00", "a10"
+        "row 2", "a01", "a11"
+        "row 3", "a02", "a12"
+
+    .. csv-table:: mat3x2
+        :header: "", "column 1", "column 2", "column 3"
+
+        "row 1", "b00", "b10", "b20"
+        "row 2", "b01", "b11", "b21"
+
+    .. csv-table:: mat3x3
+        :header: "", "column 1", "column 2", "column 3"
+
+        "row 1", "a00 * b00 + a10 * b01", "a00 * b10 + a10 * b11", "a00 * b20 + a10 * b21"
+        "row 2", "a01 * b00 + a11 * b01", "a01 * b10 + a11 * b11", "a01 * b20 + a11 * b21"
+        "row 3", "a02 * b00 + a12 * b01", "a02 * b10 + a12 * b11", "a02 * b20 + a12 * b21"
+
+.. c:function:: void glm_mat2x3_mulv(mat2x3 m, vec2 v, vec3 dest)
+
+    multiply mat2x3 with vec2 (column vector) and store in dest column vector
+
+    Parameters:
+      | *[in]*  **m**     mat2x3 (left)
+      | *[in]*  **v**     vec3 (right, column vector)
+      | *[out]* **dest**  destination (result, column vector)
+
+    .. csv-table:: mat2x3
+        :header: "", "column 1", "column 2"
+
+        "row 1", "m00", "m10"
+        "row 2", "m01", "m11"
+        "row 3", "m02", "m12"
+
+    .. csv-table:: column vec2 (1x2)
+        :header: "", "column 1"
+
+        "row 1", "v0"
+        "row 2", "v1"
+
+    .. csv-table:: column vec3 (1x3)
+        :header: "", "column 1"
+
+        "row 1", "m00 * v0 + m10 * v1"
+        "row 2", "m01 * v0 + m11 * v1"
+        "row 3", "m02 * v0 + m12 * v1"
+
+.. c:function:: void glm_mat2x3_transpose(mat2x3 m, mat3x2 dest)
+
+    transpose matrix and store in dest
+
+    Parameters:
+      | *[in]*  **m**     matrix
+      | *[out]* **dest**  destination
+
+.. c:function:: void  glm_mat2x3_scale(mat2x3 m, float s)
+
+    multiply matrix with scalar
+
+    Parameters:
+      | *[in, out]* **m** matrix
+      | *[in]*      **s** scalar

docs/source/mat2x4.rst

@@ -0,0 +1,145 @@
+.. default-domain:: C
+
+mat2x4
+======
+
+Header: cglm/mat2x4.h
+
+Table of contents (click to go):
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+Macros:
+
+1. GLM_MAT2X4_ZERO_INIT
+#. GLM_MAT2X4_ZERO
+
+Functions:
+
+1. :c:func:`glm_mat2x4_copy`
+#. :c:func:`glm_mat2x4_zero`
+#. :c:func:`glm_mat2x4_make`
+#. :c:func:`glm_mat2x4_mul`
+#. :c:func:`glm_mat2x4_mulv`
+#. :c:func:`glm_mat2x4_transpose`
+#. :c:func:`glm_mat2x4_scale`
+
+Represented
+~~~~~~~~~~~
+
+.. csv-table:: mat2x4
+   :header: "", "column 1", "column 2"
+
+   "row 1", "m00", "m10"
+   "row 2", "m01", "m11"
+   "row 3", "m02", "m12"
+   "row 4", "m03", "m13"
+
+Functions documentation
+~~~~~~~~~~~~~~~~~~~~~~~
+
+.. c:function:: void glm_mat2x4_copy(mat2x4 mat, mat2x4 dest)
+
+    copy mat2x4 to another one (dest).
+
+    Parameters:
+      | *[in]*  **mat**   source
+      | *[out]* **dest**  destination
+
+.. c:function:: void glm_mat2x4_zero(mat2x4 mat)
+
+    make given matrix zero
+
+    Parameters:
+      | *[in,out]* **mat**  matrix
+
+.. c:function:: void glm_mat2x4_make(const float * __restrict src, mat2x4 dest)
+
+    Create mat2x4 matrix from pointer
+
+    .. note:: **@src** must contain at least 8 elements.
+
+    Parameters:
+      | *[in]*  **src**  pointer to an array of floats
+      | *[out]* **dest** destination matrix2x4
+
+.. c:function:: void glm_mat2x4_mul(mat2x4 m1, mat4x2 m2, mat4 dest)
+
+    multiply m1 and m2 to dest
+
+    .. code-block:: c
+
+       glm_mat2x4_mul(mat2x4, mat4x2, mat4);
+
+    Parameters:
+      | *[in]*  **m1**    left matrix (mat2x4)
+      | *[in]*  **m2**    right matrix (mat4x2)
+      | *[out]* **dest**  destination matrix (mat4)
+
+    .. csv-table:: mat2x4
+        :header: "", "column 1", "column 2"
+
+        "row 1", "a00", "a10"
+        "row 2", "a01", "a11"
+        "row 3", "a02", "a12"
+        "row 4", "a03", "a13"
+
+    .. csv-table:: mat4x2
+        :header: "", "column 1", "column 2", "column 3", "column 4"
+
+        "row 1", "b00", "b10", "b20", "b30"
+        "row 2", "b01", "b11", "b21", "b31"
+
+    .. csv-table:: mat4x4
+        :header: "", "column 1", "column 2", "column 3", "column 4"
+
+        "row 1", "a00 * b00 + a10 * b01", "a00 * b10 + a10 * b11", "a00 * b20 + a10 * b21", "a00 * b30 + a10 * b31"
+        "row 2", "a01 * b00 + a11 * b01", "a01 * b10 + a11 * b11", "a01 * b20 + a11 * b21", "a01 * b30 + a11 * b31"
+        "row 3", "a02 * b00 + a12 * b01", "a02 * b10 + a12 * b11", "a02 * b20 + a12 * b21", "a02 * b30 + a12 * b31"
+        "row 4", "a03 * b00 + a13 * b01", "a03 * b10 + a13 * b11", "a03 * b20 + a13 * b21", "a03 * b30 + a13 * b31"
+
+.. c:function:: void glm_mat2x4_mulv(mat2x4 m, vec2 v, vec4 dest)
+
+    multiply mat2x4 with vec2 (column vector) and store in dest column vector
+
+    Parameters:
+      | *[in]*  **m**     mat2x4 (left)
+      | *[in]*  **v**     vec2 (right, column vector)
+      | *[out]* **dest**  destination (result, column vector)
+
+    .. csv-table:: mat2x4
+        :header: "", "column 1", "column 2"
+
+        "row 1", "m00", "m10"
+        "row 2", "m01", "m11"
+        "row 3", "m02", "m12"
+        "row 4", "m03", "m13"
+
+    .. csv-table:: column vec2 (1x2)
+        :header: "", "column 1"
+
+        "row 1", "v0"
+        "row 2", "v1"
+
+    .. csv-table:: column vec4 (1x4)
+        :header: "", "column 1"
+
+        "row 1", "m00 * v0 + m10 * v1"
+        "row 2", "m01 * v0 + m11 * v1"
+        "row 3", "m02 * v0 + m12 * v1"
+        "row 4", "m03 * v0 + m13 * v1"
+
+.. c:function:: void glm_mat2x4_transpose(mat2x4 m, mat4x2 dest)
+
+    transpose matrix and store in dest
+
+    Parameters:
+      | *[in]*  **m**     matrix
+      | *[out]* **dest**  destination
+
+.. c:function:: void  glm_mat2x4_scale(mat2x4 m, float s)
+
+    multiply matrix with scalar
+
+    Parameters:
+      | *[in, out]* **m** matrix
+      | *[in]*      **s** scalar

docs/source/mat3.rst

@@ -34,6 +34,7 @@ Functions:
 #. :c:func:`glm_mat3_swap_col`
 #. :c:func:`glm_mat3_swap_row`
 #. :c:func:`glm_mat3_rmc`
+#. :c:func:`glm_mat3_make`
 
 Functions documentation
 ~~~~~~~~~~~~~~~~~~~~~~~
@@ -48,7 +49,7 @@ Functions documentation
 
 .. c:function:: void  glm_mat3_identity(mat3 mat)
 
-    copy identity mat3 to mat, or makes mat to identiy
+    copy identity mat3 to mat, or makes mat to identity
 
     Parameters:
       | *[out]* **mat**  matrix
@@ -71,6 +72,7 @@ Functions documentation
 .. c:function:: void  glm_mat3_mul(mat3 m1, mat3 m2, mat3 dest)
 
     multiply m1 and m2 to dest
+
     m1, m2 and dest matrices can be same matrix, it is possible to write this:
 
     .. code-block:: c
@@ -94,7 +96,7 @@ Functions documentation
 
 .. c:function:: void  glm_mat3_transpose(mat3 m)
 
-    tranpose mat3 and store result in same matrix
+    transpose mat3 and store result in same matrix
 
     Parameters:
       | *[in]*  **mat**   source
@@ -102,10 +104,10 @@ Functions documentation
 
 .. c:function:: void  glm_mat3_mulv(mat3 m, vec3 v, vec3 dest)
 
-    multiply mat4 with vec4 (column vector) and store in dest vector
+    multiply mat3 with vec3 (column vector) and store in dest vector
 
     Parameters:
-      | *[in]*  **mat**   mat3 (left)
+      | *[in]*  **m**     mat3 (left)
       | *[in]*  **v**     vec3 (right, column vector)
       | *[out]* **dest**  destination (result, column vector)
 
@@ -122,8 +124,8 @@ Functions documentation
     multiply matrix with scalar
 
     Parameters:
-      | *[in, out]* **mat**   matrix
-      | *[in]*      **dest**  scalar
+      | *[in, out]* **m** matrix
+      | *[in]*      **s** scalar
 
 .. c:function:: float  glm_mat3_det(mat3 mat)
 
@@ -187,3 +189,13 @@ Functions documentation
 
     Returns:
         scalar value e.g. Matrix1x1
+
+.. c:function:: void glm_mat3_make(const float * __restrict src, mat3 dest)
+
+    Create mat3 matrix from pointer
+
+    .. note:: **@src** must contain at least 9 elements.
+
+    Parameters:
+      | *[in]*  **src**  pointer to an array of floats
+      | *[out]* **dest** destination matrix3x3

docs/source/mat3x2.rst

@@ -0,0 +1,137 @@
+.. default-domain:: C
+
+mat3x2
+======
+
+Header: cglm/mat3x2.h
+
+Table of contents (click to go):
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+Macros:
+
+1. GLM_MAT3X2_ZERO_INIT
+#. GLM_MAT3X2_ZERO
+
+Functions:
+
+1. :c:func:`glm_mat3x2_copy`
+#. :c:func:`glm_mat3x2_zero`
+#. :c:func:`glm_mat3x2_make`
+#. :c:func:`glm_mat3x2_mul`
+#. :c:func:`glm_mat3x2_mulv`
+#. :c:func:`glm_mat3x2_transpose`
+#. :c:func:`glm_mat3x2_scale`
+
+Represented
+~~~~~~~~~~~
+
+.. csv-table:: mat3x2
+   :header: "", "column 1", "column 2", "column 3"
+
+   "row 1", "m00", "m10", "m20"
+   "row 2", "m01", "m11", "m21"
+
+Functions documentation
+~~~~~~~~~~~~~~~~~~~~~~~
+
+.. c:function:: void glm_mat3x2_copy(mat3x2 mat, mat3x2 dest)
+
+    copy mat3x2 to another one (dest).
+
+    Parameters:
+      | *[in]*  **mat**   source
+      | *[out]* **dest**  destination
+
+.. c:function:: void glm_mat3x2_zero(mat3x2 mat)
+
+    make given matrix zero
+
+    Parameters:
+      | *[in,out]* **mat**  matrix
+
+.. c:function:: void glm_mat3x2_make(const float * __restrict src, mat3x2 dest)
+
+    Create mat3x2 matrix from pointer
+
+    .. note:: **@src** must contain at least 6 elements.
+
+    Parameters:
+      | *[in]*  **src**  pointer to an array of floats
+      | *[out]* **dest** destination matrix3x2
+
+.. c:function:: void glm_mat3x2_mul(mat3x2 m1, mat2x3 m2, mat2 dest)
+
+    multiply m1 and m2 to dest
+
+    .. code-block:: c
+
+       glm_mat3x2_mul(mat3x2, mat2x3, mat2);
+
+    Parameters:
+      | *[in]*  **m1**    left matrix (mat3x2)
+      | *[in]*  **m2**    right matrix (mat2x3)
+      | *[out]* **dest**  destination matrix (mat2)
+
+    .. csv-table:: mat3x2
+        :header: "", "column 1", "column 2", "column 3"
+
+        "row 1", "a00", "a10", "a20"
+        "row 2", "a01", "a11", "a21"
+
+    .. csv-table:: mat2x3
+        :header: "", "column 1", "column 2"
+
+        "row 1", "b00", "b10"
+        "row 2", "b01", "b11"
+        "row 3", "b02", "b12"
+
+    .. csv-table:: mat2x2
+        :header: "", "column 1", "column 2"
+
+        "row 1", "a00 * b00 + a10 * b01 + a20 * b02", "a00 * b10 + a10 * b11 + a20 * b12"
+        "row 2", "a01 * b00 + a11 * b01 + a21 * b02", "a01 * b10 + a11 * b11 + a21 * b12"
+
+.. c:function:: void glm_mat3x2_mulv(mat3x2 m, vec3 v, vec2 dest)
+
+    multiply mat3x2 with vec3 (column vector) and store in dest vector
+
+    Parameters:
+      | *[in]*  **m**     mat3x2 (left)
+      | *[in]*  **v**     vec3 (right, column vector)
+      | *[out]* **dest**  destination (result, column vector)
+
+    .. csv-table:: mat3x2
+        :header: "", "column 1", "column 2", "column 3"
+
+        "row 1", "m00", "m10", "m20"
+        "row 2", "m01", "m11", "m21"
+
+    .. csv-table:: column vec3 (1x3)
+        :header: "", "column 1"
+
+        "row 1", "v0"
+        "row 2", "v1"
+        "row 3", "v2"
+
+    .. csv-table:: column vec2 (1x2)
+        :header: "", "column 1"
+
+        "row 1", "m00 * v0 + m10 * v1 + m20 * v2"
+        "row 2", "m01 * v0 + m11 * v1 + m21 * v2"
+
+.. c:function:: void glm_mat3x2_transpose(mat3x2 m, mat2x3 dest)
+
+    transpose matrix and store in dest
+
+    Parameters:
+      | *[in]*  **m**     matrix
+      | *[out]* **dest**  destination
+
+.. c:function:: void  glm_mat3x2_scale(mat3x2 m, float s)
+
+    multiply matrix with scalar
+
+    Parameters:
+      | *[in, out]* **m** matrix
+      | *[in]*      **s** scalar

docs/source/mat3x4.rst

@@ -0,0 +1,147 @@
+.. default-domain:: C
+
+mat3x4
+======
+
+Header: cglm/mat3x4.h
+
+Table of contents (click to go):
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+Macros:
+
+1. GLM_MAT3X4_ZERO_INIT
+#. GLM_MAT3X4_ZERO
+
+Functions:
+
+1. :c:func:`glm_mat3x4_copy`
+#. :c:func:`glm_mat3x4_zero`
+#. :c:func:`glm_mat3x4_make`
+#. :c:func:`glm_mat3x4_mul`
+#. :c:func:`glm_mat3x4_mulv`
+#. :c:func:`glm_mat3x4_transpose`
+#. :c:func:`glm_mat3x4_scale`
+
+Represented
+~~~~~~~~~~~
+
+.. csv-table:: mat3x4
+   :header: "", "column 1", "column 2", "column 3"
+
+   "row 1", "m00", "m10", "m20"
+   "row 2", "m01", "m11", "m21"
+   "row 3", "m02", "m12", "m22"
+   "row 4", "m03", "m13", "m23"
+
+Functions documentation
+~~~~~~~~~~~~~~~~~~~~~~~
+
+.. c:function:: void glm_mat3x4_copy(mat3x4 mat, mat3x4 dest)
+
+    copy mat3x4 to another one (dest).
+
+    Parameters:
+      | *[in]*  **mat**   source
+      | *[out]* **dest**  destination
+
+.. c:function:: void glm_mat3x4_zero(mat3x4 mat)
+
+    make given matrix zero
+
+    Parameters:
+      | *[in,out]* **mat**  matrix
+
+.. c:function:: void glm_mat3x4_make(const float * __restrict src, mat3x4 dest)
+
+    Create mat3x4 matrix from pointer
+
+    .. note::: **@src** must contain at least 12 elements.
+
+    Parameters:
+      | *[in]*  **src**  pointer to an array of floats
+      | *[out]* **dest** destination matrix3x4
+
+.. c:function:: void glm_mat3x4_mul(mat3x4 m1, mat4x3 m2, mat4 dest)
+
+    multiply m1 and m2 to dest
+
+    .. code-block:: c
+
+       glm_mat3x4_mul(mat3x4, mat4x3, mat4);
+
+    Parameters:
+      | *[in]*  **m1**    left matrix (mat3x4)
+      | *[in]*  **m2**    right matrix (mat4x3)
+      | *[out]* **dest**  destination matrix (mat4)
+
+    .. csv-table:: mat3x4
+        :header: "", "column 1", "column 2", "column 3"
+
+        "row 1", "a00", "a10", "a20"
+        "row 2", "a01", "a11", "a21"
+        "row 3", "a02", "a12", "a22"
+        "row 4", "a03", "a13", "a23"
+
+    .. csv-table:: mat4x3
+        :header: "", "column 1", "column 2", "column 3", "column 4"
+
+        "row 1", "b00", "b10", "b20", "b30"
+        "row 2", "b01", "b11", "b21", "b31"
+        "row 3", "b02", "b12", "b22", "b32"
+
+    .. csv-table:: mat4x4
+        :header: "", "column 1", "column 2", "column 3", "column 4"
+
+        "row 1", "a00 * b00 + a10 * b01 + a20 * b02", "a00 * b10 + a10 * b11 + a20 * b12", "a00 * b20 + a10 * b21 + a20 * b22", "a00 * b30 + a10 * b31 + a20 * b32"
+        "row 2", "a01 * b00 + a11 * b01 + a21 * b02", "a01 * b10 + a11 * b11 + a21 * b12", "a01 * b20 + a11 * b21 + a21 * b22", "a01 * b30 + a11 * b31 + a21 * b32"
+        "row 3", "a02 * b00 + a12 * b01 + a22 * b02", "a02 * b10 + a12 * b11 + a22 * b12", "a02 * b20 + a12 * b21 + a22 * b22", "a02 * b30 + a12 * b31 + a22 * b32"
+        "row 4", "a03 * b00 + a13 * b01 + a23 * b02", "a03 * b10 + a13 * b11 + a23 * b12", "a03 * b20 + a13 * b21 + a23 * b22", "a03 * b30 + a13 * b31 + a23 * b32"
+
+.. c:function:: void glm_mat3x4_mulv(mat3x4 m, vec3 v, vec4 dest)
+
+    multiply mat3x4 with vec3 (column vector) and store in dest vector
+
+    Parameters:
+      | *[in]*  **m**     mat3x4 (left)
+      | *[in]*  **v**     vec3 (right, column vector)
+      | *[out]* **dest**  destination (result, column vector)
+
+    .. csv-table:: mat3x4
+        :header: "", "column 1", "column 2", "column 3"
+
+        "row 1", "m00", "m10", "m20"
+        "row 2", "m01", "m11", "m21"
+        "row 3", "m02", "m12", "m22"
+        "row 4", "m03", "m13", "m23"
+
+    .. csv-table:: column vec3 (1x3)
+        :header: "", "column 1"
+
+        "row 1", "v0"
+        "row 2", "v1"
+        "row 3", "v2"
+
+    .. csv-table:: column vec4 (1x4)
+        :header: "", "column 1"
+
+        "row 1", "m00 * v0 + m10 * v1 + m20 * v2"
+        "row 2", "m01 * v0 + m11 * v1 + m21 * v2"
+        "row 3", "m02 * v0 + m12 * v1 + m22 * v2"
+        "row 4", "m03 * v0 + m13 * v1 + m23 * v2"
+
+.. c:function:: void glm_mat3x4_transpose(mat3x4 m, mat4x3 dest)
+
+    transpose matrix and store in dest
+
+    Parameters:
+      | *[in]*  **m**     matrix
+      | *[out]* **dest**  destination
+
+.. c:function:: void  glm_mat3x4_scale(mat3x4 m, float s)
+
+    multiply matrix with scalar
+
+    Parameters:
+      | *[in, out]* **m** matrix
+      | *[in]*      **s** scalar

docs/source/mat4.rst

@@ -34,8 +34,8 @@ Functions:
 #. :c:func:`glm_mat4_mulN`
 #. :c:func:`glm_mat4_mulv`
 #. :c:func:`glm_mat4_mulv3`
-#. :c:func:`glm_mat3_trace`
-#. :c:func:`glm_mat3_trace3`
+#. :c:func:`glm_mat4_trace`
+#. :c:func:`glm_mat4_trace3`
 #. :c:func:`glm_mat4_quat`
 #. :c:func:`glm_mat4_transpose_to`
 #. :c:func:`glm_mat4_transpose`
@@ -47,6 +47,7 @@ Functions:
 #. :c:func:`glm_mat4_swap_col`
 #. :c:func:`glm_mat4_swap_row`
 #. :c:func:`glm_mat4_rmc`
+#. :c:func:`glm_mat4_make`
 
 Functions documentation
 ~~~~~~~~~~~~~~~~~~~~~~~
@@ -69,7 +70,7 @@ Functions documentation
 
 .. c:function:: void  glm_mat4_identity(mat4 mat)
 
-    copy identity mat4 to mat, or makes mat to identiy
+    copy identity mat4 to mat, or makes mat to identity
 
     Parameters:
       | *[out]* **mat**  matrix
@@ -118,6 +119,7 @@ Functions documentation
 .. c:function:: void  glm_mat4_mul(mat4 m1, mat4 m2, mat4 dest)
 
     multiply m1 and m2 to dest
+
     m1, m2 and dest matrices can be same matrix, it is possible to write this:
 
     .. code-block:: c
@@ -156,7 +158,6 @@ Functions documentation
     Parameters:
       | *[in]*  **m**     mat4 (left)
       | *[in]*  **v**     vec4 (right, column vector)
-      | *[in]*  **last**  4th item to make it vec4
       | *[out]* **dest**  vec4 (result, column vector)
 
 .. c:function:: void  glm_mat4_mulv3(mat4 m, vec3 v, float last, vec3 dest)
@@ -213,7 +214,7 @@ Functions documentation
 
 .. c:function:: void  glm_mat4_transpose(mat4 m)
 
-    tranpose mat4 and store result in same matrix
+    transpose mat4 and store result in same matrix
 
     Parameters:
       | *[in]*  **m**     source
@@ -262,7 +263,7 @@ Functions documentation
     | e.g Newton-Raphson. this should work faster than normal,
     | to get more precise use glm_mat4_inv version.
 
-    | NOTE: You will lose precision, glm_mat4_inv is more accurate
+    .. note:: You will lose precision, glm_mat4_inv is more accurate
 
     Parameters:
       | *[in]*  **mat**   source
@@ -302,3 +303,13 @@ Functions documentation
 
     Returns:
         scalar value e.g. Matrix1x1
+
+.. c:function:: void glm_mat4_make(const float * __restrict src, mat4 dest)
+
+    Create mat4 matrix from pointer
+
+    .. note:: **@src** must contain at least 16 elements.
+
+    Parameters:
+      | *[in]*  **src**  pointer to an array of floats
+      | *[out]* **dest** destination matrix4x4

docs/source/mat4x2.rst

@@ -0,0 +1,141 @@
+.. default-domain:: C
+
+mat4x2
+======
+
+Header: cglm/mat4x2.h
+
+Table of contents (click to go):
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+Macros:
+
+1. GLM_MAT4X2_ZERO_INIT
+#. GLM_MAT4X2_ZERO
+
+Functions:
+
+1. :c:func:`glm_mat4x2_copy`
+#. :c:func:`glm_mat4x2_zero`
+#. :c:func:`glm_mat4x2_make`
+#. :c:func:`glm_mat4x2_mul`
+#. :c:func:`glm_mat4x2_mulv`
+#. :c:func:`glm_mat4x2_transpose`
+#. :c:func:`glm_mat4x2_scale`
+
+Represented
+~~~~~~~~~~~
+
+.. csv-table:: mat4x2
+   :header: "", "column 1", "column 2", "column 3", "column4"
+
+   "row 1", "m00", "m10", "m20", "m30"
+   "row 2", "m01", "m11", "m21", "m31"
+
+Functions documentation
+~~~~~~~~~~~~~~~~~~~~~~~
+
+.. c:function:: void glm_mat4x2_copy(mat4x2 mat, mat4x2 dest)
+
+    copy mat4x2 to another one (dest).
+
+    Parameters:
+      | *[in]*  **mat**   source
+      | *[out]* **dest**  destination
+
+.. c:function:: void glm_mat4x2_zero(mat4x2 mat)
+
+    make given matrix zero
+
+    Parameters:
+      | *[in,out]* **mat**  matrix
+
+.. c:function:: void glm_mat4x2_make(const float * __restrict src, mat4x2 dest)
+
+    Create mat4x2 matrix from pointer
+
+    .. note:: **@src** must contain at least 8 elements.
+
+    Parameters:
+      | *[in]*  **src**  pointer to an array of floats
+      | *[out]* **dest** destination matrix4x2
+
+.. c:function:: void glm_mat4x2_mul(mat4x2 m1, mat2x4 m2, mat2 dest)
+
+    multiply m1 and m2 to dest
+
+    .. code-block:: c
+
+       glm_mat4x2_mul(mat4x2, mat2x4, mat2);
+
+    Parameters:
+      | *[in]*  **m1**    left matrix (mat4x2)
+      | *[in]*  **m2**    right matrix (mat2x4)
+      | *[out]* **dest**  destination matrix (mat2)
+
+    .. csv-table:: mat4x2
+        :header: "", "column 1", "column 2", "column 3", "column 4"
+
+        "row 1", "a00", "a10", "a20", "a30"
+        "row 2", "a01", "a11", "a21", "a31"
+
+    .. csv-table:: mat2x4
+        :header: "", "column 1", "column 2"
+
+        "row 1", "b00", "b10"
+        "row 2", "b01", "b11"
+        "row 3", "b02", "b12"
+        "row 4", "b03", "b13"
+
+    .. csv-table:: mat2x2
+        :header: "", "column 1", "column 2"
+
+        "row 1", "a00 * b00 + a10 * b01 + a20 * b02 + a30 * b03", "a00 * b10 + a10 * b11 + a20 * b12 + a30 * b13"
+        "row 2", "a01 * b00 + a11 * b01 + a21 * b02 + a31 * b03", "a01 * b10 + a11 * b11 + a21 * b12 + a31 * b13"
+
+.. c:function:: void glm_mat4x2_mulv(mat4x2 m, vec4 v, vec2 dest)
+
+    multiply mat4x2 with vec4 (column vector) and store in dest vector
+
+    Parameters:
+      | *[in]*  **m**     mat4x2 (left)
+      | *[in]*  **v**     vec4 (right, column vector)
+      | *[out]* **dest**  destination (result, column vector)
+
+    .. csv-table:: mat4x2
+        :header: "", "column 1", "column 2", "column 3"
+
+        "row 1", "m00", "m10", "m20"
+        "row 2", "m01", "m11", "m21"
+        "row 3", "m02", "m12", "m22"
+        "row 4", "m03", "m13", "m23"
+
+    .. csv-table:: column vec4 (1x4)
+        :header: "", "column 1"
+
+        "row 1", "v0"
+        "row 2", "v1"
+        "row 3", "v2"
+        "row 4", "v3"
+
+    .. csv-table:: column vec2 (1x2)
+        :header: "", "column 1"
+
+        "row 1", "m00 * v0 + m10 * v1 + m20 * v2 + m30 * v3"
+        "row 2", "m01 * v0 + m11 * v1 + m21 * v2 + m31 * v3"
+
+.. c:function:: void glm_mat4x2_transpose(mat4x2 m, mat2x4 dest)
+
+    transpose matrix and store in dest
+
+    Parameters:
+      | *[in]*  **m**     matrix
+      | *[out]* **dest**  destination
+
+.. c:function:: void  glm_mat4x2_scale(mat4x2 m, float s)
+
+    multiply matrix with scalar
+
+    Parameters:
+      | *[in, out]* **m** matrix
+      | *[in]*      **s** scalar

docs/source/mat4x3.rst

@@ -0,0 +1,144 @@
+.. default-domain:: C
+
+mat4x3
+======
+
+Header: cglm/mat4x3.h
+
+Table of contents (click to go):
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+Macros:
+
+1. GLM_MAT4X3_ZERO_INIT
+#. GLM_MAT4X3_ZERO
+
+Functions:
+
+1. :c:func:`glm_mat4x3_copy`
+#. :c:func:`glm_mat4x3_zero`
+#. :c:func:`glm_mat4x3_make`
+#. :c:func:`glm_mat4x3_mul`
+#. :c:func:`glm_mat4x3_mulv`
+#. :c:func:`glm_mat4x3_transpose`
+#. :c:func:`glm_mat4x3_scale`
+
+Represented
+~~~~~~~~~~~
+
+.. csv-table:: mat4x3
+   :header: "", "column 1", "column 2", "column 3", "column4"
+
+   "row 1", "m00", "m10", "m20", "m30"
+   "row 2", "m01", "m11", "m21", "m31"
+   "row 3", "m02", "m12", "m22", "m32"
+
+Functions documentation
+~~~~~~~~~~~~~~~~~~~~~~~
+
+.. c:function:: void glm_mat4x3_copy(mat4x3 mat, mat4x3 dest)
+
+    copy mat4x3 to another one (dest).
+
+    Parameters:
+      | *[in]*  **mat**   source
+      | *[out]* **dest**  destination
+
+.. c:function:: void glm_mat4x3_zero(mat4x3 mat)
+
+    make given matrix zero
+
+    Parameters:
+      | *[in,out]* **mat**  matrix
+
+.. c:function:: void glm_mat4x3_make(const float * __restrict src, mat4x3 dest)
+
+    Create mat4x3 matrix from pointer
+
+    .. note:: **@src** must contain at least 12 elements.
+
+    Parameters:
+      | *[in]*  **src**  pointer to an array of floats
+      | *[out]* **dest** destination matrix4x3
+
+.. c:function:: void glm_mat4x3_mul(mat4x3 m1, mat3x4 m2, mat3 dest)
+
+    multiply m1 and m2 to dest
+
+    .. code-block:: c
+
+       glm_mat4x3_mul(mat4x3, mat3x4, mat3);
+
+    Parameters:
+      | *[in]*  **m1**    left matrix (mat4x3)
+      | *[in]*  **m2**    right matrix (mat3x4)
+      | *[out]* **dest**  destination matrix (mat3)
+
+    .. csv-table:: mat4x3
+        :header: "", "column 1", "column 2", "column 3", "column 4"
+
+        "row 1", "a00", "a10", "a20", "a30"
+        "row 2", "a01", "a11", "a21", "a31"
+        "row 3", "a02", "a12", "a22", "a32"
+
+    .. csv-table:: mat3x4
+        :header: "", "column 1", "column 2", "column 3"
+
+        "row 1", "b00", "b10", "b20"
+        "row 2", "b01", "b11", "b21"
+        "row 3", "b02", "b12", "b22"
+        "row 4", "b03", "b13", "b23"
+
+    .. csv-table:: mat3x3
+        :header: "", "column 1", "column 2", "column 3"
+
+        "row 1", "a00 * b00 + a10 * b01 + a20 * b02 + a30 * b03", "a00 * b10 + a10 * b11 + a20 * b12 + a30 * b13", "a00 * b20 + a10 * b21 + a20 * b22 + a30 * b23"
+        "row 2", "a01 * b00 + a11 * b01 + a21 * b02 + a31 * b03", "a01 * b10 + a11 * b11 + a21 * b12 + a31 * b13", "a01 * b20 + a11 * b21 + a21 * b22 + a31 * b23"
+        "row 3", "a02 * b00 + a12 * b01 + a22 * b02 + a32 * b03", "a02 * b10 + a12 * b11 + a22 * b12 + a32 * b13", "a02 * b20 + a12 * b21 + a22 * b22 + a32 * b23"
+
+.. c:function:: void glm_mat4x3_mulv(mat4x3 m, vec4 v, vec3 dest)
+
+    multiply mat4x3 with vec4 (column vector) and store in dest column vector
+
+    Parameters:
+      | *[in]*  **m**     mat4x3 (left)
+      | *[in]*  **v**     vec4 (right, column vector)
+      | *[out]* **dest**  destination (result, column vector)
+
+    .. csv-table:: mat4x3
+        :header: "", "column 1", "column 2", "column 3", "column 4"
+
+        "row 1", "m00", "m10", "m20", "m30"
+        "row 2", "m01", "m11", "m21", "m31"
+        "row 3", "m02", "m12", "m22", "m32"
+
+    .. csv-table:: column vec4 (1x4)
+        :header: "", "column 1"
+
+        "row 1", "v0"
+        "row 2", "v1"
+        "row 3", "v2"
+        "row 4", "v3"
+
+    .. csv-table:: column vec3 (1x3)
+        :header: "", "column 1"
+
+        "row 1", "m00 * v0 + m10 * v1 + m20 * v2 + m30 * v3"
+        "row 2", "m01 * v0 + m11 * v1 + m21 * v2 + m31 * v3"
+        "row 3", "m02 * v0 + m12 * v1 + m22 * v2 + m32 * v3"
+
+.. c:function:: void glm_mat4x3_transpose(mat4x3 m, mat3x4 dest)
+
+    transpose matrix and store in dest
+
+    Parameters:
+      | *[in]*  **m**     matrix
+      | *[out]* **dest**  destination
+
+.. c:function:: void  glm_mat4x3_scale(mat4x3 m, float s)
+
+    multiply matrix with scalar
+
+    Parameters:
+      | *[in, out]* **m** matrix
+      | *[in]*      **s** scalar

docs/source/opengl.rst

@@ -43,14 +43,14 @@ array of matrices:
    /* ... */
    glUniformMatrix4fv(location, count, GL_FALSE, (float *)matrix);
 
-in this way, passing aray of matrices is same
+in this way, passing array of matrices is same
 
 Passing / Uniforming Vectors to OpenGL:
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
 You don't need to do extra thing when passing cglm vectors to OpengL or other APIs.
 Because a function like **glUniform4fv** accepts vector as pointer. cglm's vectors
-are array of floats. So you can pass it directly ot those functions:
+are array of floats. So you can pass it directly to those functions:
 
 .. code-block:: c
 

docs/source/opt.rst

@@ -1,39 +1,37 @@
 .. default-domain:: C
 
-Options
+🛠️ Options
 ===============================================================================
 
 A few options are provided via macros.
 
+❗️ IMPORTANT ❗️
+
+It's a good idea to set up your config macros in build settings like CMake, Xcode, or Visual Studio. This is especially important if you're using features like Modules in Xcode, where adding macros directly before the **cglm** headers might not work.
+
 Alignment Option
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
-As default, cglm requires types to be aligned. Alignment requirements:
+By default, **cglm** requires types to be aligned with specific byte requirements:
 
-vec3:   8 byte
-vec4:   16 byte
-mat4:   16 byte
-versor: 16 byte
+- vec3: 8 bytes
+- vec4: 16 bytes
+- mat4: 16 bytes (32 on AVX)
+- versor: 16 bytes
 
-By starting **v0.4.5** cglm provides an option to disable alignment requirement.
-To enable this option define **CGLM_ALL_UNALIGNED** macro before all headers.
-You can define it in Xcode, Visual Studio (or other IDEs) or you can also prefer
-to define it in build system. If you use pre-compiled versions then you
-have to compile cglm with **CGLM_ALL_UNALIGNED** macro.
+Starting with **v0.4.5**, **cglm** offers an option to relax these alignment requirements. To use this option, define the **CGLM_ALL_UNALIGNED** macro before including any headers. This definition can be made within Xcode, Visual Studio, other IDEs, or directly in your build system. If using pre-compiled versions of **cglm**, you'll need to compile them with the **CGLM_ALL_UNALIGNED** macro.
 
-**VERY VERY IMPORTANT:** If you use cglm in multiple projects and
- those projects are depends on each other, then
+**❗️NOTE:❗️** If you're using **cglm** across multiple interdependent projects:
 
- | *ALWAYS* or *NEVER USE* **CGLM_ALL_UNALIGNED** macro in linked projects
+- Always or never use the **CGLM_ALL_UNALIGNED** macro in all linked projects to avoid configuration conflicts. A **cglm** header from one project could require alignment, while a header from another might not, leading to **cglm** functions accessing invalid memory locations.
 
- if you do not know what you are doing. Because a cglm header included
- via 'project A' may force types to be aligned and another cglm header
- included via 'project B' may not require alignment. In this case
- cglm functions will read from and write to **INVALID MEMORY LOCATIONs**.
+- **Key Point:** Maintain the same **cglm** configuration across all your projects. For example, if you activate **CGLM_ALL_UNALIGNED** in one project, ensure it's set in the others too.
 
- ALWAYS USE SAME CONFIGURATION / OPTION for **cglm** if you have multiple projects.
+**❗️NOTE:❗️**
 
- For instance if you set CGLM_ALL_UNALIGNED in a project then set it in other projects too
+While **CGLM_ALL_UNALIGNED** allows for flexibility in alignment, it doesn't override C's fundamental alignment rules. For example, an array like *vec4* decays to a pointer (float*) in functions, which must adhere to the alignment requirements of a float pointer (4 bytes). This adherence is crucial because **cglm** directly dereferences these pointers instead of copying data, and failing to meet alignment requirements can lead to unpredictable errors, such as crashes.
+
+You can use `CGLM_ALIGN` and `CGLM_ALIGN_MAT` macros for aligning local variables or struct members. However, when dealing with dynamic memory allocation or custom memory locations, you'll need to ensure alignment requirements are met appropriately for those cases
 
 Clipspace Option[s]
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
@@ -43,7 +41,7 @@ By starting **v0.8.3** cglm provides options to switch between clipspace configu
 Clipspace related files are located at `include/cglm/[struct]/clipspace.h` but 
 these are included in related files like `cam.h`. If you don't want to change your existing 
 clipspace configuration and want to use different clipspace function like `glm_lookat_zo` or `glm_lookat_lh_zo`...
-then you can include individual headers or just define `CGLM_CLIPSPACE_INCLUDE_ALL` which will iclude all headers for you.
+then you can include individual headers or just define `CGLM_CLIPSPACE_INCLUDE_ALL` which will include all headers for you.
 
 1. **CGLM_CLIPSPACE_INCLUDE_ALL**
 2. **CGLM_FORCE_DEPTH_ZERO_TO_ONE**
@@ -90,6 +88,16 @@ You have to extra options for dot product: **CGLM_SSE4_DOT** and **CGLM_SSE3_DOT
 
 otherwise cglm will use custom cglm's hadd functions which are optimized too.
 
+Struct API Options
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+To configure the Struct API namespace, you can define the following macros before including the cglm/struct.h header:
+
+- **CGLM_OMIT_NS_FROM_STRUCT_API**: omits CGLM_STRUCT_API_NS (`glms_`) namespace completely if there is sub namespace e.g `mat4_`, `vec4_` ... DEFAULT is not defined
+- **CGLM_STRUCT_API_NS**:           define name space for struct api, DEFAULT is **glms**
+- **CGLM_STRUCT_API_NAME_SUFFIX**:  define name suffix, DEFAULT is **empty** e.g defining it as #define CGLM_STRUCT_API_NAME_SUFFIX  s will add s suffix to mat4_mul -> mat4s_mul
+
+
 Print Options
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 

docs/source/quat.rst

@@ -62,6 +62,7 @@ Functions:
 #. :c:func:`glm_quat_rotate`
 #. :c:func:`glm_quat_rotate_at`
 #. :c:func:`glm_quat_rotate_atm`
+#. :c:func:`glm_quat_make`
 
 Functions documentation
 ~~~~~~~~~~~~~~~~~~~~~~~
@@ -420,3 +421,13 @@ Functions documentation
       | *[in, out]*  **m**      existing transform matrix to rotate
       | *[in]*       **q**      quaternion
       | *[in]*       **pivot**  pivot
+
+.. c:function:: void glm_quat_make(const float * __restrict src, versor dest)
+
+    Create quaternion from pointer
+
+    .. note:: **@src** must contain at least 4 elements. cglm store quaternions as [x, y, z, w].
+
+    Parameters:
+      | *[in]*  **src**  pointer to an array of floats
+      | *[out]* **dest** destination quaternion

docs/source/ray.rst

@@ -13,6 +13,8 @@ Table of contents (click to go):
 Functions:
 
 1. :c:func:`glm_ray_triangle`
+#. :c:func:`glm_ray_sphere`
+#. :c:func:`glm_ray_at`
 
 Functions documentation
 ~~~~~~~~~~~~~~~~~~~~~~~
@@ -29,3 +31,39 @@ Functions documentation
       | *[in]*       **v2**            third vertex of triangle
       | *[in, out]*  **d**             float pointer to save distance to intersection
       | *[out]*      **intersection**  whether there is intersection
+
+.. c:function:: bool glm_ray_sphere(vec3 origin, vec3 dir, vec4 s, float * __restrict t1, float * __restrict t2)
+
+    ray sphere intersection
+
+    returns false if there is no intersection if true:
+
+    - t1 > 0, t2 > 0: ray intersects the sphere at t1 and t2 both ahead of the origin
+    - t1 < 0, t2 > 0: ray starts inside the sphere, exits at t2
+    - t1 < 0, t2 < 0: no intersection ahead of the ray ( returns false )
+    - the caller can check if the intersection points (t1 and t2) fall within a
+      specific range (for example, tmin < t1, t2 < tmax) to determine if the
+      intersections are within a desired segment of the ray
+
+    Parameters:
+      | *[in]*   **origin**  ray origin
+      | *[in]*   **dir**     normalized ray direction
+      | *[in]*   **s**       sphere  [center.x, center.y, center.z, radii]
+      | *[out]*  **t1**      near point1 (closer to origin)
+      | *[out]*  **t2**      far point2 (farther from origin)
+
+    Return:
+      | whether there is intersection
+
+.. c:function:: bool  glm_ray_at(vec3 orig, vec3 dir, float t, vec3 point)
+
+    point using t by 𝐏(𝑡)=𝐀+𝑡𝐛
+
+    Parameters:
+      | *[in]*   **origin**  ray origin
+      | *[in]*   **dir**     ray direction
+      | *[out]*  **t**       parameter
+      | *[out]*  **point**   point at t
+
+    Return:
+      | point at t

docs/source/sphere.rst

@@ -9,7 +9,7 @@ Header: cglm/sphere.h
 
 Sphere Representation in cglm is *vec4*: **[center.x, center.y, center.z, radii]**
 
-You can call any vec3 function by pasing sphere. Because first three elements
+You can call any vec3 function by passing sphere. Because first three elements
 defines center of sphere.
 
 Table of contents (click to go):

docs/source/sphinx-static/theme_overrides.css

@@ -0,0 +1,12 @@
+@media screen {
+  /* content column
+   *
+   * RTD theme's default is 800px as max width for the content, but we have
+   * tables with tons of columns, which need the full width of the view-port.
+   *
+   * Comment from yocto project theme_overrides.css
+   */
+
+  .wy-nav-content{ max-width: none; }
+
+}

docs/source/troubleshooting.rst

@@ -3,30 +3,30 @@
 Troubleshooting
 ================================================================================
 
-It is possible that sometimes you may get crashes or wrong results.
+It is possible that you may sometimes get crashes or wrong results.
 Follow these topics
 
 Memory Allocation:
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
-Again, **cglm** doesn't alloc any memory on heap.
-cglm functions works like memcpy; it copies data from src,
-makes calculations then copy the result to dest.
+Recall that **cglm** does not alloc any memory on the heap.
+cglm functions work like memcpy; they copy data from src,
+make calculations, then copy the result to dest.
 
 You are responsible for allocation of **src** and **dest** parameters.
 
 Alignment:
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
-**vec4** and **mat4** types requires 16 byte alignment.
-These types are marked with align attribute to let compiler know about this
+**vec4** and **mat4** types require 16 byte alignment.
+These types are marked with the align attribute to let the compiler know about this
 requirement.
 
-But since MSVC (Windows) throws the error:
+Since MSVC (Windows) throws this error:
 
 **"formal parameter with requested alignment of 16 won't be aligned"**
 
-The alignment attribute has been commented for MSVC
+The alignment attribute has been commented out for MSVC
 
 .. code-block:: c
 
@@ -37,61 +37,61 @@ The alignment attribute has been commented for MSVC
    #endif.
 
 So MSVC may not know about alignment requirements when creating variables.
-The interesting thing is that, if I remember correctly Visual Studio 2017
+The interesting thing is that, if I remember correctly, Visual Studio 2017
 doesn't throw the above error. So we may uncomment that line for Visual Studio 2017,
 you may do it yourself.
 
 **This MSVC issue is still in TODOs.**
 
-**UPDATE:** By starting v0.4.5 cglm provides an option to disable alignment requirement.
-Also alignment is disabled for older msvc verisons as default. Now alignment is only required in Visual Studio 2017 version 15.6+ if CGLM_ALL_UNALIGNED macro is not defined.
+**UPDATE:** Starting with v0.4.5, cglm provides an option to disable the alignment requirement.
+Also, alignment is disabled for older msvc versions by default. Now alignment is only required in Visual Studio 2017 version 15.6+ if the CGLM_ALL_UNALIGNED macro is not defined.
 
 Crashes, Invalid Memory Access:
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
-Probably you are trying to write to invalid memory location.
+Most likely, you are trying to write to an invalid memory location.
 
-You may used wrong function for what you want to do.
+You may have used a wrong function for what you want to do.
 
-For instance you may called **glm_vec4_** functions for **vec3** data type.
-It will try to write 32 byte but since **vec3** is 24 byte it should throw
-memory access error or exit the app without saying anything.
+For example, you may have called a **glm_vec4_** function for a **vec3** data type.
+It will try to write 32 bytes, but since **vec3** is 24 bytes, it should throw
+a memory access error or exit the app without saying anything.
 
 **UPDATE - IMPORTANT:** 
 
   | On MSVC or some other compilers, if alignment is enabled (default) then double check alignment requirements if you got a crash.
 
-  | If you send GLM_VEC4_ONE or similar macros directly to a function, it may be crashed.
-  | Because compiler may not apply alignment as defined on **typedef** to that macro while passing it (on stack) to a function.
+  | If you send GLM_VEC4_ONE or similar macros directly to a function, it may crash.
+  | Because the compiler may not apply alignment as defined on **typedef** to that macro while passing it (on stack) to a function.
 
 Wrong Results:
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
-Again, you may used wrong function.
+Again, you may have used a wrong function.
 
-For instance if you use **glm_normalize()** or **glm_vec3_normalize()** for **vec4**,
-it will assume that passed param is **vec3** and will normalize it for **vec3**.
-Since you need to **vec4** to be normalized in your case, you will get wrong results.
+For instance if you use **glm_normalize()** or **glm_vec3_normalize()** for a **vec4**,
+it will assume that the passed param is a **vec3**, and will normalize it for **vec3**.
+Since you need a **vec4** to be normalized in your case, you will get wrong results.
 
-Accessing vec4 type with vec3 functions is valid, you will not get any error, exception or crash.
+Accessing a vec4 type with vec3 functions is valid, you will not get any error, exception or crash.
 You only get wrong results if you don't know what you are doing!
 
-So be carefull, when your IDE (Xcode, Visual Studio ...) tried to autocomplete function names, READ IT :)
+So be careful, when your IDE (Xcode, Visual Studio ...) tries to autocomplete function names, READ IT :)
 
-**Also implementation may be wrong please let us know by creating an issue on Github.**
+**Also implementation may be wrong, please let us know by creating an issue on Github.**
 
 BAD_ACCESS : Thread 1: EXC_BAD_ACCESS (code=EXC_I386_GPFLT) or Similar Errors/Crashes
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
-This is similar issue with alignment. For instance if you compiled **cglm** with 
+This is a similar issue with alignment. For instance if you compiled **cglm** with 
 AVX (**-mavx**, intentionally or not) and if you use **cglm** in an environment that doesn't 
 support AVX (or if AVX is disabled intentionally) e.g. environment that max support SSE2/3/4, 
 then you probably get **BAD ACCESS** or similar...
 
 Because if you compile **cglm** with AVX it aligns **mat4** with 32 byte boundary, 
-and your project aligns that as 16 byte boundary...
+and your project aligns that as a 16 byte boundary...
 
-Check alignment, supported vector extension or simd in **cglm** and linked projects...
+Check alignment, supported vector extension, or simd in **cglm** and linked projects...
 
 Other Issues?
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

docs/source/util.rst

@@ -68,14 +68,14 @@ Functions documentation
 
 .. c:function:: void  glm_make_rad(float *degm)
 
-    | convert exsisting degree to radians. this will override degrees value
+    | convert existing degree to radians. this will override degrees value
 
     Parameters:
       | *[in, out]*  **deg**      pointer to angle in degrees
 
 .. c:function:: void  glm_make_deg(float *rad)
 
-    | convert exsisting radians to degree. this will override radians value
+    | convert existing radians to degree. this will override radians value
 
     Parameters:
       | *[in, out]*  **rad**      pointer to angle in radians

docs/source/vec2.rst

@@ -10,10 +10,10 @@ Table of contents (click to go):
 
 Macros:
 
-1. GLM_vec2_ONE_INIT
-#. GLM_vec2_ZERO_INIT
-#. GLM_vec2_ONE
-#. GLM_vec2_ZERO
+1. GLM_VEC2_ONE_INIT
+#. GLM_VEC2_ZERO_INIT
+#. GLM_VEC2_ONE
+#. GLM_VEC2_ZERO
 
 Functions:
 
@@ -45,12 +45,16 @@ Functions:
 #. :c:func:`glm_vec2_normalize`
 #. :c:func:`glm_vec2_normalize_to`
 #. :c:func:`glm_vec2_rotate`
+#. :c:func:`glm_vec2_center`
 #. :c:func:`glm_vec2_distance2`
 #. :c:func:`glm_vec2_distance`
 #. :c:func:`glm_vec2_maxv`
 #. :c:func:`glm_vec2_minv`
 #. :c:func:`glm_vec2_clamp`
 #. :c:func:`glm_vec2_lerp`
+#. :c:func:`glm_vec2_make`
+#. :c:func:`glm_vec2_reflect`
+#. :c:func:`glm_vec2_refract`
 
 Functions documentation
 ~~~~~~~~~~~~~~~~~~~~~~~
@@ -115,7 +119,7 @@ Functions documentation
     norm * norm (magnitude) of vector
 
     we can use this func instead of calling norm * norm, because it would call
-    sqrtf fuction twice but with this func we can avoid func call, maybe this is
+    sqrtf function twice but with this func we can avoid func call, maybe this is
     not good name for this func
 
     Parameters:
@@ -313,6 +317,15 @@ Functions documentation
       | *[in]*  **axis**   axis vector
       | *[out]* **dest**   destination
 
+.. c:function:: void  glm_vec2_center(vec2 v1, vec2 v2, vec2 dest)
+
+    find center point of two vector
+
+    Parameters:
+      | *[in]*  **v1**    vector1
+      | *[in]*  **v2**    vector2
+      | *[out]* **dest**  center point
+
 .. c:function:: float glm_vec2_distance2(vec2 v1, vec2 v2)
 
     squared distance between two vectors
@@ -373,3 +386,39 @@ Functions documentation
       | *[in]*  **to**     to value
       | *[in]*  **t**      interpolant (amount) clamped between 0 and 1
       | *[out]* **dest**   destination
+
+.. c:function:: void glm_vec2_make(const float * __restrict src, vec2 dest)
+
+    Create two dimensional vector from pointer
+
+    .. note:: **@src** must contain at least 2 elements.
+
+    Parameters:
+      | *[in]*  **src**  pointer to an array of floats
+      | *[out]* **dest** destination vector
+
+.. c:function:: void glm_vec2_reflect(vec2 v, vec2 n, vec2 dest)
+
+    Reflection vector using an incident ray and a surface normal
+
+    Parameters:
+      | *[in]*  **v**     incident vector
+      | *[in]*  **n**     *❗️ normalized ❗️* normal vector
+      | *[out]* **dest**  destination: reflection result
+
+.. c:function:: bool glm_vec2_refract(vec2 v, vec2 n, float eta, vec2 dest)
+
+    Computes refraction vector for an incident vector and a surface normal.
+   
+    Calculates the refraction vector based on Snell's law. If total internal reflection
+    occurs (angle too great given eta), dest is set to zero and returns false.
+    Otherwise, computes refraction vector, stores it in dest, and returns true.
+
+    Parameters:
+      | *[in]*  **v**     *❗️ normalized ❗️* incident vector
+      | *[in]*  **n**     *❗️ normalized ❗️* normal vector
+      | *[in]*  **eta**   ratio of indices of refraction (incident/transmitted)
+      | *[out]* **dest**  refraction vector if refraction occurs; zero vector otherwise
+
+    Returns:
+      returns true if refraction occurs; false if total internal reflection occurs.

docs/source/vec3-ext.rst

@@ -50,6 +50,14 @@ Functions documentation
       | *[in]*  **val**   value
       | *[out]* **dest**  destination
 
+.. c:function:: void  glm_vec3_fill(vec3 v, float val)
+
+    fill a vector with specified value
+
+    Parameters:
+      | *[out]*  **v**    vector
+      | *[in]*   **val**  value
+
 .. c:function:: bool  glm_vec3_eq(vec3 v, float val)
 
     check if vector is equal to value (without epsilon)

docs/source/vec3.rst

@@ -13,7 +13,7 @@ Header: cglm/vec3.h
 We mostly use vectors in graphics math, to make writing code faster
 and easy to read, some *vec3* functions are aliased in global namespace.
 For instance :c:func:`glm_dot` is alias of :c:func:`glm_vec3_dot`,
-alias means inline wrapper here. There is no call verison of alias functions
+alias means inline wrapper here. There is no call version of alias functions
 
 There are also functions for rotating *vec3* vector. **_m4**, **_m3** prefixes
 rotate *vec3* with matrix.
@@ -79,6 +79,10 @@ Functions:
 #. :c:func:`glm_vec3_ortho`
 #. :c:func:`glm_vec3_clamp`
 #. :c:func:`glm_vec3_lerp`
+#. :c:func:`glm_vec3_make`
+#. :c:func:`glm_vec3_faceforward`
+#. :c:func:`glm_vec3_reflect`
+#. :c:func:`glm_vec3_refract`
 
 Functions documentation
 ~~~~~~~~~~~~~~~~~~~~~~~
@@ -147,7 +151,7 @@ Functions documentation
     norm * norm (magnitude) of vector
 
     we can use this func instead of calling norm * norm, because it would call
-    sqrtf fuction twice but with this func we can avoid func call, maybe this is
+    sqrtf function twice but with this func we can avoid func call, maybe this is
     not good name for this func
 
     Parameters:
@@ -308,7 +312,7 @@ Functions documentation
 
 .. c:function:: void  glm_vec3_flipsign(vec3 v)
 
-    **DEPRACATED!**
+    **DEPRECATED!**
 
     use :c:func:`glm_vec3_negate`
 
@@ -317,7 +321,7 @@ Functions documentation
 
 .. c:function:: void  glm_vec3_flipsign_to(vec3 v, vec3 dest)
 
-    **DEPRACATED!**
+    **DEPRECATED!**
 
     use :c:func:`glm_vec3_negate_to`
 
@@ -327,7 +331,7 @@ Functions documentation
 
 .. c:function:: void  glm_vec3_inv(vec3 v)
 
-    **DEPRACATED!**
+    **DEPRECATED!**
 
     use :c:func:`glm_vec3_negate`
 
@@ -336,7 +340,7 @@ Functions documentation
 
 .. c:function:: void  glm_vec3_inv_to(vec3 v, vec3 dest)
 
-    **DEPRACATED!**
+    **DEPRECATED!**
 
     use :c:func:`glm_vec3_negate_to`
 
@@ -376,7 +380,7 @@ Functions documentation
 
 .. c:function:: float  glm_vec3_angle(vec3 v1, vec3 v2)
 
-    angle betwen two vector
+    angle between two vector
 
     Parameters:
       | *[in]*  **v1**   vector1
@@ -501,3 +505,50 @@ Functions documentation
       | *[in]*  **to**     to value
       | *[in]*  **t**      interpolant (amount) clamped between 0 and 1
       | *[out]* **dest**   destination
+
+.. c:function:: void glm_vec3_make(const float * __restrict src, vec3 dest)
+
+    Create three dimensional vector from pointer
+
+    .. note::: **@src** must contain at least 3 elements.
+
+    Parameters:
+      | *[in]*  **src**  pointer to an array of floats
+      | *[out]* **dest** destination vector
+
+.. c:function:: void glm_vec3_faceforward(vec3 n, vec3 v, vec3 nref, vec3 dest)
+
+    A vector pointing in the same direction as another
+
+    Parameters:
+      | *[in]*  **n**     vector to orient
+      | *[in]*  **v**     incident vector
+      | *[in]*  **nref**  reference vector
+      | *[out]* **dest**  destination: oriented vector, pointing away from the surface.
+
+.. c:function:: void glm_vec3_reflect(vec3 v, vec3 n, vec3 dest)
+
+    Reflection vector using an incident ray and a surface normal
+
+    Parameters:
+      | *[in]*  **v**     incident vector
+      | *[in]*  **n**     *❗️ normalized ❗️* normal vector
+      | *[out]* **dest**  destination: reflection result
+
+.. c:function:: bool glm_vec3_refract(vec3 v, vec3 n, float eta, vec3 dest)
+
+    
+    Computes refraction vector for an incident vector and a surface normal.
+   
+    Calculates the refraction vector based on Snell's law. If total internal reflection
+    occurs (angle too great given eta), dest is set to zero and returns false.
+    Otherwise, computes refraction vector, stores it in dest, and returns true.
+
+    Parameters:
+      | *[in]*  **v**     *❗️ normalized ❗️* incident vector
+      | *[in]*  **n**     *❗️ normalized ❗️* normal vector
+      | *[in]*  **eta**   ratio of indices of refraction (incident/transmitted)
+      | *[out]* **dest**  refraction vector if refraction occurs; zero vector otherwise
+
+    Returns:
+      returns true if refraction occurs; false if total internal reflection occurs.

docs/source/vec4.rst

@@ -59,6 +59,9 @@ Functions:
 #. :c:func:`glm_vec4_clamp`
 #. :c:func:`glm_vec4_lerp`
 #. :c:func:`glm_vec4_cubic`
+#. :c:func:`glm_vec4_make`
+#. :c:func:`glm_vec4_reflect`
+#. :c:func:`glm_vec4_refract`
 
 Functions documentation
 ~~~~~~~~~~~~~~~~~~~~~~~
@@ -107,6 +110,13 @@ Functions documentation
     Parameters:
       | *[in, out]*  **v**     vector
 
+.. c:function:: void  glm_vec4_one(vec4 v)
+
+    makes all members one
+
+    Parameters:
+      | *[in, out]*  **v**     vector
+
 .. c:function:: float  glm_vec4_dot(vec4 a, vec4 b)
 
     dot product of vec4
@@ -123,7 +133,7 @@ Functions documentation
     norm * norm (magnitude) of vector
 
     we can use this func instead of calling norm * norm, because it would call
-    sqrtf fuction twice but with this func we can avoid func call, maybe this is
+    sqrtf function twice but with this func we can avoid func call, maybe this is
     not good name for this func
 
     Parameters:
@@ -283,7 +293,7 @@ Functions documentation
 
 .. c:function:: void  glm_vec4_flipsign(vec4 v)
 
-    **DEPRACATED!**
+    **DEPRECATED!**
 
     use :c:func:`glm_vec4_negate`
 
@@ -292,7 +302,7 @@ Functions documentation
 
 .. c:function:: void  glm_vec4_flipsign_to(vec4 v, vec4 dest)
 
-    **DEPRACATED!**
+    **DEPRECATED!**
 
     use :c:func:`glm_vec4_negate_to`
 
@@ -302,7 +312,7 @@ Functions documentation
 
 .. c:function:: void  glm_vec4_inv(vec4 v)
 
-    **DEPRACATED!**
+    **DEPRECATED!**
 
     use :c:func:`glm_vec4_negate`
 
@@ -311,7 +321,7 @@ Functions documentation
 
 .. c:function:: void  glm_vec4_inv_to(vec4 v, vec4 dest)
 
-    **DEPRACATED!**
+    **DEPRECATED!**
 
     use :c:func:`glm_vec4_negate_to`
 
@@ -406,3 +416,43 @@ Functions documentation
     Parameters:
       | *[in]*  **s**      parameter
       | *[out]* **dest**   destination
+
+.. c:function:: void glm_vec4_make(const float * __restrict src, vec4 dest)
+
+    Create four dimensional vector from pointer
+
+    .. note:: **@src** must contain at least 4 elements.
+
+    Parameters:
+      | *[in]*  **src**  pointer to an array of floats
+      | *[out]* **dest** destination vector
+
+.. c:function:: bool glm_vec4_reflect(vec4 v, vec4 n, vec4 dest)
+
+    Reflection vector using an incident ray and a surface normal
+
+    Parameters:
+      | *[in]*  **v**     incident vector
+      | *[in]*  **n**     *❗️ normalized ❗️* normal vector
+      | *[out]* **dest**  destination: reflection result
+
+.. c:function:: bool glm_vec4_refract(vec4 v, vec4 n, float eta, vec4 dest)
+
+    computes refraction vector for an incident vector and a surface normal.
+   
+    Calculates the refraction vector based on Snell's law. If total internal reflection
+    occurs (angle too great given eta), dest is set to zero and returns false.
+    Otherwise, computes refraction vector, stores it in dest, and returns true.
+   
+    This implementation does not explicitly preserve the 'w' component of the
+    incident vector 'I' in the output 'dest', users requiring the preservation of
+    the 'w' component should manually adjust 'dest' after calling this function.
+
+    Parameters:
+      | *[in]*  **v**     *❗️ normalized ❗️* incident vector
+      | *[in]*  **n**     *❗️ normalized ❗️* normal vector
+      | *[in]*  **eta**   ratio of indices of refraction (incident/transmitted)
+      | *[out]* **dest**  refraction vector if refraction occurs; zero vector otherwise
+
+    Returns:
+      returns true if refraction occurs; false if total internal reflection occurs.

include/cglm/aabb2d.h

@@ -0,0 +1,271 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglm_aabb2d_h
+#define cglm_aabb2d_h
+
+#include "common.h"
+#include "vec2.h"
+#include "util.h"
+
+/* DEPRECATED! use _diag */
+#define glm_aabb2d_size(aabb)         glm_aabb2d_diag(aabb)
+
+/*!
+ * @brief make [aabb] zero
+ *
+ * @param[in, out]  aabb aabb
+ */
+CGLM_INLINE
+void
+glm_aabb2d_zero(vec2 aabb[2]) {
+  glm_vec2_zero(aabb[0]);
+  glm_vec2_zero(aabb[1]);
+}
+
+/*!
+ * @brief copy all members of [aabb] to [dest]
+ *
+ * @param[in]  aabb source
+ * @param[out] dest destination
+ */
+CGLM_INLINE
+void
+glm_aabb2d_copy(vec2 aabb[2], vec2 dest[2]) {
+  glm_vec2_copy(aabb[0], dest[0]);
+  glm_vec2_copy(aabb[1], dest[1]);
+}
+
+/*!
+ * @brief apply transform to Axis-Aligned Bounding aabb
+ *
+ * @param[in]  aabb  bounding aabb
+ * @param[in]  m    transform matrix
+ * @param[out] dest transformed bounding aabb
+ */
+CGLM_INLINE
+void
+glm_aabb2d_transform(vec2 aabb[2], mat3 m, vec2 dest[2]) {
+  vec2 v[2], xa, xb, ya, yb;
+
+  glm_vec2_scale(m[0], aabb[0][0], xa);
+  glm_vec2_scale(m[0], aabb[1][0], xb);
+
+  glm_vec2_scale(m[1], aabb[0][1], ya);
+  glm_vec2_scale(m[1], aabb[1][1], yb);
+
+  /* translation + min(xa, xb) + min(ya, yb) */
+  glm_vec2(m[2], v[0]);
+  glm_vec2_minadd(xa, xb, v[0]);
+  glm_vec2_minadd(ya, yb, v[0]);
+
+  /* translation + max(xa, xb) + max(ya, yb) */
+  glm_vec2(m[2], v[1]);
+  glm_vec2_maxadd(xa, xb, v[1]);
+  glm_vec2_maxadd(ya, yb, v[1]);
+
+  glm_vec2_copy(v[0], dest[0]);
+  glm_vec2_copy(v[1], dest[1]);
+}
+
+/*!
+ * @brief merges two AABB bounding aabb and creates new one
+ *
+ * two aabb must be in same space, if one of aabb is in different space then
+ * you should consider to convert it's space by glm_aabb_space
+ *
+ * @param[in]  aabb1 bounding aabb 1
+ * @param[in]  aabb2 bounding aabb 2
+ * @param[out] dest merged bounding aabb
+ */
+CGLM_INLINE
+void
+glm_aabb2d_merge(vec2 aabb1[2], vec2 aabb2[2], vec2 dest[2]) {
+  dest[0][0] = glm_min(aabb1[0][0], aabb2[0][0]);
+  dest[0][1] = glm_min(aabb1[0][1], aabb2[0][1]);
+
+  dest[1][0] = glm_max(aabb1[1][0], aabb2[1][0]);
+  dest[1][1] = glm_max(aabb1[1][1], aabb2[1][1]);
+}
+
+/*!
+ * @brief crops a bounding aabb with another one.
+ *
+ * this could be useful for gettng a baabb which fits with view frustum and
+ * object bounding aabbes. In this case you crop view frustum aabb with objects
+ * aabb
+ *
+ * @param[in]  aabb     bounding aabb 1
+ * @param[in]  cropAabb crop aabb
+ * @param[out] dest     cropped bounding aabb
+ */
+CGLM_INLINE
+void
+glm_aabb2d_crop(vec2 aabb[2], vec2 cropAabb[2], vec2 dest[2]) {
+  dest[0][0] = glm_max(aabb[0][0], cropAabb[0][0]);
+  dest[0][1] = glm_max(aabb[0][1], cropAabb[0][1]);
+
+  dest[1][0] = glm_min(aabb[1][0], cropAabb[1][0]);
+  dest[1][1] = glm_min(aabb[1][1], cropAabb[1][1]);
+}
+
+/*!
+ * @brief crops a bounding aabb with another one.
+ *
+ * this could be useful for gettng a baabb which fits with view frustum and
+ * object bounding aabbes. In this case you crop view frustum aabb with objects
+ * aabb
+ *
+ * @param[in]  aabb      bounding aabb
+ * @param[in]  cropAabb  crop aabb
+ * @param[in]  clampAabb minimum aabb
+ * @param[out] dest      cropped bounding aabb
+ */
+CGLM_INLINE
+void
+glm_aabb2d_crop_until(vec2 aabb[2],
+                      vec2 cropAabb[2],
+                      vec2 clampAabb[2],
+                      vec2 dest[2]) {
+  glm_aabb2d_crop(aabb, cropAabb, dest);
+  glm_aabb2d_merge(clampAabb, dest, dest);
+}
+
+/*!
+ * @brief invalidate AABB min and max values
+ *
+ * @param[in, out]  aabb bounding aabb
+ */
+CGLM_INLINE
+void
+glm_aabb2d_invalidate(vec2 aabb[2]) {
+  glm_vec2_fill(aabb[0], FLT_MAX);
+  glm_vec2_fill(aabb[1], -FLT_MAX);
+}
+
+/*!
+ * @brief check if AABB is valid or not
+ *
+ * @param[in]  aabb bounding aabb
+ */
+CGLM_INLINE
+bool
+glm_aabb2d_isvalid(vec2 aabb[2]) {
+  return glm_vec2_max(aabb[0]) != FLT_MAX
+         && glm_vec2_min(aabb[1]) != -FLT_MAX;
+}
+
+/*!
+ * @brief distance between of min and max
+ *
+ * @param[in]  aabb bounding aabb
+ */
+CGLM_INLINE
+float
+glm_aabb2d_diag(vec2 aabb[2]) {
+  return glm_vec2_distance(aabb[0], aabb[1]);
+}
+
+/*!
+ * @brief size of aabb
+ *
+ * @param[in]  aabb bounding aabb
+ * @param[out]  dest size
+ */
+CGLM_INLINE
+void
+glm_aabb2d_sizev(vec2 aabb[2], vec2 dest) {
+  glm_vec2_sub(aabb[1], aabb[0], dest); 
+}
+
+/*!
+ * @brief radius of sphere which surrounds AABB
+ *
+ * @param[in]  aabb bounding aabb
+ */
+CGLM_INLINE
+float
+glm_aabb2d_radius(vec2 aabb[2]) {
+  return glm_aabb2d_diag(aabb) * 0.5f;
+}
+
+/*!
+ * @brief computes center point of AABB
+ *
+ * @param[in]   aabb  bounding aabb
+ * @param[out]  dest center of bounding aabb
+ */
+CGLM_INLINE
+void
+glm_aabb2d_center(vec2 aabb[2], vec2 dest) {
+  glm_vec2_center(aabb[0], aabb[1], dest);
+}
+
+/*!
+ * @brief check if two AABB intersects
+ *
+ * @param[in]   aabb    bounding aabb
+ * @param[in]   other  other bounding aabb
+ */
+CGLM_INLINE
+bool
+glm_aabb2d_aabb(vec2 aabb[2], vec2 other[2]) {
+  return (aabb[0][0] <= other[1][0] && aabb[1][0] >= other[0][0])
+      && (aabb[0][1] <= other[1][1] && aabb[1][1] >= other[0][1]);
+}
+
+/*!
+ * @brief check if AABB intersects with a circle
+ *
+ * Circle Representation in cglm: [center.x, center.y, radii]
+ *
+ * @param[in]   aabb   solid bounding aabb
+ * @param[in]   c      solid circle
+ */
+CGLM_INLINE
+bool
+glm_aabb2d_circle(vec2 aabb[2], vec3 c) {
+  float dmin;
+  int   a, b;
+
+  a = (c[0] < aabb[0][0]) + (c[0] > aabb[1][0]);
+  b = (c[1] < aabb[0][1]) + (c[1] > aabb[1][1]);
+
+  dmin  = glm_pow2((c[0] - aabb[!(a - 1)][0]) * (a != 0))
+        + glm_pow2((c[1] - aabb[!(b - 1)][1]) * (b != 0));
+
+  return dmin <= glm_pow2(c[2]);
+}
+
+/*!
+ * @brief check if point is inside of AABB
+ *
+ * @param[in]   aabb    bounding aabb
+ * @param[in]   point  point
+ */
+CGLM_INLINE
+bool
+glm_aabb2d_point(vec2 aabb[2], vec2 point) {
+  return (point[0] >= aabb[0][0] && point[0] <= aabb[1][0])
+      && (point[1] >= aabb[0][1] && point[1] <= aabb[1][1]);
+}
+
+/*!
+ * @brief check if AABB contains other AABB
+ *
+ * @param[in]   aabb    bounding aabb
+ * @param[in]   other  other bounding aabb
+ */
+CGLM_INLINE
+bool
+glm_aabb2d_contains(vec2 aabb[2], vec2 other[2]) {
+  return (aabb[0][0] <= other[0][0] && aabb[1][0] >= other[1][0])
+      && (aabb[0][1] <= other[0][1] && aabb[1][1] >= other[1][1]);
+}
+
+#endif /* cglm_aabb2d_h */
+

include/cglm/affine-mat.h

@@ -8,6 +8,7 @@
 /*
  Functions:
    CGLM_INLINE void glm_mul(mat4 m1, mat4 m2, mat4 dest);
+   CGLM_INLINE void glm_mul_rot(mat4 m1, mat4 m2, mat4 dest);
    CGLM_INLINE void glm_inv_tr(mat4 mat);
  */
 
@@ -30,6 +31,10 @@
 #  include "simd/neon/affine.h"
 #endif
 
+#ifdef CGLM_SIMD_WASM
+#  include "simd/wasm/affine.h"
+#endif
+
 /*!
  * @brief this is similar to glm_mat4_mul but specialized to affine transform
  *
@@ -49,7 +54,9 @@
 CGLM_INLINE
 void
 glm_mul(mat4 m1, mat4 m2, mat4 dest) {
-#ifdef __AVX__
+#if defined(__wasm__) && defined(__wasm_simd128__)
+  glm_mul_wasm(m1, m2, dest);
+#elif defined(__AVX__)
   glm_mul_avx(m1, m2, dest);
 #elif defined( __SSE__ ) || defined( __SSE2__ )
   glm_mul_sse2(m1, m2, dest);
@@ -107,7 +114,9 @@ glm_mul(mat4 m1, mat4 m2, mat4 dest) {
 CGLM_INLINE
 void
 glm_mul_rot(mat4 m1, mat4 m2, mat4 dest) {
-#if defined( __SSE__ ) || defined( __SSE2__ )
+#if defined(__wasm__) && defined(__wasm_simd128__)
+  glm_mul_rot_wasm(m1, m2, dest);
+#elif defined( __SSE__ ) || defined( __SSE2__ )
   glm_mul_rot_sse2(m1, m2, dest);
 #elif defined(CGLM_NEON_FP)
   glm_mul_rot_neon(m1, m2, dest);
@@ -156,7 +165,9 @@ glm_mul_rot(mat4 m1, mat4 m2, mat4 dest) {
 CGLM_INLINE
 void
 glm_inv_tr(mat4 mat) {
-#if defined( __SSE__ ) || defined( __SSE2__ )
+#if defined(__wasm__) && defined(__wasm_simd128__)
+  glm_inv_tr_wasm(mat);
+#elif defined( __SSE__ ) || defined( __SSE2__ )
   glm_inv_tr_sse2(mat);
 #elif defined(CGLM_NEON_FP)
   glm_inv_tr_neon(mat);

include/cglm/affine-post.h

@@ -34,9 +34,9 @@
  * @brief translate existing transform matrix by v vector
  *        and stores result in same matrix
  *
- *  this is POST transform, applies to existing transform as last transfrom
+ *  this is POST transform, applies to existing transform as last transform
  *
- * @param[in, out]  m  affine transfrom
+ * @param[in, out]  m  affine transform
  * @param[in]       v  translate vector [x, y, z]
  */
 CGLM_INLINE
@@ -51,9 +51,9 @@ glm_translated(mat4 m, vec3 v) {
  *
  * source matrix will remain same
  *
- *  this is POST transform, applies to existing transform as last transfrom
+ *  this is POST transform, applies to existing transform as last transform
  *
- * @param[in]  m    affine transfrom
+ * @param[in]  m    affine transform
  * @param[in]  v    translate vector [x, y, z]
  * @param[out] dest translated matrix
  */
@@ -67,9 +67,9 @@ glm_translated_to(mat4 m, vec3 v, mat4 dest) {
 /*!
  * @brief translate existing transform matrix by x factor
  *
- *  this is POST transform, applies to existing transform as last transfrom
+ *  this is POST transform, applies to existing transform as last transform
  *
- * @param[in, out]  m  affine transfrom
+ * @param[in, out]  m  affine transform
  * @param[in]       x  x factor
  */
 CGLM_INLINE
@@ -81,9 +81,9 @@ glm_translated_x(mat4 m, float x) {
 /*!
  * @brief translate existing transform matrix by y factor
  *
- *  this is POST transform, applies to existing transform as last transfrom
+ *  this is POST transform, applies to existing transform as last transform
  *
- * @param[in, out]  m  affine transfrom
+ * @param[in, out]  m  affine transform
  * @param[in]       y  y factor
  */
 CGLM_INLINE
@@ -95,9 +95,9 @@ glm_translated_y(mat4 m, float y) {
 /*!
  * @brief translate existing transform matrix by z factor
  *
- *  this is POST transform, applies to existing transform as last transfrom
+ *  this is POST transform, applies to existing transform as last transform
  *
- * @param[in, out]  m  affine transfrom
+ * @param[in, out]  m  affine transform
  * @param[in]       z  z factor
  */
 CGLM_INLINE
@@ -110,9 +110,9 @@ glm_translated_z(mat4 m, float z) {
  * @brief rotate existing transform matrix around X axis by angle
  *        and store result in dest
  *
- *  this is POST transform, applies to existing transform as last transfrom
+ *  this is POST transform, applies to existing transform as last transform
  *
- * @param[in]   m      affine transfrom
+ * @param[in]   m      affine transform
  * @param[in]   angle  angle (radians)
  * @param[out]  dest   rotated matrix
  */
@@ -137,9 +137,9 @@ glm_rotated_x(mat4 m, float angle, mat4 dest) {
  * @brief rotate existing transform matrix around Y axis by angle
  *        and store result in dest
  *
- *  this is POST transform, applies to existing transform as last transfrom
+ *  this is POST transform, applies to existing transform as last transform
  *
- * @param[in]   m      affine transfrom
+ * @param[in]   m      affine transform
  * @param[in]   angle  angle (radians)
  * @param[out]  dest   rotated matrix
  */
@@ -164,9 +164,9 @@ glm_rotated_y(mat4 m, float angle, mat4 dest) {
  * @brief rotate existing transform matrix around Z axis by angle
  *        and store result in dest
  *
- *  this is POST transform, applies to existing transform as last transfrom
+ *  this is POST transform, applies to existing transform as last transform
  *
- * @param[in]   m      affine transfrom
+ * @param[in]   m      affine transform
  * @param[in]   angle  angle (radians)
  * @param[out]  dest   rotated matrix
  */
@@ -190,9 +190,9 @@ glm_rotated_z(mat4 m, float angle, mat4 dest) {
 /*!
  * @brief rotate existing transform matrix around given axis by angle
  *
- *  this is POST transform, applies to existing transform as last transfrom
+ *  this is POST transform, applies to existing transform as last transform
  *
- * @param[in, out]  m      affine transfrom
+ * @param[in, out]  m      affine transform
  * @param[in]       angle  angle (radians)
  * @param[in]       axis   axis
  */
@@ -208,9 +208,9 @@ glm_rotated(mat4 m, float angle, vec3 axis) {
  * @brief rotate existing transform
  *        around given axis by angle at given pivot point (rotation center)
  *
- *  this is POST transform, applies to existing transform as last transfrom
+ *  this is POST transform, applies to existing transform as last transform
  *
- * @param[in, out]  m      affine transfrom
+ * @param[in, out]  m      affine transform
  * @param[in]       pivot  rotation center
  * @param[in]       angle  angle (radians)
  * @param[in]       axis   axis
@@ -230,9 +230,9 @@ glm_rotated_at(mat4 m, vec3 pivot, float angle, vec3 axis) {
 /*!
  * @brief rotate existing transform matrix around given axis by angle around self (doesn't affected by position)
  *
- *  this is POST transform, applies to existing transform as last transfrom
+ *  this is POST transform, applies to existing transform as last transform
  *
- * @param[in, out]  m      affine transfrom
+ * @param[in, out]  m      affine transform
  * @param[in]       angle  angle (radians)
  * @param[in]       axis   axis
  */

include/cglm/affine-pre.h

@@ -35,7 +35,7 @@
  * @brief translate existing transform matrix by v vector
  *        and stores result in same matrix
  *
- * @param[in, out]  m  affine transfrom
+ * @param[in, out]  m  affine transform
  * @param[in]       v  translate vector [x, y, z]
  */
 CGLM_INLINE
@@ -66,7 +66,7 @@ glm_translate(mat4 m, vec3 v) {
  *
  * source matrix will remain same
  *
- * @param[in]  m    affine transfrom
+ * @param[in]  m    affine transform
  * @param[in]  v    translate vector [x, y, z]
  * @param[out] dest translated matrix
  */
@@ -80,7 +80,7 @@ glm_translate_to(mat4 m, vec3 v, mat4 dest) {
 /*!
  * @brief translate existing transform matrix by x factor
  *
- * @param[in, out]  m  affine transfrom
+ * @param[in, out]  m  affine transform
  * @param[in]       x  x factor
  */
 CGLM_INLINE
@@ -98,7 +98,7 @@ glm_translate_x(mat4 m, float x) {
 /*!
  * @brief translate existing transform matrix by y factor
  *
- * @param[in, out]  m  affine transfrom
+ * @param[in, out]  m  affine transform
  * @param[in]       y  y factor
  */
 CGLM_INLINE
@@ -116,7 +116,7 @@ glm_translate_y(mat4 m, float y) {
 /*!
  * @brief translate existing transform matrix by z factor
  *
- * @param[in, out]  m  affine transfrom
+ * @param[in, out]  m  affine transform
  * @param[in]       z  z factor
  */
 CGLM_INLINE
@@ -135,7 +135,7 @@ glm_translate_z(mat4 m, float z) {
  * @brief rotate existing transform matrix around X axis by angle
  *        and store result in dest
  *
- * @param[in]   m      affine transfrom
+ * @param[in]   m      affine transform
  * @param[in]   angle  angle (radians)
  * @param[out]  dest   rotated matrix
  */
@@ -160,7 +160,7 @@ glm_rotate_x(mat4 m, float angle, mat4 dest) {
  * @brief rotate existing transform matrix around Y axis by angle
  *        and store result in dest
  *
- * @param[in]   m      affine transfrom
+ * @param[in]   m      affine transform
  * @param[in]   angle  angle (radians)
  * @param[out]  dest   rotated matrix
  */
@@ -185,7 +185,7 @@ glm_rotate_y(mat4 m, float angle, mat4 dest) {
  * @brief rotate existing transform matrix around Z axis by angle
  *        and store result in dest
  *
- * @param[in]   m      affine transfrom
+ * @param[in]   m      affine transform
  * @param[in]   angle  angle (radians)
  * @param[out]  dest   rotated matrix
  */
@@ -207,9 +207,27 @@ glm_rotate_z(mat4 m, float angle, mat4 dest) {
 }
 
 /*!
- * @brief rotate existing transform matrix around given axis by angle
+ * @brief rotate existing transform matrix 
+ *        around given axis by angle at ORIGIN (0,0,0)
  *
- * @param[in, out]  m      affine transfrom
+ *   **❗️IMPORTANT ❗️**
+ *
+ *   If you need to rotate object around itself e.g. center of object or at
+ *   some point [of object] then `glm_rotate_at()` would be better choice to do so.
+ *
+ *   Even if object's model transform is identiy, rotation may not be around
+ *   center of object if object does not lay out at ORIGIN perfectly.
+ *
+ *   Using `glm_rotate_at()` with center of bounding shape ( AABB, Sphere ... )
+ *   would be an easy option to rotate around object if object is not at origin.
+ *
+ *   One another option to rotate around itself at any point is `glm_spin()`
+ *   which is perfect if only rotating around model position is desired e.g. not
+ *   specific point on model for instance center of geometry or center of mass,
+ *   again if geometry is not perfectly centered at origin at identity transform,
+ *   rotation may not be around geometry.
+ *
+ * @param[in, out]  m      affine transform
  * @param[in]       angle  angle (radians)
  * @param[in]       axis   axis
  */
@@ -225,7 +243,7 @@ glm_rotate(mat4 m, float angle, vec3 axis) {
  * @brief rotate existing transform
  *        around given axis by angle at given pivot point (rotation center)
  *
- * @param[in, out]  m      affine transfrom
+ * @param[in, out]  m      affine transform
  * @param[in]       pivot  rotation center
  * @param[in]       angle  angle (radians)
  * @param[in]       axis   axis
@@ -250,7 +268,7 @@ glm_rotate_at(mat4 m, vec3 pivot, float angle, vec3 axis) {
  * this should work faster than glm_rotate_at because it reduces
  * one glm_translate.
  *
- * @param[out] m      affine transfrom
+ * @param[out] m      affine transform
  * @param[in]  pivot  rotation center
  * @param[in]  angle  angle (radians)
  * @param[in]  axis   axis
@@ -268,9 +286,10 @@ glm_rotate_atm(mat4 m, vec3 pivot, float angle, vec3 axis) {
 }
 
 /*!
- * @brief rotate existing transform matrix around given axis by angle around self (doesn't affected by position)
+ * @brief rotate existing transform matrix 
+ *        around given axis by angle around self (doesn't affected by position)
  *
- * @param[in, out]  m      affine transfrom
+ * @param[in, out]  m      affine transform
  * @param[in]       angle  angle (radians)
  * @param[in]       axis   axis
  */

include/cglm/affine.h

@@ -44,7 +44,7 @@
 /*!
  * @brief creates NEW translate transform matrix by v vector
  *
- * @param[out]  m  affine transfrom
+ * @param[out]  m  affine transform
  * @param[in]   v  translate vector [x, y, z]
  */
 CGLM_INLINE
@@ -58,7 +58,7 @@ glm_translate_make(mat4 m, vec3 v) {
  * @brief scale existing transform matrix by v vector
  *        and store result in dest
  *
- * @param[in]  m    affine transfrom
+ * @param[in]  m    affine transform
  * @param[in]  v    scale vector [x, y, z]
  * @param[out] dest scaled matrix
  */
@@ -75,7 +75,7 @@ glm_scale_to(mat4 m, vec3 v, mat4 dest) {
 /*!
  * @brief creates NEW scale matrix by v vector
  *
- * @param[out]  m  affine transfrom
+ * @param[out]  m  affine transform
  * @param[in]   v  scale vector [x, y, z]
  */
 CGLM_INLINE
@@ -91,7 +91,7 @@ glm_scale_make(mat4 m, vec3 v) {
  * @brief scales existing transform matrix by v vector
  *        and stores result in same matrix
  *
- * @param[in, out]  m  affine transfrom
+ * @param[in, out]  m  affine transform
  * @param[in]       v  scale vector [x, y, z]
  */
 CGLM_INLINE
@@ -104,7 +104,7 @@ glm_scale(mat4 m, vec3 v) {
  * @brief applies uniform scale to existing transform matrix v = [s, s, s]
  *        and stores result in same matrix
  *
- * @param[in, out]  m  affine transfrom
+ * @param[in, out]  m  affine transform
  * @param[in]       s  scale factor
  */
 CGLM_INLINE
@@ -119,7 +119,7 @@ glm_scale_uni(mat4 m, float s) {
  *
  * axis will be normalized so you don't need to normalize it
  *
- * @param[out] m     affine transfrom
+ * @param[out] m     affine transform
  * @param[in]  angle angle (radians)
  * @param[in]  axis  axis
  */
@@ -220,7 +220,7 @@ glm_decompose_rs(mat4 m, mat4 r, vec3 s) {
  * @brief decompose affine transform, TODO: extract shear factors.
  *        DON'T pass projected matrix here
  *
- * @param[in]  m affine transfrom
+ * @param[in]  m affine transform
  * @param[out] t translation vector
  * @param[out] r rotation matrix (mat4)
  * @param[out] s scaling vector [X, Y, Z]

include/cglm/affine2d.h

@@ -33,7 +33,7 @@
  * @brief translate existing 2d transform matrix by v vector
  *        and stores result in same matrix
  *
- * @param[in, out]  m  affine transfrom
+ * @param[in, out]  m  affine transform
  * @param[in]       v  translate vector [x, y]
  */
 CGLM_INLINE
@@ -50,7 +50,7 @@ glm_translate2d(mat3 m, vec2 v) {
  *
  * source matrix will remain same
  *
- * @param[in]  m    affine transfrom
+ * @param[in]  m    affine transform
  * @param[in]  v    translate vector [x, y]
  * @param[out] dest translated matrix
  */
@@ -64,7 +64,7 @@ glm_translate2d_to(mat3 m, vec2 v, mat3 dest) {
 /*!
  * @brief translate existing 2d transform matrix by x factor
  *
- * @param[in, out]  m  affine transfrom
+ * @param[in, out]  m  affine transform
  * @param[in]       x  x factor
  */
 CGLM_INLINE
@@ -78,7 +78,7 @@ glm_translate2d_x(mat3 m, float x) {
 /*!
  * @brief translate existing 2d transform matrix by y factor
  *
- * @param[in, out]  m  affine transfrom
+ * @param[in, out]  m  affine transform
  * @param[in]       y  y factor
  */
 CGLM_INLINE
@@ -92,7 +92,7 @@ glm_translate2d_y(mat3 m, float y) {
 /*!
  * @brief creates NEW translate 2d transform matrix by v vector
  *
- * @param[out]  m  affine transfrom
+ * @param[out]  m  affine transform
  * @param[in]   v  translate vector [x, y]
  */
 CGLM_INLINE
@@ -107,7 +107,7 @@ glm_translate2d_make(mat3 m, vec2 v) {
  * @brief scale existing 2d transform matrix by v vector
  *        and store result in dest
  *
- * @param[in]  m    affine transfrom
+ * @param[in]  m    affine transform
  * @param[in]  v    scale vector [x, y]
  * @param[out] dest scaled matrix
  */
@@ -130,7 +130,7 @@ glm_scale2d_to(mat3 m, vec2 v, mat3 dest) {
 /*!
  * @brief creates NEW 2d scale matrix by v vector
  *
- * @param[out]  m  affine transfrom
+ * @param[out]  m  affine transform
  * @param[in]   v  scale vector [x, y]
  */
 CGLM_INLINE
@@ -145,7 +145,7 @@ glm_scale2d_make(mat3 m, vec2 v) {
  * @brief scales existing 2d transform matrix by v vector
  *        and stores result in same matrix
  *
- * @param[in, out]  m  affine transfrom
+ * @param[in, out]  m  affine transform
  * @param[in]       v  scale vector [x, y]
  */
 CGLM_INLINE
@@ -164,7 +164,7 @@ glm_scale2d(mat3 m, vec2 v) {
  * @brief applies uniform scale to existing 2d transform matrix v = [s, s]
  *        and stores result in same matrix
  *
- * @param[in, out]  m  affine transfrom
+ * @param[in, out]  m  affine transform
  * @param[in]       s  scale factor
  */
 CGLM_INLINE
@@ -182,7 +182,7 @@ glm_scale2d_uni(mat3 m, float s) {
 /*!
  * @brief creates NEW rotation matrix by angle around Z axis
  *
- * @param[out] m     affine transfrom
+ * @param[out] m     affine transform
  * @param[in]  angle angle (radians)
  */
 CGLM_INLINE
@@ -210,7 +210,7 @@ glm_rotate2d_make(mat3 m, float angle) {
  * @brief rotate existing 2d transform matrix around Z axis by angle
  *         and store result in same matrix
  *
- * @param[in, out]  m      affine transfrom
+ * @param[in, out]  m      affine transform
  * @param[in]       angle  angle (radians)
  */
 CGLM_INLINE
@@ -237,7 +237,7 @@ glm_rotate2d(mat3 m, float angle) {
  * @brief rotate existing 2d transform matrix around Z axis by angle
  *        and store result in dest
  *
- * @param[in]  m      affine transfrom
+ * @param[in]  m      affine transform
  * @param[in]  angle  angle (radians)
  * @param[out] dest   destination
  */

include/cglm/bezier.h

@@ -23,7 +23,7 @@
 #define GLM_HERMITE_MAT ((mat4)GLM_HERMITE_MAT_INIT)
 
 #define CGLM_DECASTEL_EPS   1e-9f
-#define CGLM_DECASTEL_MAX   1000.0f
+#define CGLM_DECASTEL_MAX   1000
 #define CGLM_DECASTEL_SMALL 1e-20f
 
 /*!

include/cglm/box.h

@@ -104,7 +104,7 @@ glm_aabb_crop(vec3 box[2], vec3 cropBox[2], vec3 dest[2]) {
  *
  * @param[in]  box      bounding box
  * @param[in]  cropBox  crop box
- * @param[in]  clampBox miniumum box
+ * @param[in]  clampBox minimum box
  * @param[out] dest     cropped bounding box
  */
 CGLM_INLINE

include/cglm/call.h

@@ -19,14 +19,21 @@ extern "C" {
 #include "call/ivec3.h"
 #include "call/ivec4.h"
 #include "call/mat2.h"
+#include "call/mat2x3.h"
+#include "call/mat2x4.h"
 #include "call/mat3.h"
+#include "call/mat3x2.h"
+#include "call/mat3x4.h"
 #include "call/mat4.h"
+#include "call/mat4x2.h"
+#include "call/mat4x3.h"
 #include "call/affine.h"
 #include "call/cam.h"
 #include "call/quat.h"
 #include "call/euler.h"
 #include "call/plane.h"
 #include "call/frustum.h"
+#include "call/aabb2d.h"
 #include "call/box.h"
 #include "call/io.h"
 #include "call/project.h"

include/cglm/call/aabb2d.h

@@ -0,0 +1,91 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_aabb2d_h
+#define cglmc_aabb2d_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../cglm.h"
+
+/* DEPRECATED! use _diag */
+#define glmc_aabb2d_size(aabb) glmc_aabb2d_diag(aabb)
+
+CGLM_EXPORT
+void
+glmc_aabb2d_zero(vec2 aabb[2]);
+
+CGLM_EXPORT
+void
+glmc_aabb2d_copy(vec2 aabb[2], vec2 dest[2]);
+
+CGLM_EXPORT
+void
+glmc_aabb2d_transform(vec2 aabb[2], mat3 m, vec2 dest[2]);
+
+CGLM_EXPORT
+void
+glmc_aabb2d_merge(vec2 aabb1[2], vec2 aabb2[2], vec2 dest[2]);
+
+CGLM_EXPORT
+void
+glmc_aabb2d_crop(vec2 aabb[2], vec2 cropAabb[2], vec2 dest[2]);
+
+CGLM_EXPORT
+void
+glmc_aabb2d_crop_until(vec2 aabb[2],
+                     vec2 cropAabb[2],
+                     vec2 clampAabb[2],
+                     vec2 dest[2]);
+
+CGLM_EXPORT
+void
+glmc_aabb2d_invalidate(vec2 aabb[2]);
+
+CGLM_EXPORT
+bool
+glmc_aabb2d_isvalid(vec2 aabb[2]);
+
+CGLM_EXPORT
+float
+glmc_aabb2d_diag(vec2 aabb[2]);
+
+CGLM_EXPORT
+void
+glmc_aabb2d_sizev(vec2 aabb[2], vec2 dest);
+
+CGLM_EXPORT
+float
+glmc_aabb2d_radius(vec2 aabb[2]);
+
+CGLM_EXPORT
+void
+glmc_aabb2d_center(vec2 aabb[2], vec2 dest);
+
+CGLM_EXPORT
+bool
+glmc_aabb2d_aabb(vec2 aabb[2], vec2 other[2]);
+
+CGLM_EXPORT
+bool
+glmc_aabb2d_point(vec2 aabb[2], vec2 point);
+
+CGLM_EXPORT
+bool
+glmc_aabb2d_contains(vec2 aabb[2], vec2 other[2]);
+
+CGLM_EXPORT
+bool
+glmc_aabb2d_circle(vec2 aabb[2], vec3 s);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_aabb2d_h */
+
+

include/cglm/call/euler.h

@@ -49,6 +49,31 @@ CGLM_EXPORT
 void
 glmc_euler_by_order(vec3 angles, glm_euler_seq axis, mat4 dest);
 
+CGLM_EXPORT
+void
+glmc_euler_xyz_quat(vec3 angles, versor dest);
+
+CGLM_EXPORT
+void
+glmc_euler_xzy_quat(vec3 angles, versor dest);
+
+CGLM_EXPORT
+void
+glmc_euler_yxz_quat(vec3 angles, versor dest);
+
+CGLM_EXPORT
+void
+glmc_euler_yzx_quat(vec3 angles, versor dest);
+
+CGLM_EXPORT
+void
+glmc_euler_zxy_quat(vec3 angles, versor dest);
+
+CGLM_EXPORT
+void
+glmc_euler_zyx_quat(vec3 angles, versor dest);
+
+
 #ifdef __cplusplus
 }
 #endif

include/cglm/call/ivec2.h

@@ -29,6 +29,14 @@ CGLM_EXPORT
 void
 glmc_ivec2_one(ivec2 v);
 
+CGLM_EXPORT
+int
+glmc_ivec2_dot(ivec2 a, ivec2 b);
+
+CGLM_EXPORT
+int
+glmc_ivec2_cross(ivec2 a, ivec2 b);
+
 CGLM_EXPORT
 void
 glmc_ivec2_add(ivec2 a, ivec2 b, ivec2 dest);
@@ -53,6 +61,82 @@ CGLM_EXPORT
 void
 glmc_ivec2_scale(ivec2 v, int s, ivec2 dest);
 
+CGLM_EXPORT
+void
+glmc_ivec2_div(ivec2 a, ivec2 b, ivec2 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec2_divs(ivec2 v, int s, ivec2 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec2_mod(ivec2 a, ivec2 b, ivec2 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec2_addadd(ivec2 a, ivec2 b, ivec2 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec2_addadds(ivec2 a, int s, ivec2 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec2_subadd(ivec2 a, ivec2 b, ivec2 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec2_subadds(ivec2 a, int s, ivec2 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec2_muladd(ivec2 a, ivec2 b, ivec2 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec2_muladds(ivec2 a, int s, ivec2 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec2_maxadd(ivec2 a, ivec2 b, ivec2 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec2_minadd(ivec2 a, ivec2 b, ivec2 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec2_subsub(ivec2 a, ivec2 b, ivec2 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec2_subsubs(ivec2 a, int s, ivec2 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec2_addsub(ivec2 a, ivec2 b, ivec2 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec2_addsubs(ivec2 a, int s, ivec2 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec2_mulsub(ivec2 a, ivec2 b, ivec2 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec2_mulsubs(ivec2 a, int s, ivec2 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec2_maxsub(ivec2 a, ivec2 b, ivec2 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec2_minsub(ivec2 a, ivec2 b, ivec2 dest);
+
 CGLM_EXPORT
 int
 glmc_ivec2_distance2(ivec2 a, ivec2 b);
@@ -61,6 +145,18 @@ CGLM_EXPORT
 float
 glmc_ivec2_distance(ivec2 a, ivec2 b);
 
+CGLM_EXPORT
+void
+glmc_ivec2_fill(ivec2 v, int val);
+
+CGLM_EXPORT
+bool
+glmc_ivec2_eq(ivec2 v, int val);
+
+CGLM_EXPORT
+bool
+glmc_ivec2_eqv(ivec2 a, ivec2 b);
+
 CGLM_EXPORT
 void
 glmc_ivec2_maxv(ivec2 a, ivec2 b, ivec2 dest);

include/cglm/call/ivec3.h

@@ -29,6 +29,18 @@ CGLM_EXPORT
 void
 glmc_ivec3_one(ivec3 v);
 
+CGLM_EXPORT
+int
+glmc_ivec3_dot(ivec3 a, ivec3 b);
+
+CGLM_EXPORT
+int
+glmc_ivec3_norm2(ivec3 v);
+
+CGLM_EXPORT
+int
+glmc_ivec3_norm(ivec3 v);
+
 CGLM_EXPORT
 void
 glmc_ivec3_add(ivec3 a, ivec3 b, ivec3 dest);
@@ -53,6 +65,82 @@ CGLM_EXPORT
 void
 glmc_ivec3_scale(ivec3 v, int s, ivec3 dest);
 
+CGLM_EXPORT
+void
+glmc_ivec3_div(ivec3 a, ivec3 b, ivec3 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec3_divs(ivec3 v, int s, ivec3 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec3_mod(ivec3 a, ivec3 b, ivec3 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec3_addadd(ivec3 a, ivec3 b, ivec3 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec3_addadds(ivec3 a, int s, ivec3 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec3_subadd(ivec3 a, ivec3 b, ivec3 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec3_subadds(ivec3 a, int s, ivec3 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec3_muladd(ivec3 a, ivec3 b, ivec3 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec3_muladds(ivec3 a, int s, ivec3 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec3_maxadd(ivec3 a, ivec3 b, ivec3 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec3_minadd(ivec3 a, ivec3 b, ivec3 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec3_subsub(ivec3 a, ivec3 b, ivec3 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec3_subsubs(ivec3 a, int s, ivec3 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec3_addsub(ivec3 a, ivec3 b, ivec3 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec3_addsubs(ivec3 a, int s, ivec3 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec3_mulsub(ivec3 a, ivec3 b, ivec3 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec3_mulsubs(ivec3 a, int s, ivec3 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec3_maxsub(ivec3 a, ivec3 b, ivec3 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec3_minsub(ivec3 a, ivec3 b, ivec3 dest);
+
 CGLM_EXPORT
 int
 glmc_ivec3_distance2(ivec3 a, ivec3 b);
@@ -61,6 +149,18 @@ CGLM_EXPORT
 float
 glmc_ivec3_distance(ivec3 a, ivec3 b);
 
+CGLM_EXPORT
+void
+glmc_ivec3_fill(ivec3 v, int val);
+
+CGLM_EXPORT
+bool
+glmc_ivec3_eq(ivec3 v, int val);
+
+CGLM_EXPORT
+bool
+glmc_ivec3_eqv(ivec3 a, ivec3 b);
+
 CGLM_EXPORT
 void
 glmc_ivec3_maxv(ivec3 a, ivec3 b, ivec3 dest);

include/cglm/call/ivec4.h

@@ -53,6 +53,70 @@ CGLM_EXPORT
 void
 glmc_ivec4_scale(ivec4 v, int s, ivec4 dest);
 
+CGLM_EXPORT 
+void 
+glmc_ivec4_addadd(ivec4 a, ivec4 b, ivec4 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec4_addadds(ivec4 a, int s, ivec4 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec4_subadd(ivec4 a, ivec4 b, ivec4 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec4_subadds(ivec4 a, int s, ivec4 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec4_muladd(ivec4 a, ivec4 b, ivec4 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec4_muladds(ivec4 a, int s, ivec4 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec4_maxadd(ivec4 a, ivec4 b, ivec4 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec4_minadd(ivec4 a, ivec4 b, ivec4 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec4_subsub(ivec4 a, ivec4 b, ivec4 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec4_subsubs(ivec4 a, int s, ivec4 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec4_addsub(ivec4 a, ivec4 b, ivec4 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec4_addsubs(ivec4 a, int s, ivec4 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec4_mulsub(ivec4 a, ivec4 b, ivec4 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec4_mulsubs(ivec4 a, int s, ivec4 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec4_maxsub(ivec4 a, ivec4 b, ivec4 dest);
+
+CGLM_EXPORT 
+void 
+glmc_ivec4_minsub(ivec4 a, ivec4 b, ivec4 dest);
+
 CGLM_EXPORT
 int
 glmc_ivec4_distance2(ivec4 a, ivec4 b);

include/cglm/call/mat2.h

@@ -73,6 +73,10 @@ CGLM_EXPORT
 float
 glmc_mat2_rmc(vec2 r, mat2 m, vec2 c);
 
+CGLM_EXPORT
+void
+glmc_mat2_make(const float * __restrict src, mat2 dest);
+
 #ifdef __cplusplus
 }
 #endif

include/cglm/call/mat2x3.h

@@ -0,0 +1,47 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_mat2x3_h
+#define cglmc_mat2x3_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_mat2x3_copy(mat2x3 mat, mat2x3 dest);
+
+CGLM_EXPORT
+void
+glmc_mat2x3_zero(mat2x3 mat);
+
+CGLM_EXPORT
+void
+glmc_mat2x3_make(const float * __restrict src, mat2x3 dest);
+
+CGLM_EXPORT
+void
+glmc_mat2x3_mul(mat2x3 m1, mat3x2 m2, mat3 dest);
+
+CGLM_EXPORT
+void
+glmc_mat2x3_mulv(mat2x3 m, vec2 v, vec3 dest);
+
+CGLM_EXPORT
+void
+glmc_mat2x3_transpose(mat2x3 m, mat3x2 dest);
+
+CGLM_EXPORT
+void
+glmc_mat2x3_scale(mat2x3 m, float s);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_mat2x3_h */

include/cglm/call/mat2x4.h

@@ -0,0 +1,47 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_mat2x4_h
+#define cglmc_mat2x4_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_mat2x4_copy(mat2x4 mat, mat2x4 dest);
+
+CGLM_EXPORT
+void
+glmc_mat2x4_zero(mat2x4 mat);
+
+CGLM_EXPORT
+void
+glmc_mat2x4_make(const float * __restrict src, mat2x4 dest);
+
+CGLM_EXPORT
+void
+glmc_mat2x4_mul(mat2x4 m1, mat4x2 m2, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_mat2x4_mulv(mat2x4 m, vec2 v, vec4 dest);
+
+CGLM_EXPORT
+void
+glmc_mat2x4_transpose(mat2x4 m, mat4x2 dest);
+
+CGLM_EXPORT
+void
+glmc_mat2x4_scale(mat2x4 m, float s);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_mat2x4_h */

include/cglm/call/mat3.h

@@ -80,6 +80,10 @@ CGLM_EXPORT
 float
 glmc_mat3_rmc(vec3 r, mat3 m, vec3 c);
 
+CGLM_EXPORT
+void
+glmc_mat3_make(const float * __restrict src, mat3 dest);
+
 #ifdef __cplusplus
 }
 #endif

include/cglm/call/mat3x2.h

@@ -0,0 +1,47 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_mat3x2_h
+#define cglmc_mat3x2_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_mat3x2_copy(mat3x2 mat, mat3x2 dest);
+
+CGLM_EXPORT
+void
+glmc_mat3x2_zero(mat3x2 mat);
+
+CGLM_EXPORT
+void
+glmc_mat3x2_make(const float * __restrict src, mat3x2 dest);
+
+CGLM_EXPORT
+void
+glmc_mat3x2_mul(mat3x2 m1, mat2x3 m2, mat2 dest);
+
+CGLM_EXPORT
+void
+glmc_mat3x2_mulv(mat3x2 m, vec3 v, vec2 dest);
+
+CGLM_EXPORT
+void
+glmc_mat3x2_transpose(mat3x2 m, mat2x3 dest);
+
+CGLM_EXPORT
+void
+glmc_mat3x2_scale(mat3x2 m, float s);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_mat3x2_h */

include/cglm/call/mat3x4.h

@@ -0,0 +1,47 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_mat3x4_h
+#define cglmc_mat3x4_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_mat3x4_copy(mat3x4 mat, mat3x4 dest);
+
+CGLM_EXPORT
+void
+glmc_mat3x4_zero(mat3x4 mat);
+
+CGLM_EXPORT
+void
+glmc_mat3x4_make(const float * __restrict src, mat3x4 dest);
+
+CGLM_EXPORT
+void
+glmc_mat3x4_mul(mat3x4 m1, mat4x3 m2, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_mat3x4_mulv(mat3x4 m, vec3 v, vec4 dest);
+
+CGLM_EXPORT
+void
+glmc_mat3x4_transpose(mat3x4 m, mat4x3 dest);
+
+CGLM_EXPORT
+void
+glmc_mat3x4_scale(mat3x4 m, float s);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_mat3x4_h */

include/cglm/call/mat4.h

@@ -121,6 +121,10 @@ CGLM_EXPORT
 float
 glmc_mat4_rmc(vec4 r, mat4 m, vec4 c);
 
+CGLM_EXPORT
+void
+glmc_mat4_make(const float * __restrict src, mat4 dest);
+
 #ifdef __cplusplus
 }
 #endif

include/cglm/call/mat4x2.h

@@ -0,0 +1,47 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_mat4x2_h
+#define cglmc_mat4x2_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_mat4x2_copy(mat4x2 mat, mat4x2 dest);
+
+CGLM_EXPORT
+void
+glmc_mat4x2_zero(mat4x2 mat);
+
+CGLM_EXPORT
+void
+glmc_mat4x2_make(const float * __restrict src, mat4x2 dest);
+
+CGLM_EXPORT
+void
+glmc_mat4x2_mul(mat4x2 m1, mat2x4 m2, mat2 dest);
+
+CGLM_EXPORT
+void
+glmc_mat4x2_mulv(mat4x2 m, vec4 v, vec2 dest);
+
+CGLM_EXPORT
+void
+glmc_mat4x2_transpose(mat4x2 m, mat2x4 dest);
+
+CGLM_EXPORT
+void
+glmc_mat4x2_scale(mat4x2 m, float s);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_mat4x2_h */

include/cglm/call/mat4x3.h

@@ -0,0 +1,47 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_mat4x3_h
+#define cglmc_mat4x3_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_mat4x3_copy(mat4x3 mat, mat4x3 dest);
+
+CGLM_EXPORT
+void
+glmc_mat4x3_zero(mat4x3 mat);
+
+CGLM_EXPORT
+void
+glmc_mat4x3_make(const float * __restrict src, mat4x3 dest);
+
+CGLM_EXPORT
+void
+glmc_mat4x3_mul(mat4x3 m1, mat3x4 m2, mat3 dest);
+
+CGLM_EXPORT
+void
+glmc_mat4x3_mulv(mat4x3 m, vec4 v, vec3 dest);
+
+CGLM_EXPORT
+void
+glmc_mat4x3_transpose(mat4x3 m, mat3x4 dest);
+
+CGLM_EXPORT
+void
+glmc_mat4x3_scale(mat4x3 m, float s);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_mat4x3_h */

include/cglm/call/quat.h

@@ -161,6 +161,10 @@ CGLM_EXPORT
 void
 glmc_quat_rotate_atm(mat4 m, versor q, vec3 pivot);
 
+CGLM_EXPORT
+void
+glmc_quat_make(const float * __restrict src, versor dest);
+
 #ifdef __cplusplus
 }
 #endif

include/cglm/call/ray.h

@@ -20,7 +20,19 @@ glmc_ray_triangle(vec3   origin,
                   vec3   v1,
                   vec3   v2,
                   float *d);
-    
+
+CGLM_EXPORT
+bool
+glmc_ray_sphere(vec3 origin,
+                vec3 dir,
+                vec4 s,
+                float * __restrict t1,
+                float * __restrict t2);
+
+CGLM_EXPORT
+void
+glmc_ray_at(vec3 orig, vec3 dir, float t, vec3 point);
+
 #ifdef __cplusplus
 }
 #endif

include/cglm/call/vec2.h

@@ -105,6 +105,30 @@ CGLM_EXPORT
 void
 glmc_vec2_minadd(vec2 a, vec2 b, vec2 dest);
 
+CGLM_EXPORT
+void
+glmc_vec2_subsub(vec2 a, vec2 b, vec2 dest);
+
+CGLM_EXPORT
+void
+glmc_vec2_addsub(vec2 a, vec2 b, vec2 dest);
+
+CGLM_EXPORT
+void
+glmc_vec2_mulsub(vec2 a, vec2 b, vec2 dest);
+
+CGLM_EXPORT
+void
+glmc_vec2_mulsubs(vec2 a, float s, vec2 dest);
+
+CGLM_EXPORT
+void
+glmc_vec2_maxsub(vec2 a, vec2 b, vec2 dest);
+
+CGLM_EXPORT
+void
+glmc_vec2_minsub(vec2 a, vec2 b, vec2 dest);
+
 CGLM_EXPORT
 void
 glmc_vec2_negate_to(vec2 v, vec2 dest);
@@ -125,6 +149,10 @@ CGLM_EXPORT
 void
 glmc_vec2_rotate(vec2 v, float angle, vec2 dest);
 
+CGLM_EXPORT
+void
+glmc_vec2_center(vec2 a, vec2 b, vec2 dest);
+
 CGLM_EXPORT
 float
 glmc_vec2_distance2(vec2 a, vec2 b);
@@ -165,6 +193,18 @@ CGLM_EXPORT
 void
 glmc_vec2_complex_conjugate(vec2 a, vec2 dest);
 
+CGLM_EXPORT
+void
+glmc_vec2_make(const float * __restrict src, vec2 dest);
+
+CGLM_EXPORT
+void
+glmc_vec2_reflect(vec2 v, vec2 n, vec2 dest);
+
+CGLM_EXPORT
+bool
+glmc_vec2_refract(vec2 v, vec2 n, float eta, vec2 dest);
+
 #ifdef __cplusplus
 }
 #endif

include/cglm/call/vec3.h

@@ -132,6 +132,30 @@ CGLM_EXPORT
 void
 glmc_vec3_minadd(vec3 a, vec3 b, vec3 dest);
 
+CGLM_EXPORT
+void
+glmc_vec3_subsub(vec3 a, vec3 b, vec3 dest);
+
+CGLM_EXPORT
+void
+glmc_vec3_addsub(vec3 a, vec3 b, vec3 dest);
+
+CGLM_EXPORT
+void
+glmc_vec3_mulsub(vec3 a, vec3 b, vec3 dest);
+
+CGLM_EXPORT
+void
+glmc_vec3_mulsubs(vec3 a, float s, vec3 dest);
+
+CGLM_EXPORT
+void
+glmc_vec3_maxsub(vec3 a, vec3 b, vec3 dest);
+
+CGLM_EXPORT
+void
+glmc_vec3_minsub(vec3 a, vec3 b, vec3 dest);
+
 CGLM_EXPORT
 void
 glmc_vec3_negate(vec3 v);
@@ -306,6 +330,22 @@ CGLM_EXPORT
 void
 glmc_vec3_sqrt(vec3 v, vec3 dest);
 
+CGLM_EXPORT
+void
+glmc_vec3_make(const float * __restrict src, vec3 dest);
+
+CGLM_EXPORT
+void
+glmc_vec3_faceforward(vec3 n, vec3 v, vec3 nref, vec3 dest);
+
+CGLM_EXPORT
+void
+glmc_vec3_reflect(vec3 v, vec3 n, vec3 dest);
+
+CGLM_EXPORT
+bool
+glmc_vec3_refract(vec3 v, vec3 n, float eta, vec3 dest);
+
 #ifdef __cplusplus
 }
 #endif

include/cglm/call/vec4.h

@@ -133,6 +133,30 @@ CGLM_EXPORT
 void
 glmc_vec4_minadd(vec4 a, vec4 b, vec4 dest);
 
+CGLM_EXPORT
+void
+glmc_vec4_subsub(vec4 a, vec4 b, vec4 dest);
+
+CGLM_EXPORT
+void
+glmc_vec4_addsub(vec4 a, vec4 b, vec4 dest);
+
+CGLM_EXPORT
+void
+glmc_vec4_mulsub(vec4 a, vec4 b, vec4 dest);
+
+CGLM_EXPORT
+void
+glmc_vec4_mulsubs(vec4 a, float s, vec4 dest);
+
+CGLM_EXPORT
+void
+glmc_vec4_maxsub(vec4 a, vec4 b, vec4 dest);
+
+CGLM_EXPORT
+void
+glmc_vec4_minsub(vec4 a, vec4 b, vec4 dest);
+
 CGLM_EXPORT
 void
 glmc_vec4_negate(vec4 v);
@@ -283,6 +307,18 @@ CGLM_EXPORT
 void
 glmc_vec4_sqrt(vec4 v, vec4 dest);
 
+CGLM_EXPORT
+void
+glmc_vec4_make(const float * __restrict src, vec4 dest);
+
+CGLM_EXPORT
+void
+glmc_vec4_reflect(vec4 v, vec4 n, vec4 dest);
+
+CGLM_EXPORT
+bool
+glmc_vec4_refract(vec4 v, vec4 n, float eta, vec4 dest);
+
 #ifdef __cplusplus
 }
 #endif

include/cglm/cam.h

@@ -474,7 +474,7 @@ glm_persp_decomp_x(mat4 proj,
 
 /*!
  * @brief decomposes top and bottom values of perspective projection.
- *        y stands for y axis (top / botom axis)
+ *        y stands for y axis (top / bottom axis)
  *
  * @param[in]  proj   perspective projection matrix
  * @param[out] top    top

include/cglm/cglm.h

@@ -16,14 +16,21 @@
 #include "ivec3.h"
 #include "ivec4.h"
 #include "mat4.h"
+#include "mat4x2.h"
+#include "mat4x3.h"
 #include "mat3.h"
+#include "mat3x2.h"
+#include "mat3x4.h"
 #include "mat2.h"
+#include "mat2x3.h"
+#include "mat2x4.h"
 #include "affine.h"
 #include "cam.h"
 #include "frustum.h"
 #include "quat.h"
 #include "euler.h"
 #include "plane.h"
+#include "aabb2d.h"
 #include "box.h"
 #include "color.h"
 #include "util.h"

include/cglm/clipspace/persp_lh_no.h

@@ -266,7 +266,7 @@ glm_persp_decomp_x_lh_no(mat4 proj,
  * @brief decomposes top and bottom values of perspective projection
  *        with a left-hand coordinate system and a
  *        clip-space of [-1, 1].
- *        y stands for y axis (top / botom axis)
+ *        y stands for y axis (top / bottom axis)
  *
  * @param[in]  proj   perspective projection matrix
  * @param[out] top    top

include/cglm/clipspace/persp_rh_no.h

@@ -266,7 +266,7 @@ glm_persp_decomp_x_rh_no(mat4 proj,
  * @brief decomposes top and bottom values of perspective projection
  *        with a right-hand coordinate system and a
  *        clip-space of [-1, 1].
- *        y stands for y axis (top / botom axis)
+ *        y stands for y axis (top / bottom axis)
  *
  * @param[in]  proj   perspective projection matrix
  * @param[out] top    top

include/cglm/common.h

@@ -37,6 +37,18 @@
 #  define CGLM_INLINE static inline __attribute((always_inline))
 #endif
 
+#if defined(__GNUC__) || defined(__clang__)
+#  define CGLM_UNLIKELY(expr) __builtin_expect(!!(expr), 0)
+#  define CGLM_LIKELY(expr)   __builtin_expect(!!(expr), 1)
+#else
+#  define CGLM_UNLIKELY(expr) (expr)
+#  define CGLM_LIKELY(expr)   (expr)
+#endif
+
+#if defined(_M_FP_FAST) || defined(__FAST_MATH__)
+#  define CGLM_FAST_MATH
+#endif
+
 #define GLM_SHUFFLE4(z, y, x, w) (((z) << 6) | ((y) << 4) | ((x) << 2) | (w))
 #define GLM_SHUFFLE3(z, y, x)    (((z) << 4) | ((y) << 2) | (x))
 
@@ -45,7 +57,7 @@
 
 #ifndef CGLM_USE_DEFAULT_EPSILON
 #  ifndef GLM_FLT_EPSILON
-#    define GLM_FLT_EPSILON 1e-5
+#    define GLM_FLT_EPSILON 1e-5f
 #  endif
 #else
 #  define GLM_FLT_EPSILON FLT_EPSILON
@@ -82,4 +94,34 @@
 #  endif
 #endif
 
+/* struct API configurator */
+/* TODO: move struct/common.h? */
+/* WARN: dont use concant helpers outside cglm headers, because they may be changed */
+
+#define CGLM_MACRO_CONCAT_HELPER(A, B, C, D, E, ...) A ## B ## C ## D ## E ## __VA_ARGS__
+#define CGLM_MACRO_CONCAT(A, B, C, D, E, ...) CGLM_MACRO_CONCAT_HELPER(A, B, C, D, E,__VA_ARGS__)
+
+#ifndef CGLM_OMIT_NS_FROM_STRUCT_API
+#  ifndef CGLM_STRUCT_API_NS
+#    define CGLM_STRUCT_API_NS glms
+#  endif
+#  ifndef CGLM_STRUCT_API_NS_SEPERATOR
+#    define CGLM_STRUCT_API_NS_SEPERATOR _
+#  endif
+#else
+#  define CGLM_STRUCT_API_NS
+#  define CGLM_STRUCT_API_NS_SEPERATOR
+#endif
+
+#ifndef CGLM_STRUCT_API_NAME_SUFFIX
+#  define CGLM_STRUCT_API_NAME_SUFFIX
+#endif
+
+#define CGLM_STRUCTAPI(A, ...) CGLM_MACRO_CONCAT(CGLM_STRUCT_API_NS,             \
+                                                 CGLM_STRUCT_API_NS_SEPERATOR,   \
+                                                 A,                              \
+                                                 CGLM_STRUCT_API_NAME_SUFFIX,    \
+                                                 _,                              \
+                                                 __VA_ARGS__)
+
 #endif /* cglm_common_h */

include/cglm/ease.h

@@ -235,7 +235,7 @@ glm_ease_back_inout(float t) {
 
   o = 1.70158f;
   s = o * 1.525f;
-  x = 0.5;
+  x = 0.5f;
   n = t / 0.5f;
 
   if (n < 1.0f) {

include/cglm/euler.h

@@ -9,7 +9,7 @@
  NOTE:
   angles must be passed as [X-Angle, Y-Angle, Z-angle] order
   For instance you don't pass angles as [Z-Angle, X-Angle, Y-angle] to
-  glm_euler_zxy funciton, All RELATED functions accept angles same order
+  glm_euler_zxy function, All RELATED functions accept angles same order
   which is [X, Y, Z].
  */
 
@@ -30,6 +30,12 @@
    CGLM_INLINE void glm_euler_by_order(vec3         angles,
                                        glm_euler_seq ord,
                                        mat4         dest);
+   CGLM_INLINE void glm_euler_xyz_quat(vec3 angles, versor dest);
+   CGLM_INLINE void glm_euler_xzy_quat(vec3 angles, versor dest);
+   CGLM_INLINE void glm_euler_yxz_quat(vec3 angles, versor dest);
+   CGLM_INLINE void glm_euler_yzx_quat(vec3 angles, versor dest);
+   CGLM_INLINE void glm_euler_zxy_quat(vec3 angles, versor dest);
+   CGLM_INLINE void glm_euler_zyx_quat(vec3 angles, versor dest);
  */
 
 #ifndef cglm_euler_h
@@ -37,6 +43,47 @@
 
 #include "common.h"
 
+#ifdef CGLM_FORCE_LEFT_HANDED
+#  include "handed/euler_to_quat_lh.h"
+#else
+#  include "handed/euler_to_quat_rh.h"
+#endif
+
+
+#ifndef CGLM_CLIPSPACE_INCLUDE_ALL
+#  if CGLM_CONFIG_CLIP_CONTROL == CGLM_CLIP_CONTROL_LH_ZO
+#    include "clipspace/ortho_lh_zo.h"
+#    include "clipspace/persp_lh_zo.h"
+#    include "clipspace/view_lh_zo.h"
+#  elif CGLM_CONFIG_CLIP_CONTROL == CGLM_CLIP_CONTROL_LH_NO
+#    include "clipspace/ortho_lh_no.h"
+#    include "clipspace/persp_lh_no.h"
+#    include "clipspace/view_lh_no.h"
+#  elif CGLM_CONFIG_CLIP_CONTROL == CGLM_CLIP_CONTROL_RH_ZO
+#    include "clipspace/ortho_rh_zo.h"
+#    include "clipspace/persp_rh_zo.h"
+#    include "clipspace/view_rh_zo.h"
+#  elif CGLM_CONFIG_CLIP_CONTROL == CGLM_CLIP_CONTROL_RH_NO
+#    include "clipspace/ortho_rh_no.h"
+#    include "clipspace/persp_rh_no.h"
+#    include "clipspace/view_rh_no.h"
+#  endif
+#else
+#  include "clipspace/ortho_lh_zo.h"
+#  include "clipspace/persp_lh_zo.h"
+#  include "clipspace/ortho_lh_no.h"
+#  include "clipspace/persp_lh_no.h"
+#  include "clipspace/ortho_rh_zo.h"
+#  include "clipspace/persp_rh_zo.h"
+#  include "clipspace/ortho_rh_no.h"
+#  include "clipspace/persp_rh_no.h"
+#  include "clipspace/view_lh_zo.h"
+#  include "clipspace/view_lh_no.h"
+#  include "clipspace/view_rh_zo.h"
+#  include "clipspace/view_rh_no.h"
+#endif
+
+
 /*!
  * if you have axis order like vec3 orderVec = [0, 1, 2] or [0, 2, 1]...
  * vector then you can convert it to this enum by doing this:
@@ -188,7 +235,6 @@ glm_euler_xzy(vec3 angles, mat4 dest) {
   dest[3][3] =  1.0f;
 }
 
-
 /*!
  * @brief build rotation matrix from euler angles
  *
@@ -448,4 +494,108 @@ glm_euler_by_order(vec3 angles, glm_euler_seq ord, mat4 dest) {
   dest[3][3] = 1.0f;
 }
 
+
+/*!
+ * @brief creates NEW quaternion using rotation angles and does
+ *        rotations in x y z order (roll pitch yaw)
+ * 
+ * @param[in]   angles angles x y z (radians)
+ * @param[out]  dest   quaternion
+ */
+CGLM_INLINE
+void
+glm_euler_xyz_quat(vec3 angles, versor dest) {
+#ifdef CGLM_FORCE_LEFT_HANDED
+  glm_euler_xyz_quat_lh(angles, dest);
+#else
+  glm_euler_xyz_quat_rh(angles, dest);
+#endif
+}
+
+/*!
+ * @brief creates NEW quaternion using rotation angles and does
+ *        rotations in x z y order (roll yaw pitch)
+ * 
+ * @param[in]   angles angles x y z (radians)
+ * @param[out]  dest   quaternion
+ */
+CGLM_INLINE
+void
+glm_euler_xzy_quat(vec3 angles, versor dest) {
+#ifdef CGLM_FORCE_LEFT_HANDED
+  glm_euler_xzy_quat_lh(angles, dest);
+#else
+  glm_euler_xzy_quat_rh(angles, dest);
+#endif
+}
+
+/*!
+ * @brief creates NEW quaternion using rotation angles and does
+ *        rotations in y x z order (pitch roll yaw)
+ * 
+ * @param[in]   angles angles x y z (radians)
+ * @param[out]  dest   quaternion
+ */
+CGLM_INLINE
+void
+glm_euler_yxz_quat(vec3 angles, versor dest) {
+#ifdef CGLM_FORCE_LEFT_HANDED
+  glm_euler_yxz_quat_lh(angles, dest);
+#else
+  glm_euler_yxz_quat_rh(angles, dest);
+#endif
+}
+
+/*!
+ * @brief creates NEW quaternion using rotation angles and does
+ *        rotations in y z x order (pitch yaw roll)
+ * 
+ * @param[in]   angles angles x y z (radians)
+ * @param[out]  dest   quaternion
+ */
+CGLM_INLINE
+void
+glm_euler_yzx_quat(vec3 angles, versor dest) {
+#ifdef CGLM_FORCE_LEFT_HANDED
+  glm_euler_yzx_quat_lh(angles, dest);
+#else
+  glm_euler_yzx_quat_rh(angles, dest);
+#endif
+}
+
+/*!
+ * @brief creates NEW quaternion using rotation angles and does
+ *        rotations in z x y order (yaw roll pitch)
+ * 
+ * @param[in]   angles angles x y z (radians)
+ * @param[out]  dest   quaternion
+ */
+CGLM_INLINE
+void
+glm_euler_zxy_quat(vec3 angles, versor dest) {
+#ifdef CGLM_FORCE_LEFT_HANDED
+  glm_euler_zxy_quat_lh(angles, dest);
+#else
+  glm_euler_zxy_quat_rh(angles, dest);
+#endif
+}
+
+/*!
+ * @brief creates NEW quaternion using rotation angles and does
+ *        rotations in z y x order (yaw pitch roll)
+ * 
+ * @param[in]   angles angles x y z (radians)
+ * @param[out]  dest   quaternion
+ */
+CGLM_INLINE
+void
+glm_euler_zyx_quat(vec3 angles, versor dest) {
+#ifdef CGLM_FORCE_LEFT_HANDED
+  glm_euler_zyx_quat_lh(angles, dest);
+#else
+  glm_euler_zyx_quat_rh(angles, dest);
+#endif
+}
+
+
 #endif /* cglm_euler_h */

include/cglm/handed/euler_to_quat_lh.h

@@ -0,0 +1,167 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+/*
+ Functions:
+   CGLM_INLINE void glm_euler_xyz_quat_lh(vec3 angles, versor dest);
+   CGLM_INLINE void glm_euler_xzy_quat_lh(vec3 angles, versor dest);
+   CGLM_INLINE void glm_euler_yxz_quat_lh(vec3 angles, versor dest);
+   CGLM_INLINE void glm_euler_yzx_quat_lh(vec3 angles, versor dest);
+   CGLM_INLINE void glm_euler_zxy_quat_lh(vec3 angles, versor dest);
+   CGLM_INLINE void glm_euler_zyx_quat_lh(vec3 angles, versor dest);
+ */
+
+/*
+ Things to note:
+ The only difference between euler to quat rh vs lh is that the zsin part is negative
+ */
+
+#ifndef cglm_euler_to_quat_lh_h
+#define cglm_euler_to_quat_lh_h
+
+#include "../common.h"
+
+
+/*!
+ * @brief creates NEW quaternion using rotation angles and does
+ *        rotations in x y z order in left hand (roll pitch yaw)
+ * 
+ * @param[in]   angles angles x y z (radians)
+ * @param[out]  dest   quaternion
+ */
+CGLM_INLINE
+void
+glm_euler_xyz_quat_lh(vec3 angles, versor dest) {
+  float xc, yc, zc,
+        xs, ys, zs;
+
+  xs =  sinf(angles[0] * 0.5f); xc = cosf(angles[0] * 0.5f);
+  ys =  sinf(angles[1] * 0.5f); yc = cosf(angles[1] * 0.5f);
+  zs = -sinf(angles[2] * 0.5f); zc = cosf(angles[2] * 0.5f);
+
+  dest[0] = xc * ys * zs + xs * yc * zc;
+  dest[1] = xc * ys * zc - xs * yc * zs;
+  dest[2] = xc * yc * zs + xs * ys * zc;
+  dest[3] = xc * yc * zc - xs * ys * zs;
+}
+
+/*!
+ * @brief creates NEW quaternion using rotation angles and does
+ *        rotations in x z y order in left hand (roll yaw pitch)
+ * 
+ * @param[in]   angles angles x y z (radians)
+ * @param[out]  dest   quaternion
+ */
+CGLM_INLINE
+void
+glm_euler_xzy_quat_lh(vec3 angles, versor dest) {
+  float xc, yc, zc,
+        xs, ys, zs;
+
+  xs =  sinf(angles[0] * 0.5f); xc = cosf(angles[0] * 0.5f);
+  ys =  sinf(angles[1] * 0.5f); yc = cosf(angles[1] * 0.5f);
+  zs = -sinf(angles[2] * 0.5f); zc = cosf(angles[2] * 0.5f);
+
+  dest[0] = -xc * zs * ys + xs * zc * yc;
+  dest[1] =  xc * zc * ys - xs * zs * yc;
+  dest[2] =  xc * zs * yc + xs * zc * ys;
+  dest[3] =  xc * zc * yc + xs * zs * ys;  
+}
+
+/*!
+ * @brief creates NEW quaternion using rotation angles and does
+ *        rotations in y x z order in left hand (pitch roll yaw)
+ * 
+ * @param[in]   angles angles x y z (radians)
+ * @param[out]  dest   quaternion
+ */
+CGLM_INLINE
+void
+glm_euler_yxz_quat_lh(vec3 angles, versor dest) {
+  float xc, yc, zc,
+        xs, ys, zs;
+
+  xs =  sinf(angles[0] * 0.5f); xc = cosf(angles[0] * 0.5f);
+  ys =  sinf(angles[1] * 0.5f); yc = cosf(angles[1] * 0.5f);
+  zs = -sinf(angles[2] * 0.5f); zc = cosf(angles[2] * 0.5f);
+
+  dest[0] =  yc * xs * zc + ys * xc * zs;
+  dest[1] = -yc * xs * zs + ys * xc * zc;
+  dest[2] =  yc * xc * zs - ys * xs * zc;
+  dest[3] =  yc * xc * zc + ys * xs * zs;
+}
+
+/*!
+ * @brief creates NEW quaternion using rotation angles and does
+ *        rotations in y z x order in left hand (pitch yaw roll)
+ * 
+ * @param[in]   angles angles x y z (radians)
+ * @param[out]  dest   quaternion
+ */
+CGLM_INLINE
+void
+glm_euler_yzx_quat_lh(vec3 angles, versor dest) {
+  float xc, yc, zc,
+        xs, ys, zs;
+
+  xs =  sinf(angles[0] * 0.5f); xc = cosf(angles[0] * 0.5f);
+  ys =  sinf(angles[1] * 0.5f); yc = cosf(angles[1] * 0.5f);
+  zs = -sinf(angles[2] * 0.5f); zc = cosf(angles[2] * 0.5f);
+
+  dest[0] = yc * zc * xs + ys * zs * xc;
+  dest[1] = yc * zs * xs + ys * zc * xc;
+  dest[2] = yc * zs * xc - ys * zc * xs;
+  dest[3] = yc * zc * xc - ys * zs * xs;
+}
+
+/*!
+ * @brief creates NEW quaternion using rotation angles and does
+ *        rotations in z x y order in left hand (yaw roll pitch)
+ * 
+ * @param[in]   angles angles x y z (radians)
+ * @param[out]  dest   quaternion
+ */
+CGLM_INLINE
+void
+glm_euler_zxy_quat_lh(vec3 angles, versor dest) {
+  float xc, yc, zc,
+        xs, ys, zs;
+
+  xs =  sinf(angles[0] * 0.5f); xc = cosf(angles[0] * 0.5f);
+  ys =  sinf(angles[1] * 0.5f); yc = cosf(angles[1] * 0.5f);
+  zs = -sinf(angles[2] * 0.5f); zc = cosf(angles[2] * 0.5f);
+
+  dest[0] = zc * xs * yc - zs * xc * ys;
+  dest[1] = zc * xc * ys + zs * xs * yc;
+  dest[2] = zc * xs * ys + zs * xc * yc;
+  dest[3] = zc * xc * yc - zs * xs * ys;
+}
+
+/*!
+ * @brief creates NEW quaternion using rotation angles and does
+ *        rotations in z y x order in left hand (yaw pitch roll)
+ * 
+ * @param[in]   angles angles x y z (radians)
+ * @param[out]  dest   quaternion
+ */
+CGLM_INLINE
+void
+glm_euler_zyx_quat_lh(vec3 angles, versor dest) {
+  float xc, yc, zc,
+        xs, ys, zs;
+
+  xs =  sinf(angles[0] * 0.5f); xc = cosf(angles[0] * 0.5f);
+  ys =  sinf(angles[1] * 0.5f); yc = cosf(angles[1] * 0.5f);
+  zs = -sinf(angles[2] * 0.5f); zc = cosf(angles[2] * 0.5f);
+
+  dest[0] =  zc * yc * xs - zs * ys * xc;
+  dest[1] =  zc * ys * xc + zs * yc * xs;
+  dest[2] = -zc * ys * xs + zs * yc * xc;
+  dest[3] =  zc * yc * xc + zs * ys * xs;
+}
+
+#endif /*cglm_euler_to_quat_lh_h*/

include/cglm/handed/euler_to_quat_rh.h

@@ -0,0 +1,170 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+/*
+ Functions:
+   CGLM_INLINE void glm_euler_xyz_quat_rh(vec3 angles, versor dest);
+   CGLM_INLINE void glm_euler_xzy_quat_rh(vec3 angles, versor dest);
+   CGLM_INLINE void glm_euler_yxz_quat_rh(vec3 angles, versor dest);
+   CGLM_INLINE void glm_euler_yzx_quat_rh(vec3 angles, versor dest);
+   CGLM_INLINE void glm_euler_zxy_quat_rh(vec3 angles, versor dest);
+   CGLM_INLINE void glm_euler_zyx_quat_rh(vec3 angles, versor dest);
+ */
+
+/*
+ Things to note:
+ The only difference between euler to quat rh vs lh is that the zsin part is negative
+ */
+
+#ifndef cglm_euler_to_quat_rh_h
+#define cglm_euler_to_quat_rh_h
+
+#include "../common.h"
+
+/*!
+ * @brief creates NEW quaternion using rotation angles and does
+ *        rotations in x y z order in right hand (roll pitch yaw)
+ * 
+ * @param[in]   angles angles x y z (radians)
+ * @param[out]  dest   quaternion
+ */
+CGLM_INLINE
+void
+glm_euler_xyz_quat_rh(vec3 angles, versor dest) {
+  float xc, yc, zc,
+        xs, ys, zs;
+
+  xs = sinf(angles[0] * 0.5f); xc = cosf(angles[0] * 0.5f);
+  ys = sinf(angles[1] * 0.5f); yc = cosf(angles[1] * 0.5f);
+  zs = sinf(angles[2] * 0.5f); zc = cosf(angles[2] * 0.5f);
+
+  dest[0] = xc * ys * zs + xs * yc * zc;
+  dest[1] = xc * ys * zc - xs * yc * zs;
+  dest[2] = xc * yc * zs + xs * ys * zc;
+  dest[3] = xc * yc * zc - xs * ys * zs;
+
+}
+
+/*!
+ * @brief creates NEW quaternion using rotation angles and does
+ *        rotations in x z y order in right hand (roll yaw pitch)
+ * 
+ * @param[in]   angles angles x y z (radians)
+ * @param[out]  dest   quaternion
+ */
+CGLM_INLINE
+void
+glm_euler_xzy_quat_rh(vec3 angles, versor dest) {
+  float xc, yc, zc,
+        xs, ys, zs;
+
+  xs = sinf(angles[0] * 0.5f); xc = cosf(angles[0] * 0.5f);
+  ys = sinf(angles[1] * 0.5f); yc = cosf(angles[1] * 0.5f);
+  zs = sinf(angles[2] * 0.5f); zc = cosf(angles[2] * 0.5f);
+
+  dest[0] = -xc * zs * ys + xs * zc * yc;
+  dest[1] =  xc * zc * ys - xs * zs * yc;
+  dest[2] =  xc * zs * yc + xs * zc * ys;
+  dest[3] =  xc * zc * yc + xs * zs * ys;  
+
+}
+
+/*!
+ * @brief creates NEW quaternion using rotation angles and does
+ *        rotations in y x z order in right hand (pitch roll yaw)
+ * 
+ * @param[in]   angles angles x y z (radians)
+ * @param[out]  dest   quaternion
+ */
+CGLM_INLINE
+void
+glm_euler_yxz_quat_rh(vec3 angles, versor dest) {
+  float xc, yc, zc,
+        xs, ys, zs;
+
+  xs = sinf(angles[0] * 0.5f); xc = cosf(angles[0] * 0.5f);
+  ys = sinf(angles[1] * 0.5f); yc = cosf(angles[1] * 0.5f);
+  zs = sinf(angles[2] * 0.5f); zc = cosf(angles[2] * 0.5f);
+
+  dest[0] =  yc * xs * zc + ys * xc * zs;
+  dest[1] = -yc * xs * zs + ys * xc * zc;
+  dest[2] =  yc * xc * zs - ys * xs * zc;
+  dest[3] =  yc * xc * zc + ys * xs * zs;
+}
+
+/*!
+ * @brief creates NEW quaternion using rotation angles and does
+ *        rotations in y z x order in right hand (pitch yaw roll)
+ * 
+ * @param[in]   angles angles x y z (radians)
+ * @param[out]  dest   quaternion
+ */
+CGLM_INLINE
+void
+glm_euler_yzx_quat_rh(vec3 angles, versor dest) {
+  float xc, yc, zc,
+        xs, ys, zs;
+
+  xs = sinf(angles[0] * 0.5f); xc = cosf(angles[0] * 0.5f);
+  ys = sinf(angles[1] * 0.5f); yc = cosf(angles[1] * 0.5f);
+  zs = sinf(angles[2] * 0.5f); zc = cosf(angles[2] * 0.5f);
+
+  dest[0] = yc * zc * xs + ys * zs * xc;
+  dest[1] = yc * zs * xs + ys * zc * xc;
+  dest[2] = yc * zs * xc - ys * zc * xs;
+  dest[3] = yc * zc * xc - ys * zs * xs;
+
+}
+
+/*!
+ * @brief creates NEW quaternion using rotation angles and does
+ *        rotations in z x y order in right hand (yaw roll pitch)
+ * 
+ * @param[in]   angles angles x y z (radians)
+ * @param[out]  dest   quaternion
+ */
+CGLM_INLINE
+void
+glm_euler_zxy_quat_rh(vec3 angles, versor dest) {
+  float xc, yc, zc,
+        xs, ys, zs;
+
+  xs = sinf(angles[0] * 0.5f); xc = cosf(angles[0] * 0.5f);
+  ys = sinf(angles[1] * 0.5f); yc = cosf(angles[1] * 0.5f);
+  zs = sinf(angles[2] * 0.5f); zc = cosf(angles[2] * 0.5f);
+
+  dest[0] = zc * xs * yc - zs * xc * ys;
+  dest[1] = zc * xc * ys + zs * xs * yc;
+  dest[2] = zc * xs * ys + zs * xc * yc;
+  dest[3] = zc * xc * yc - zs * xs * ys;
+}
+
+/*!
+ * @brief creates NEW quaternion using rotation angles and does
+ *        rotations in z y x order in right hand (yaw pitch roll)
+ * 
+ * @param[in]   angles angles x y z (radians)
+ * @param[out]  dest   quaternion
+ */
+CGLM_INLINE
+void
+glm_euler_zyx_quat_rh(vec3 angles, versor dest) {
+  float xc, yc, zc,
+        xs, ys, zs;
+
+  xs = sinf(angles[0] * 0.5f); xc = cosf(angles[0] * 0.5f);
+  ys = sinf(angles[1] * 0.5f); yc = cosf(angles[1] * 0.5f);
+  zs = sinf(angles[2] * 0.5f); zc = cosf(angles[2] * 0.5f);
+
+  dest[0] =  zc * yc * xs - zs * ys * xc;
+  dest[1] =  zc * ys * xc + zs * yc * xs;
+  dest[2] = -zc * ys * xs + zs * yc * xc;
+  dest[3] =  zc * yc * xc + zs * ys * xs;
+}
+
+
+#endif /*cglm_euler_to_quat_rh_h*/