docs: clarify some rotation rules

Ndebug

.github/FUNDING.yml

@@ -1,7 +1,7 @@
 # These are supported funding model platforms
 
-github: # Replace with up to 4 GitHub Sponsors-enabled usernames e.g., [user1, user2]
+github: [recp]
-patreon: # Replace with a single Patreon username
+patreon: recp
 open_collective: cglm
 ko_fi: # Replace with a single Ko-fi username
 tidelift: # Replace with a single Tidelift platform-name/package-name e.g., npm/babel

.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

.gitignore

@@ -76,3 +76,7 @@ build/
 conftest.dir/*
 confdefs.h
 *.xcuserdatad
+.idea
+cmake-build-debug
+*.o.tmp
+xcode/*

.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,5 +1,10 @@
 cmake_minimum_required(VERSION 3.8.2)
-project(cglm VERSION 0.7.9 LANGUAGES C)
+project(cglm
+  VERSION 0.9.3
+  HOMEPAGE_URL https://github.com/recp/cglm
+  DESCRIPTION "OpenGL Mathematics (glm) for C"
+  LANGUAGES C
+)
 
 set(CMAKE_C_STANDARD 11)
 set(CMAKE_C_STANDARD_REQUIRED YES)
@@ -11,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)
@@ -18,25 +28,32 @@ else(CGLM_STATIC)
 endif()
 
 if(CGLM_USE_C99)
-  set(C_STANDARD 99)
+  set(CMAKE_C_STANDARD 99)
 endif()
 
 if(MSVC)
-  add_definitions(-DNDEBUG -D_WINDOWS -D_USRDLL)
+  add_definitions(-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
+  if(NOT CMAKE_BUILD_TYPE MATCHES Debug)
-  # Make sure /RTC1 is disabled, otherwise it will use functions from the CRT
+    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
@@ -47,7 +64,10 @@ include(GNUInstallDirs)
 
 set(CPACK_PROJECT_NAME ${PROJECT_NAME})
 set(CPACK_PROJECT_VERSION ${PROJECT_VERSION})
-include(CPack)
+
+if(NOT CPack_CMake_INCLUDED)
+  include(CPack)
+endif()
 
 # Target Start
 add_library(${PROJECT_NAME}
@@ -58,11 +78,20 @@ add_library(${PROJECT_NAME}
   src/quat.c
   src/cam.c
   src/vec2.c
+  src/ivec2.c
   src/vec3.c
+  src/ivec3.c
   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
@@ -73,6 +102,20 @@ add_library(${PROJECT_NAME}
   src/bezier.c
   src/ray.c
   src/affine2d.c
+  src/clipspace/ortho_lh_no.c
+  src/clipspace/ortho_lh_zo.c
+  src/clipspace/ortho_rh_no.c
+  src/clipspace/ortho_rh_zo.c
+  src/clipspace/persp_lh_no.c
+  src/clipspace/persp_lh_zo.c
+  src/clipspace/persp_rh_no.c
+  src/clipspace/persp_rh_zo.c
+  src/clipspace/view_lh_no.c
+  src/clipspace/view_lh_zo.c
+  src/clipspace/view_rh_no.c
+  src/clipspace/view_rh_zo.c
+  src/clipspace/project_no.c
+  src/clipspace/project_zo.c
   )
 
 if(CGLM_SHARED)
@@ -85,6 +128,12 @@ set_target_properties(${PROJECT_NAME} PROPERTIES
                               VERSION ${PROJECT_VERSION} 
                             SOVERSION ${PROJECT_VERSION_MAJOR})
 
+if(WIN32)
+  # Because SOVERSION has no effect to file naming on Windows
+  set_target_properties(${PROJECT_NAME} PROPERTIES
+    RUNTIME_OUTPUT_NAME ${PROJECT_NAME}-${PROJECT_VERSION_MAJOR})
+endif()
+
 target_include_directories(${PROJECT_NAME}
     PUBLIC 
         $<INSTALL_INTERFACE:include>    
@@ -93,6 +142,11 @@ target_include_directories(${PROJECT_NAME}
         ${CMAKE_CURRENT_SOURCE_DIR}/src
 )
 
+# Target for header-only usage
+add_library(${PROJECT_NAME}_headers INTERFACE)
+target_include_directories(${PROJECT_NAME}_headers INTERFACE
+  ${CMAKE_CURRENT_SOURCE_DIR}/include)
+
 # Test Configuration
 if(CGLM_USE_TEST)
   include(CTest)
@@ -105,7 +159,7 @@ install(TARGETS ${PROJECT_NAME}
         EXPORT  ${PROJECT_NAME}
         LIBRARY DESTINATION ${CMAKE_INSTALL_LIBDIR}
         ARCHIVE DESTINATION ${CMAKE_INSTALL_LIBDIR}
-        RUNTIME DESTINATION ${CMAKE_INSTALL_LIBDIR})
+        RUNTIME DESTINATION ${CMAKE_INSTALL_BINDIR})
 
 install(DIRECTORY include/${PROJECT_NAME} DESTINATION ${CMAKE_INSTALL_INCLUDEDIR}
         PATTERN ".*" EXCLUDE)
@@ -117,6 +171,25 @@ export(TARGETS ${PROJECT_NAME}
 )
 
 install(EXPORT      ${PROJECT_NAME}
+        FILE        "${PROJECT_NAME}Config.cmake"
         NAMESPACE   ${PROJECT_NAME}::
-        DESTINATION ${CMAKE_INSTALL_LIBDIR}/${PROJECT_NAME}/cmake)
+        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
+  DESTINATION ${CMAKE_INSTALL_LIBDIR}/pkgconfig)

CREDITS

@@ -74,3 +74,11 @@ Link to paper: http://webserver2.tecgraf.puc-rio.br/~mgattass/cg/trbRR/Fast%20Mi
 
 14. ARM NEON: Matrix Vector Multiplication
 https://stackoverflow.com/a/57793352/2676533
+
+16. ARM NEON Div
+
+http://github.com/microsoft/DirectXMath
+
+17. Pick Matrix
+
+glu project -> project.c

Makefile.am

@@ -11,8 +11,7 @@ AM_CFLAGS = -Wall \
             -std=gnu11 \
             -O3 \
             -Wstrict-aliasing=2 \
-            -fstrict-aliasing \
+            -fstrict-aliasing
-            -Werror=strict-prototypes
 
 lib_LTLIBRARIES = libcglm.la
 libcglm_la_LDFLAGS = -no-undefined -version-info 0:1:0
@@ -43,22 +42,34 @@ 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 \
+               include/cglm/affine-mat.h \
                include/cglm/vec2.h \
                include/cglm/vec2-ext.h \
+               include/cglm/ivec2.h \
                include/cglm/vec3.h \
                include/cglm/vec3-ext.h \
+               include/cglm/ivec3.h \
                include/cglm/vec4.h \
                include/cglm/vec4-ext.h \
+               include/cglm/ivec4.h \
                include/cglm/euler.h \
                include/cglm/util.h \
                include/cglm/quat.h \
-               include/cglm/affine-mat.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 \
@@ -69,14 +80,41 @@ cglm_HEADERS = include/cglm/version.h \
                include/cglm/ray.h \
                include/cglm/affine2d.h
 
+cglm_clipspacedir=$(includedir)/cglm/clipspace
+cglm_clipspace_HEADERS = include/cglm/clipspace/persp.h \
+                         include/cglm/clipspace/persp_lh_no.h \
+                         include/cglm/clipspace/persp_lh_zo.h \
+                         include/cglm/clipspace/persp_rh_no.h \
+                         include/cglm/clipspace/persp_rh_zo.h \
+                         include/cglm/clipspace/ortho_lh_no.h \
+                         include/cglm/clipspace/ortho_lh_zo.h \
+                         include/cglm/clipspace/ortho_rh_no.h \
+                         include/cglm/clipspace/ortho_rh_zo.h \
+                         include/cglm/clipspace/view_lh.h \
+                         include/cglm/clipspace/view_rh.h \
+                         include/cglm/clipspace/view_lh_no.h \
+                         include/cglm/clipspace/view_lh_zo.h \
+                         include/cglm/clipspace/view_rh_no.h \
+                         include/cglm/clipspace/view_rh_zo.h \
+                         include/cglm/clipspace/project_no.h \
+                         include/cglm/clipspace/project_zo.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 \
-                    include/cglm/call/affine.h \
+                    include/cglm/call/ivec2.h \
+                    include/cglm/call/ivec3.h \
+                    include/cglm/call/ivec4.h \
                     include/cglm/call/io.h \
                     include/cglm/call/cam.h \
                     include/cglm/call/quat.h \
@@ -90,12 +128,30 @@ cglm_call_HEADERS = include/cglm/call/mat4.h \
                     include/cglm/call/curve.h \
                     include/cglm/call/bezier.h \
                     include/cglm/call/ray.h \
+                    include/cglm/call/affine.h \
                     include/cglm/call/affine2d.h
 
+cglm_call_clipspacedir=$(includedir)/cglm/call/clipspace
+cglm_call_clipspace_HEADERS = include/cglm/call/clipspace/persp_lh_no.h \
+                              include/cglm/call/clipspace/persp_lh_zo.h \
+                              include/cglm/call/clipspace/persp_rh_no.h \
+                              include/cglm/call/clipspace/persp_rh_zo.h \
+                              include/cglm/call/clipspace/ortho_lh_no.h \
+                              include/cglm/call/clipspace/ortho_lh_zo.h \
+                              include/cglm/call/clipspace/ortho_rh_no.h \
+                              include/cglm/call/clipspace/ortho_rh_zo.h \
+                              include/cglm/call/clipspace/view_lh_no.h \
+                              include/cglm/call/clipspace/view_lh_zo.h \
+                              include/cglm/call/clipspace/view_rh_no.h \
+                              include/cglm/call/clipspace/view_rh_zo.h \
+                              include/cglm/call/clipspace/project_no.h \
+                              include/cglm/call/clipspace/project_zo.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 \
@@ -109,19 +165,46 @@ cglm_simd_avx_HEADERS = include/cglm/simd/avx/mat4.h \
                         include/cglm/simd/avx/affine.h
 
 cglm_simd_neondir=$(includedir)/cglm/simd/neon
-cglm_simd_neon_HEADERS = include/cglm/simd/neon/mat4.h
+cglm_simd_neon_HEADERS = include/cglm/simd/neon/affine.h \
+                         include/cglm/simd/neon/mat2.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/affine.h \
+                      include/cglm/struct/ivec4.h \
                       include/cglm/struct/io.h \
                       include/cglm/struct/cam.h \
                       include/cglm/struct/quat.h \
@@ -129,11 +212,31 @@ 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/affine2d.h
+
+cglm_struct_clipspacedir=$(includedir)/cglm/struct/clipspace
+cglm_struct_clipspace_HEADERS = include/cglm/struct/clipspace/persp_lh_no.h \
+                                include/cglm/struct/clipspace/persp_lh_zo.h \
+                                include/cglm/struct/clipspace/persp_rh_no.h \
+                                include/cglm/struct/clipspace/persp_rh_zo.h \
+                                include/cglm/struct/clipspace/ortho_lh_no.h \
+                                include/cglm/struct/clipspace/ortho_lh_zo.h \
+                                include/cglm/struct/clipspace/ortho_rh_no.h \
+                                include/cglm/struct/clipspace/ortho_rh_zo.h \
+                                include/cglm/struct/clipspace/view_lh_no.h \
+                                include/cglm/struct/clipspace/view_lh_zo.h \
+                                include/cglm/struct/clipspace/view_rh_no.h \
+                                include/cglm/struct/clipspace/view_rh_zo.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 \
@@ -142,11 +245,20 @@ libcglm_la_SOURCES=\
     src/quat.c \
     src/cam.c \
     src/vec2.c \
+    src/ivec2.c \
     src/vec3.c \
+    src/ivec3.c \
     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 \
@@ -156,13 +268,26 @@ libcglm_la_SOURCES=\
     src/curve.c \
     src/bezier.c \
     src/ray.c \
-    src/affine2d.c
+    src/affine2d.c \
+    src/clipspace/ortho_lh_no.c \
+    src/clipspace/ortho_lh_zo.c \
+    src/clipspace/ortho_rh_no.c \
+    src/clipspace/ortho_rh_zo.c \
+    src/clipspace/persp_lh_no.c \
+    src/clipspace/persp_lh_zo.c \
+    src/clipspace/persp_rh_no.c \
+    src/clipspace/persp_rh_zo.c \
+    src/clipspace/view_lh_no.c \
+    src/clipspace/view_lh_zo.c \
+    src/clipspace/view_rh_no.c \
+    src/clipspace/view_rh_zo.c \
+    src/clipspace/project_no.c \
+    src/clipspace/project_zo.c
     
 test_tests_SOURCES=\
     test/runner.c \
     test/src/test_common.c \
     test/src/tests.c \
-    test/src/test_cam.c \
     test/src/test_clamp.c \
     test/src/test_euler.c \
     test/src/test_bezier.c \

README.md

@@ -1,19 +1,55 @@
 # 🎥 OpenGL Mathematics (glm) for `C`
-[![Build Status](https://travis-ci.org/recp/cglm.svg?branch=master)](https://travis-ci.org/recp/cglm)
- [![Build status](https://ci.appveyor.com/api/projects/status/av7l3gc0yhfex8y4/branch/master?svg=true)](https://ci.appveyor.com/project/recp/cglm/branch/master)
-[![Documentation Status](https://readthedocs.org/projects/cglm/badge/?version=latest)](http://cglm.readthedocs.io/en/latest/?badge=latest)
-[![Coverage Status](https://coveralls.io/repos/github/recp/cglm/badge.svg?branch=master)](https://coveralls.io/github/recp/cglm?branch=master)
-[![codecov](https://codecov.io/gh/recp/cglm/branch/master/graph/badge.svg)](https://codecov.io/gh/recp/cglm)
-[![Codacy Badge](https://api.codacy.com/project/badge/Grade/6a62b37d5f214f178ebef269dc4a6bf1)](https://www.codacy.com/app/recp/cglm?utm_source=github.com&utm_medium=referral&utm_content=recp/cglm&utm_campaign=Badge_Grade)
-[![Backers on Open Collective](https://opencollective.com/cglm/backers/badge.svg)](#backers)
-[![Sponsors on Open Collective](https://opencollective.com/cglm/sponsors/badge.svg)](#sponsors)
 
-#### Documentation
+<p align="center">
+   <img alt="" src="cglm.png" width="550" />
+</p>
+<br>
+<p align="center">
+    <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="http://cglm.readthedocs.io/en/latest/?badge=latest">
+        <img src="https://readthedocs.org/projects/cglm/badge/?version=latest"
+             alt="Documentation Status">
+    </a>
+    <a href="https://www.codacy.com/app/recp/cglm?utm_source=github.com&amp;utm_medium=referral&amp;utm_content=recp/cglm&amp;utm_campaign=Badge_Grade">
+        <img src="https://api.codacy.com/project/badge/Grade/6a62b37d5f214f178ebef269dc4a6bf1"
+             alt="Codacy Badge"/>
+    </a>
+    <a href="https://coveralls.io/github/recp/cglm?branch=master">
+        <img src="https://coveralls.io/repos/github/recp/cglm/badge.svg?branch=master"
+             alt="Coverage Status"/>
+    </a>
+    <a href="https://codecov.io/gh/recp/cglm">
+        <img src="https://codecov.io/gh/recp/cglm/branch/master/graph/badge.svg"
+             alt="Coverage Status"/>
+    </a>
+    <br /><br />
+    <a href="#sponsors">
+        <img src="https://opencollective.com/cglm/sponsors/badge.svg"
+             alt="Sponsors on Open Collective"/>
+    </a>
+    <a href="#backers">
+        <img src="https://opencollective.com/cglm/backers/badge.svg"
+             alt="Backers on Open Collective"/>
+    </a>
+</p>
+
+<br>
+
+<p align="center">
+Highly optimized 2D|3D math library, also known as <b>OpenGL Mathematics (glm) for `C`</b>. <b>cglm</b> provides lot of utils to help math operations to be fast and quick to write. It is community friendly, feel free to bring any issues, bugs you faced. 
+</p>
+
+---
+
+#### 📚 Documentation
 
 Almost all functions (inline versions) and parameters are documented inside the corresponding headers. <br />
 Complete documentation: http://cglm.readthedocs.io
 
-#### Note for previous versions:
+#### 📌 Note for previous versions:
 
 - _dup (duplicate) is changed to _copy. For instance `glm_vec_dup -> glm_vec3_copy`
 - OpenGL related functions are dropped to make this lib platform/third-party independent
@@ -26,31 +62,29 @@ 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 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:
+#### 📌 Note for C++ developers:
 If you are not aware of the original GLM library yet, you may also want to look at:
 https://github.com/g-truc/glm
 
-#### Note for new comers (Important):
+#### 📌 Note for new comers (Important):
 - `vec4` and `mat4` variables must be aligned. (There will be unaligned versions later)
 - **in** and **[in, out]** parameters must be initialized (please). But **[out]** parameters not, initializing out param is  also redundant
 - All functions are inline if you don't want to use pre-compiled versions with glmc_ prefix, you can ignore build process. Just include headers.
 - if your debugger takes you to cglm headers then make sure you are not trying to copy vec4 to vec3 or alig issues...
 - Welcome!
 
-#### Note for experienced developers:
+#### 📌 Note for experienced developers:
 - Since I'm testing this library in my projects, sometimes bugs occurs; finding that bug[s] and making improvements would be more easy with multiple developer/contributor and their projects or knowledge. Consider to make some tests if you suspect something is wrong and any feedbacks, contributions and bug reports are always welcome.
 
-#### Allocations?
+#### 📌 Allocations?
 `cglm` doesn't alloc any memory on heap. So it doesn't provide any allocator. You should alloc memory for **out** parameters too if you pass pointer of memory location. Don't forget that **vec4** (also quat/**versor**) and **mat4** must be aligned (16-bytes), because *cglm* uses SIMD instructions to optimize most operations if available.
 
-#### Returning vector or matrix... ?
+#### 📌 Returning vector or matrix... ?
 
 **cglm** supports both *ARRAY API* and *STRUCT API*, so you can return structs if you utilize struct api (`glms_`).
 
-#### Other APIs like Vulkan, Metal, Dx?
-Currently *cglm* uses default clip space configuration (-1, 1) for camera functions (perspective, extract corners...), in the future other clip space configurations will be supported
-
 <hr/>
 
 <table>
@@ -68,7 +102,9 @@ Currently *cglm* uses default clip space configuration (-1, 1) for camera functi
   </tbody>
 </table>
 
-## Features
+## 🚀 Features
+- **scalar** and **simd** (sse, avx, neon...) 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)
@@ -95,9 +131,9 @@ Currently *cglm* uses default clip space configuration (-1, 1) for camera functi
 
 <hr />
 
-You have two option to call a function/operation: inline or library call (link)
+You have two options to call a function/operation: inline or library call (link)
 Almost all functions are marked inline (always_inline) so compiler will probably inline.
-To call pre-compiled version, just use `glmc_` (c stands for 'call') instead of `glm_`.
+To call pre-compiled versions, just use `glmc_` (c stands for 'call') instead of `glm_`.
 
 ```C
   #include <cglm/cglm.h>   /* for inline */
@@ -108,7 +144,7 @@ To call pre-compiled version, 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.
 
@@ -148,7 +184,7 @@ Struct functions generally take their parameters as *values* and *return* their
 
 The types used are actually unions that allow access to the same data multiple ways. One of those ways involves anonymous structures, available since C11. MSVC also supports it for earlier C versions out of the box and GCC/Clang do if you enable `-fms-extensions`. To explicitly enable these anonymous structures, `#define CGLM_USE_ANONYMOUS_STRUCT` to `1`, to disable them, to `0`. For backward compatibility, you can also `#define CGLM_NO_ANONYMOUS_STRUCT` (value is irrelevant) to disable them. If you don't specify explicitly, cglm will do a best guess based on your compiler and the C version you're using.
 
-## Build
+## 🔨 Build
 
 ### CMake (All platforms)
 ```bash
@@ -168,6 +204,24 @@ option(CGLM_USE_C99 "" OFF) # C11
 option(CGLM_USE_TEST "Enable Tests" OFF) # for make check - make test
 ```
 
+#### Use as header-only library with your CMake project
+
+This requires no building or installation of cglm.
+
+* Example:
+
+``` cmake
+cmake_minimum_required(VERSION 3.8.2)
+
+project(<Your Project Name>)
+
+add_executable(${PROJECT_NAME} src/main.c)
+target_link_libraries(${LIBRARY_NAME} PRIVATE
+  cglm_headers)
+
+add_subdirectory(external/cglm/ EXCLUDE_FROM_ALL)
+```
+
 #### Use with your CMake project
 * Example:
 ```cmake
@@ -184,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
@@ -199,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:
@@ -373,7 +459,7 @@ You can pass matrices the same way to other APIs e.g. Vulkan, DX...
 
 ## Notes
 
-- This library uses float types only, does not support Integers, Double... yet
+- This library does not support double type... yet
 - If headers are not working properly with your compiler, IDE please open an issue, because I'm using GCC and clang to test it maybe sometimes MSVC
 
 **TODO:**
@@ -382,7 +468,7 @@ You can pass matrices the same way to other APIs e.g. Vulkan, DX...
 - [x] Add version info
 - [ ] Unaligned operations (e.g. `glm_umat4_mul`)
 - [x] Extra documentation
-- [ ] ARM Neon Arch (In Progress)
+- [x] ARM Neon Arch
 
 
 ## Contributors

appveyor.yml

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

cglm.podspec

@@ -2,10 +2,10 @@ Pod::Spec.new do |s|
 
   # Description
   s.name         = "cglm"
-  s.version      = "0.7.2"
+  s.version      = "0.9.2"
-  s.summary      = "📽 Optimized OpenGL/Graphics Math (glm) for C"
+  s.summary      = "📽 Highly Optimized Graphics Math (glm) for C"
   s.description  = <<-DESC
-cglm is math library for graphics programming for C. It is similar to original glm but it is written for C instead of C++ (you can use here too). See the documentation or README for all features.
+cglm is math library for graphics programming for C. See the documentation or README for all features.
                    DESC
 
   s.documentation_url = "http://cglm.readthedocs.io"
@@ -25,4 +25,13 @@ cglm is math library for graphics programming for C. It is similar to original g
 
   # Linking
   s.library = "m"
+
+  # Configuration
+  s.pod_target_xcconfig = {
+    'CLANG_ENABLE_MODULES' => 'NO',
+    'CLANG_ALLOW_NON_MODULAR_INCLUDES_IN_FRAMEWORK_MODULES' => 'YES',
+    'CLANG_WARN_DOCUMENTATION_COMMENTS' => 'NO',
+    'GCC_C_LANGUAGE_STANDARD' => 'gnu11',
+    'GCC_PREPROCESSOR_DEFINITIONS' => '$(inherited) GLM_TESTS_NO_COLORFUL_OUTPUT'
+  }
 end

configure.ac

@@ -7,8 +7,8 @@
 #*****************************************************************************
 
 AC_PREREQ([2.69])
-AC_INIT([cglm], [0.7.9], [info@recp.me])
+AC_INIT([cglm], [0.9.3], [info@recp.me])
-AM_INIT_AUTOMAKE([-Wall -Werror foreign subdir-objects serial-tests])
+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

@@ -0,0 +1,129 @@
+.. default-domain:: C
+
+3D Affine Transforms (common)
+================================================================================
+
+Common transform functions.
+
+Table of contents (click to go):
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+Functions:
+
+1. :c:func:`glm_translate_make`
+#. :c:func:`glm_scale_to`
+#. :c:func:`glm_scale_make`
+#. :c:func:`glm_scale`
+#. :c:func:`glm_scale_uni`
+#. :c:func:`glm_rotate_make`
+#. :c:func:`glm_rotate_atm`
+#. :c:func:`glm_decompose_scalev`
+#. :c:func:`glm_uniscaled`
+#. :c:func:`glm_decompose_rs`
+#. :c:func:`glm_decompose`
+
+Functions documentation
+~~~~~~~~~~~~~~~~~~~~~~~
+
+.. c:function:: void  glm_translate_make(mat4 m, vec3 v)
+
+    creates NEW translate transform matrix by *v* vector.
+
+    Parameters:
+      | *[in, out]* **m**  affine transform
+      | *[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 transform
+      | *[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 transform
+      | *[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 transform
+      | *[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 transform
+      | *[in]*      **v** scale factor
+
+.. 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 transform
+      | *[in]*  **axis** 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 transform
+      | *[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 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

@@ -0,0 +1,127 @@
+.. default-domain:: C
+
+3D Affine Transforms (post)
+================================================================================
+
+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):
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+Functions:
+
+1. :c:func:`glm_translated_to`
+#. :c:func:`glm_translated`
+#. :c:func:`glm_translated_x`
+#. :c:func:`glm_translated_y`
+#. :c:func:`glm_translated_z`
+#. :c:func:`glm_rotated_x`
+#. :c:func:`glm_rotated_y`
+#. :c:func:`glm_rotated_z`
+#. :c:func:`glm_rotated`
+#. :c:func:`glm_rotated_at`
+#. :c:func:`glm_spinned`
+
+Functions documentation
+~~~~~~~~~~~~~~~~~~~~~~~
+
+.. c:function:: void  glm_translated_to(mat4 m, vec3 v, mat4 dest)
+
+    translate existing transform matrix by *v* vector and store result in dest
+
+    Parameters:
+      | *[in]*  **m**    affine transform
+      | *[in]*  **v**    translate vector [x, y, z]
+      | *[out]* **dest** translated matrix
+
+.. c:function:: void  glm_translated(mat4 m, vec3 v)
+
+    translate existing transform matrix by *v* vector
+    and stores result in same matrix
+
+    Parameters:
+      | *[in, out]* **m**  affine transform
+      | *[in]*      **v**  translate vector [x, y, z]
+
+.. c:function:: void  glm_translated_x(mat4 m, float x)
+
+    translate existing transform matrix by x factor
+
+    Parameters:
+      | *[in, out]* **m**  affine transform
+      | *[in]*      **v**  x factor
+
+.. c:function:: void  glm_translated_y(mat4 m, float y)
+
+    translate existing transform matrix by *y* factor
+
+    Parameters:
+      | *[in, out]* **m**  affine transform
+      | *[in]*      **v**  y factor
+
+.. c:function:: void  glm_translated_z(mat4 m, float z)
+
+    translate existing transform matrix by *z* factor
+
+    Parameters:
+      | *[in, out]* **m**  affine transform
+      | *[in]*      **v**  z factor
+
+.. c:function:: void  glm_rotated_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 transform
+      | *[in]*  **angle** angle (radians)
+      | *[out]* **dest**  rotated matrix
+
+.. c:function:: void  glm_rotated_y(mat4 m, float angle, mat4 dest)
+
+    rotate existing transform matrix around Y axis by angle
+    and store result in dest
+
+    Parameters:
+      | *[in]*  **m**     affine transform
+      | *[in]*  **angle** angle (radians)
+      | *[out]* **dest**  rotated matrix
+
+.. c:function:: void  glm_rotated_z(mat4 m, float angle, mat4 dest)
+
+    rotate existing transform matrix around Z axis by angle
+    and store result in dest
+
+    Parameters:
+      | *[in]*  **m**     affine transform
+      | *[in]*  **angle** angle (radians)
+      | *[out]* **dest**  rotated matrix
+
+.. c:function:: void  glm_rotated(mat4 m, float angle, vec3 axis)
+
+    rotate existing transform matrix around Z axis by angle and axis
+
+    Parameters:
+      | *[in, out]* **m**     affine transform
+      | *[in]*      **angle** angle (radians)
+      | *[in]*      **axis**  axis
+
+.. c:function:: void  glm_rotated_at(mat4 m, vec3 pivot, float angle, vec3 axis)
+
+    rotate existing transform around given axis by angle at given pivot point (rotation center)
+
+    Parameters:
+      | *[in, out]* **m**     affine transform
+      | *[in]*      **pivot** pivot, anchor point, rotation center
+      | *[in]*      **angle** angle (radians)
+      | *[in]*      **axis**  axis
+
+.. c:function:: void  glm_spinned(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 transform
+      | *[in]*      **angle** angle (radians)
+      | *[in]*      **axis**  axis

docs/source/affine-pre.rst

@@ -0,0 +1,154 @@
+.. default-domain:: C
+
+3D Affine Transforms (pre)
+================================================================================
+
+Pre transform functions which are regular transform functions.
+
+Table of contents (click to go):
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+Functions:
+
+1. :c:func:`glm_translate_to`
+#. :c:func:`glm_translate`
+#. :c:func:`glm_translate_x`
+#. :c:func:`glm_translate_y`
+#. :c:func:`glm_translate_z`
+#. :c:func:`glm_translate_make`
+#. :c:func:`glm_scale_to`
+#. :c:func:`glm_scale_make`
+#. :c:func:`glm_scale`
+#. :c:func:`glm_scale_uni`
+#. :c:func:`glm_rotate_x`
+#. :c:func:`glm_rotate_y`
+#. :c:func:`glm_rotate_z`
+#. :c:func:`glm_rotate_make`
+#. :c:func:`glm_rotate`
+#. :c:func:`glm_rotate_at`
+#. :c:func:`glm_rotate_atm`
+#. :c:func:`glm_decompose_scalev`
+#. :c:func:`glm_uniscaled`
+#. :c:func:`glm_decompose_rs`
+#. :c:func:`glm_decompose`
+#. :c:func:`glm_spin`
+
+Functions documentation
+~~~~~~~~~~~~~~~~~~~~~~~
+
+.. c:function:: void  glm_translate_to(mat4 m, vec3 v, mat4 dest)
+
+    translate existing transform matrix by *v* vector and store result in dest
+
+    Parameters:
+      | *[in]*  **m**    affine transform
+      | *[in]*  **v**    translate vector [x, y, z]
+      | *[out]* **dest** translated matrix
+
+.. c:function:: void  glm_translate(mat4 m, vec3 v)
+
+    translate existing transform matrix by *v* vector
+    and stores result in same matrix
+
+    Parameters:
+      | *[in, out]* **m**  affine transform
+      | *[in]*      **v**  translate vector [x, y, z]
+
+.. c:function:: void  glm_translate_x(mat4 m, float x)
+
+    translate existing transform matrix by x factor
+
+    Parameters:
+      | *[in, out]* **m**  affine transform
+      | *[in]*      **v**  x factor
+
+.. c:function:: void  glm_translate_y(mat4 m, float y)
+
+    translate existing transform matrix by *y* factor
+
+    Parameters:
+      | *[in, out]* **m**  affine transform
+      | *[in]*      **v**  y factor
+
+.. c:function:: void  glm_translate_z(mat4 m, float z)
+
+    translate existing transform matrix by *z* factor
+
+    Parameters:
+      | *[in, out]* **m**  affine transform
+      | *[in]*      **v**  z 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 transform
+      | *[in]*  **angle** angle (radians)
+      | *[out]* **dest**  rotated matrix
+
+.. c:function:: void  glm_rotate_y(mat4 m, float angle, mat4 dest)
+
+    rotate existing transform matrix around Y axis by angle
+    and store result in dest
+
+    Parameters:
+      | *[in]*  **m**     affine transform
+      | *[in]*  **angle** angle (radians)
+      | *[out]* **dest**  rotated matrix
+
+.. c:function:: void  glm_rotate_z(mat4 m, float angle, mat4 dest)
+
+    rotate existing transform matrix around Z axis by angle
+    and store result in dest
+
+    Parameters:
+      | *[in]*  **m**     affine transform
+      | *[in]*  **angle** angle (radians)
+      | *[out]* **dest**  rotated matrix
+
+.. c:function:: void  glm_rotate(mat4 m, float angle, vec3 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 transform
+      | *[in]*      **angle** angle (radians)
+      | *[in]*      **axis**  axis
+
+.. c:function:: void  glm_rotate_at(mat4 m, vec3 pivot, float angle, vec3 axis)
+
+    rotate existing transform around given axis by angle at given pivot point (rotation center)
+
+    Parameters:
+      | *[in, out]* **m**     affine transform
+      | *[in]*      **pivot** pivot, anchor point, rotation center
+      | *[in]*      **angle** angle (radians)
+      | *[in]*      **axis**  axis
+
+.. 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 transform
+      | *[in]*      **angle** angle (radians)
+      | *[in]*      **axis**  axis

docs/source/affine.rst

@@ -5,6 +5,18 @@
 
 Header: cglm/affine.h
 
+Before starting, **cglm** provides two kind of transform functions; pre and post. 
+
+Pre functions (`T' = Tnew * T`) are like `glm_translate`, `glm_rotate` which means it will translate the vector first and then apply the model transformation.
+Post functions (`T' = T * Tnew`) are like `glm_translated`, `glm_rotated` which means it will apply the model transformation first and then translate the vector.
+
+`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 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.
+
 Initialize Transform Matrices
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 Functions with **_make** prefix expect you don't have a matrix and they create
@@ -12,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.
@@ -23,18 +35,22 @@ 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. 
 
 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)
 3. Move model back from origin to pivot (reverse of step-1): **translate(pivot.x, pivot.y, pivot.z)**
 
 **glm_rotate_at**, **glm_quat_rotate_at** and their helper functions works that way. 
+So if you use them you don't need to do these steps manually which are done by **cglm**.
 
 The implementation would be:
 
@@ -45,6 +61,13 @@ The implementation would be:
   glm_rotate(m, angle, axis);
   glm_translate(m, pivotInv); /* pivotInv = -pivot */
 
+or just:
+
+.. code-block:: c
+  :linenos:
+
+  glm_rotate_at(m, pivot, angle, axis);
+
 .. _TransformsOrder:
 
 Transforms Order
@@ -54,16 +77,16 @@ It is important to understand this part especially if you call transform
 functions multiple times
 
 `glm_translate`, `glm_rotate`, `glm_scale` and `glm_quat_rotate` and their
-helpers functions works like this (cglm may provide reverse order too as alternative in the future):
+helpers functions works like this (cglm provides reverse order as `ed` suffix e.g `glm_translated`, `glm_rotated` see post transforms):
 
 .. 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:
@@ -147,199 +170,27 @@ Functions:
 #. :c:func:`glm_decompose_rs`
 #. :c:func:`glm_decompose`
 
+Post functions (**NEW**):
+
+1. :c:func:`glm_translated_to`
+#. :c:func:`glm_translated`
+#. :c:func:`glm_translated_x`
+#. :c:func:`glm_translated_y`
+#. :c:func:`glm_translated_z`
+#. :c:func:`glm_rotated_x`
+#. :c:func:`glm_rotated_y`
+#. :c:func:`glm_rotated_z`
+#. :c:func:`glm_rotated`
+#. :c:func:`glm_rotated_at`
+#. :c:func:`glm_spinned`
+
 Functions documentation
 ~~~~~~~~~~~~~~~~~~~~~~~
 
-.. c:function:: void  glm_translate_to(mat4 m, vec3 v, mat4 dest)
+.. toctree::
+   :maxdepth: 1
+   :caption: Affine categories:
 
-    translate existing transform matrix by *v* vector and store result in dest
+   affine-common
-
+   affine-pre
-    Parameters:
+   affine-post
-      | *[in]*  **m**    affine transfrom
-      | *[in]*  **v**    translate vector [x, y, z]
-      | *[out]* **dest** translated matrix
-
-.. c:function:: void  glm_translate(mat4 m, vec3 v)
-
-    translate existing transform matrix by *v* vector
-    and stores result in same matrix
-
-    Parameters:
-      | *[in, out]* **m**  affine transfrom
-      | *[in]*      **v**  translate vector [x, y, z]
-
-.. c:function:: void  glm_translate_x(mat4 m, float x)
-
-    translate existing transform matrix by x factor
-
-    Parameters:
-      | *[in, out]* **m**  affine transfrom
-      | *[in]*      **v**  x factor
-
-.. c:function:: void  glm_translate_y(mat4 m, float y)
-
-    translate existing transform matrix by *y* factor
-
-    Parameters:
-      | *[in, out]* **m**  affine transfrom
-      | *[in]*      **v**  y factor
-
-.. c:function:: void  glm_translate_z(mat4 m, float z)
-
-    translate existing transform matrix by *z* factor
-
-    Parameters:
-      | *[in, out]* **m**  affine transfrom
-      | *[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]*  **angle** angle (radians)
-      | *[out]* **dest**  rotated matrix
-
-.. c:function:: void  glm_rotate_y(mat4 m, float angle, mat4 dest)
-
-    rotate existing transform matrix around Y axis by angle
-    and store result in dest
-
-    Parameters:
-      | *[in]*  **m**     affine transfrom
-      | *[in]*  **angle** angle (radians)
-      | *[out]* **dest**  rotated matrix
-
-.. c:function:: void  glm_rotate_z(mat4 m, float angle, mat4 dest)
-
-    rotate existing transform matrix around Z axis by angle
-    and store result in dest
-
-    Parameters:
-      | *[in]*  **m**     affine transfrom
-      | *[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
-
-    Parameters:
-      | *[in, out]* **m**     affine transfrom
-      | *[in]*      **angle** angle (radians)
-      | *[in]*      **axis**  axis
-
-.. c:function:: void  glm_rotate_at(mat4 m, vec3 pivot, float angle, vec3 axis)
-
-    rotate existing transform around given axis by angle at given pivot point (rotation center)
-
-    Parameters:
-      | *[in, out]* **m**     affine transfrom
-      | *[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]

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,56 +1,28 @@
-API documentation
+📚 API documentation
 ================================
 
-Some functions may exist twice,
+**cglm** provides a few APIs for similar functions.
-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
+* 📦 **Inline API**: All functions are inline. You can include **cglm/cglm.h** header
-of two **vec3**. Now we can also do this with :code:`glm_dot`,
+  to use this API. This is the default API. `glm_` is namespace/prefix for this API. 
-same for *_cross* and so on...
+* 📦 **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
+And also there are also sub categories:
-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
+* 📦 **Array API**: Types are raw arrays and functions takes array as argument. You can include **cglm/cglm.h** header
-of functions. Since functions in this header[s] are inline you don't need to
+  to use this API. This is the default API. `glm_` is namespace/prefix for this API. 
-build or link *cglm* against your project.
+* 📦 **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*
+📌 Since struct api and call api are built top of inline array api, follow inline array api docs for individual functions. 
-version of functions. You need to build *cglm* and link it.
-Follow the :doc:`build` documentation for this
 
 .. toctree::
    :maxdepth: 1
-   :caption: API categories:
+   :caption: API documentations:
 
-   affine
+   api_inline_array
-   affine-mat
+   api_struct
-   affine2d
+   api_call
-   cam
+   api_simd
-   frustum
-   box
-   quat
-   euler
-   mat4
-   mat3
-   mat2
-   vec3
-   vec3-ext
-   vec4
-   vec4-ext
-   vec2
-   vec2-ext
-   color
-   plane
-   project
-   util
-   io
-   call
-   sphere
-   curve
-   bezier
-   version
-   ray

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,98 @@
+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
+
+
+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:**
+.. 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.
+   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
+   Just import cglm to your project as dependency / external lib by copy-paste then use it as usual
 
 CMake (All platforms):
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
@@ -32,6 +32,22 @@ If you don't want to install **cglm** to your system's folder you can get static
   option(CGLM_USE_C99 "" OFF) # C11 
   option(CGLM_USE_TEST "Enable Tests" OFF) # for make check - make test
 
+**Use as header-only library with your CMake project example**
+This requires no building or installation of cglm.
+
+.. code-block:: CMake
+  :linenos:
+
+  cmake_minimum_required(VERSION 3.8.2)
+  
+  project(<Your Project Name>)
+  
+  add_executable(${PROJECT_NAME} src/main.c)
+  target_link_libraries(${LIBRARY_NAME} PRIVATE
+    cglm_headers)
+  
+  add_subdirectory(external/cglm/ EXCLUDE_FROM_ALL)
+
 **Use with your CMake project example**
 
 .. code-block:: CMake
@@ -50,7 +66,7 @@ If you don't want to install **cglm** to your system's folder you can get static
 Meson (All platforms):
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
-.. code-block:: bash
+.. code-block::
   :linenos:
 
   $ meson build # [Optional] --default-library=static
@@ -60,7 +76,7 @@ Meson (All platforms):
 
 **Meson Options:**
 
-.. code-block:: CMake
+.. code-block:: 
   :linenos:
 
   c_std=c11
@@ -69,9 +85,9 @@ Meson (All platforms):
   enable_tests=false # to run tests: ninja test
 
 
-**Use with your CMake project example**
+**Use with your Meson project**
 
-.. code-block:: CMake
+.. code-block::
   :linenos:
 
   # Clone cglm or create a cglm.wrap under <source_root>/subprojects

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,7 +18,7 @@ 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
+.. 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
@@ -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
@@ -140,7 +140,7 @@ Functions documentation
     | set up perspective projection matrix
 
     Parameters:
-      | *[in]*  **fovy**    field of view angle
+      | *[in]*  **fovy**    field of view angle (in radians)
       | *[in]*  **aspect**  aspect ratio ( width / height )
       | *[in]*  **nearVal** near clipping plane
       | *[in]*  **farVal**  far clipping planes
@@ -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

@@ -25,7 +25,7 @@
 
 # If your documentation needs a minimal Sphinx version, state it here.
 #
-# needs_sphinx = '1.0'
+# needs_sphinx = '3.0'
 
 # Add any Sphinx extension module names here, as strings. They can be
 # extensions coming with Sphinx (named 'sphinx.ext.*') or your custom
@@ -62,16 +62,16 @@ author = u'Recep Aslantas'
 # built documents.
 #
 # The short X.Y version.
-version = u'0.7.9'
+version = u'0.9.2'
 # The full version, including alpha/beta/rc tags.
-release = u'0.7.9'
+release = u'0.9.2'
 
 # 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.
@@ -111,7 +111,7 @@ 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 = ['_static']
 
 
 # -- Options for HTMLHelp output ------------------------------------------
@@ -197,3 +197,7 @@ epub_exclude_files = ['search.html']
 
 # If true, `todo` and `todoList` produce output, else they produce nothing.
 todo_include_todos = True
+
+# -- Options for the C domain ------------------------------------------------
+
+c_id_attributes = ['__restrict']

docs/source/features.rst

@@ -1,6 +1,8 @@
 Features
 ================================================================================
 
+* **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)
@@ -16,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.
+.. 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.
+   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**
+   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

@@ -6,13 +6,12 @@
 cglm Documentation
 ================================
 
-**cglm** is optimized 3D math library written in C99 (compatible with C89).
+**cglm** is an optimized 3D math library written in C99 (compatible with C89).
-It is similar to original **glm** library except this is mainly for **C**
+It is similar to the original **glm** library, except **cglm** is mainly for
+**C**.
 
-This library stores matrices as column-major order but in the future row-major
+**cglm** stores matrices as column-major order but in the future row-major is
-is considered to be supported as optional.
+considered to be supported as optional.
-
-Also currently only **float** type is supported for most operations.
 
 .. toctree::
    :maxdepth: 2
@@ -29,7 +28,7 @@ Also currently only **float** type is supported for most operations.
    opengl
 
 .. toctree::
-   :maxdepth: 2
+   :maxdepth: 3
    :caption: API:
 
    api
@@ -46,8 +45,8 @@ Also currently only **float** type is supported for most operations.
 
    troubleshooting
 
-Indices and tables
+Indices and Tables:
-==================
+===================
 
 * :ref:`genindex`
 * :ref:`modindex`

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
+.. note:: print functions use **%0.4f** precision if you need more
-(you probably will in some cases), you can change it temporary.
+   (you probably will in some cases), you can change it temporary.
-cglm may provide precision parameter in the future
+   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

@@ -0,0 +1,260 @@
+.. default-domain:: C
+
+ivec2
+=====
+
+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`
+#. :c:func:`glm_ivec2_minv`
+#. :c:func:`glm_ivec2_clamp`
+#. :c:func:`glm_ivec2_abs`
+
+Functions documentation
+~~~~~~~~~~~~~~~~~~~~~~~
+
+.. c:function:: void glm_ivec2(int * v, ivec2 dest)
+    
+    init ivec2 using vec3 or vec4
+    
+    Parameters:
+      | *[in]*  **v**    vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec2_copy(ivec2 a, ivec2 dest)
+    
+    copy all members of [a] to [dest]
+
+    Parameters:
+      | *[in]*  **a**    source vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec2_zero(ivec2 v)
+
+    set all members of [v] to zero
+
+    Parameters:
+      | *[out]* **v** vector
+
+.. c:function:: void glm_ivec2_one(ivec2 v)
+
+    set all members of [v] to one
+
+    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]
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec2_adds(ivec2 v, int s, ivec2 dest)
+
+    add scalar s to vector [v] and store result in [dest]
+    
+    Parameters:
+      | *[in]*  **v**    vector
+      | *[in]*  **s**    scalar
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec2_sub(ivec2 a, ivec2 b, ivec2 dest)
+
+    subtract vector [b] from vector [a] and store result in [dest]
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec2_subs(ivec2 v, int s, ivec2 dest)
+
+    subtract scalar s from vector [v] and store result in [dest]
+    
+    Parameters:
+      | *[in]*  **v**    vector
+      | *[in]*  **s**    scalar
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec2_mul(ivec2 a, ivec2 b, ivec2 dest)
+
+    multiply vector [a] with vector [b] and store result in [dest]
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec2_scale(ivec2 v, int s, ivec2 dest)
+
+    multiply vector [a] with scalar s and store result in [dest]
+    
+    Parameters:
+      | *[in]*  **v**    vector
+      | *[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
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+    
+    Returns:
+        squared distance (distance * distance)
+
+.. c:function:: float glm_ivec2_distance(ivec2 a, ivec2 b)
+
+    distance between two vectors
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+    
+    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
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec2_minv(ivec2 a, ivec2 b, ivec2 dest)
+
+    set each member of dest to lesser of vector a and b
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec2_clamp(ivec2 v, int minVal, int maxVal)
+
+    clamp each member of [v] between minVal and maxVal (inclusive)
+
+    Parameters:
+      | *[in, out]* **v**      vector
+      | *[in]*      **minVal** minimum value
+      | *[in]*      **maxVal** maximum value
+
+.. c:function:: void glm_ivec2_abs(ivec2 v, ivec2 dest)
+
+    absolute value of each vector item
+
+    Parameters:
+      | *[in]*   **v**     vector
+      | *[out]*  **dest**  destination vector

docs/source/ivec3.rst

@@ -0,0 +1,272 @@
+.. default-domain:: C
+
+ivec3
+=====
+
+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`
+#. :c:func:`glm_ivec2_abs`
+
+Functions documentation
+~~~~~~~~~~~~~~~~~~~~~~~
+
+.. c:function:: void glm_ivec3(ivec4 v4, ivec3 dest)
+    
+    init ivec3 using ivec4
+    
+    Parameters:
+      | *[in]*  **v**    vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec3_copy(ivec3 a, ivec3 dest)
+    
+    copy all members of [a] to [dest]
+
+    Parameters:
+      | *[in]*  **a**    source vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec3_zero(ivec3 v)
+
+    set all members of [v] to zero
+
+    Parameters:
+      | *[out]* **v** vector
+
+.. c:function:: void glm_ivec3_one(ivec3 v)
+
+    set all members of [v] to one
+
+    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]
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec3_adds(ivec3 v, int s, ivec3 dest)
+
+    add scalar s to vector [v] and store result in [dest]
+    
+    Parameters:
+      | *[in]*  **v**    vector
+      | *[in]*  **s**    scalar
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec3_sub(ivec3 a, ivec3 b, ivec3 dest)
+
+    subtract vector [b] from vector [a] and store result in [dest]
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec3_subs(ivec3 v, int s, ivec3 dest)
+
+    subtract scalar s from vector [v] and store result in [dest]
+    
+    Parameters:
+      | *[in]*  **v**    vector
+      | *[in]*  **s**    scalar
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec3_mul(ivec3 a, ivec3 b, ivec3 dest)
+
+    multiply vector [a] with vector [b] and store result in [dest]
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec3_scale(ivec3 v, int s, ivec3 dest)
+
+    multiply vector [a] with scalar s and store result in [dest]
+    
+    Parameters:
+      | *[in]*  **v**    vector
+      | *[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
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+    
+    Returns:
+        squared distance (distance * distance)
+
+.. c:function:: float glm_ivec3_distance(ivec3 a, ivec3 b)
+
+    distance between two vectors
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+    
+    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
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec3_minv(ivec3 a, ivec3 b, ivec3 dest)
+
+    set each member of dest to lesser of vector a and b
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec3_clamp(ivec3 v, int minVal, int maxVal)
+
+    clamp each member of [v] between minVal and maxVal (inclusive)
+
+    Parameters:
+      | *[in, out]* **v**      vector
+      | *[in]*      **minVal** minimum value
+      | *[in]*      **maxVal** maximum value
+
+.. c:function:: void glm_ivec3_abs(ivec3 v, ivec3 dest)
+
+    absolute value of each vector item
+
+    Parameters:
+      | *[in]*   **v**     vector
+      | *[out]*  **dest**  destination vector

docs/source/ivec4.rst

@@ -0,0 +1,179 @@
+.. default-domain:: C
+
+ivec4
+=====
+
+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`
+#. :c:func:`glm_ivec4_copy`
+#. :c:func:`glm_ivec4_zero`
+#. :c:func:`glm_ivec4_one`
+#. :c:func:`glm_ivec4_add`
+#. :c:func:`glm_ivec4_adds`
+#. :c:func:`glm_ivec4_sub`
+#. :c:func:`glm_ivec4_subs`
+#. :c:func:`glm_ivec4_mul`
+#. :c:func:`glm_ivec4_scale`
+#. :c:func:`glm_ivec4_distance2`
+#. :c:func:`glm_ivec4_distance`
+#. :c:func:`glm_ivec4_maxv`
+#. :c:func:`glm_ivec4_minv`
+#. :c:func:`glm_ivec4_clamp`
+#. :c:func:`glm_ivec4_abs`
+
+Functions documentation
+~~~~~~~~~~~~~~~~~~~~~~~
+
+.. c:function:: void glm_ivec4(ivec3 v3, int last, ivec4 dest)
+    
+    init ivec4 using ivec3
+    
+    Parameters:
+      | *[in]*  **v**    vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec4_copy(ivec4 a, ivec4 dest)
+    
+    copy all members of [a] to [dest]
+
+    Parameters:
+      | *[in]*  **a**    source vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec4_zero(ivec4 v)
+
+    set all members of [v] to zero
+
+    Parameters:
+      | *[out]* **v** vector
+
+.. c:function:: void glm_ivec4_one(ivec4 v)
+
+    set all members of [v] to one
+
+    Parameters:
+      | *[out]* **v** vector
+
+.. c:function:: void glm_ivec4_add(ivec4 a, ivec4 b, ivec4 dest)
+
+    add vector [a] to vector [b] and store result in [dest]
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec4_adds(ivec4 v, int s, ivec4 dest)
+
+    add scalar s to vector [v] and store result in [dest]
+    
+    Parameters:
+      | *[in]*  **v**    vector
+      | *[in]*  **s**    scalar
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec4_sub(ivec4 a, ivec4 b, ivec4 dest)
+
+    subtract vector [b] from vector [a] and store result in [dest]
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec4_subs(ivec4 v, int s, ivec4 dest)
+
+    subtract scalar s from vector [v] and store result in [dest]
+    
+    Parameters:
+      | *[in]*  **v**    vector
+      | *[in]*  **s**    scalar
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec4_mul(ivec4 a, ivec4 b, ivec4 dest)
+
+    multiply vector [a] with vector [b] and store result in [dest]
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec4_scale(ivec4 v, int s, ivec4 dest)
+
+    multiply vector [a] with scalar s and store result in [dest]
+    
+    Parameters:
+      | *[in]*  **v**    vector
+      | *[in]*  **s**    scalar
+      | *[out]* **dest** destination
+
+.. c:function:: int glm_ivec4_distance2(ivec4 a, ivec4 b)
+
+    squared distance between two vectors
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+    
+    Returns:
+        squared distance (distance * distance)
+
+.. c:function:: float glm_ivec4_distance(ivec4 a, ivec4 b)
+
+    distance between two vectors
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+    
+    Returns:
+        distance
+
+.. c:function:: void glm_ivec4_maxv(ivec4 a, ivec4 b, ivec4 dest)
+
+    set each member of dest to greater of vector a and b
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec4_minv(ivec4 a, ivec4 b, ivec4 dest)
+
+    set each member of dest to lesser of vector a and b
+
+    Parameters:
+      | *[in]*  **a**    first vector
+      | *[in]*  **b**    second vector
+      | *[out]* **dest** destination
+
+.. c:function:: void glm_ivec4_clamp(ivec4 v, int minVal, int maxVal)
+
+    clamp each member of [v] between minVal and maxVal (inclusive)
+
+    Parameters:
+      | *[in, out]* **v**      vector
+      | *[in]*      **minVal** minimum value
+      | *[in]*      **maxVal** maximum value
+
+.. c:function:: void glm_ivec4_abs(ivec4 v, ivec4 dest)
+
+    absolute value of each vector item
+
+    Parameters:
+      | *[in]*   **v**     vector
+      | *[out]*  **dest**  destination vector

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, out]* **m**  matrix
-      | *[in]*      **dest**  scalar
+      | *[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(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,90 @@
+.. 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`
+
+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(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, mat2 dest)
+
+    multiply m1 and m2 to dest
+
+    .. code-block:: c
+
+       glm_mat2x3_mul(mat2x3, mat3x2, mat2);
+
+    Parameters:
+      | *[in]*  **m1**    left matrix (mat2x3)
+      | *[in]*  **m2**    right matrix (mat3x2)
+      | *[out]* **dest**  destination matrix (mat2)
+
+.. c:function:: void glm_mat2x3_mulv(mat2x3 m, vec3 v, vec2 dest)
+
+    multiply mat2x3 with vec3 (column vector) and store in dest vector
+
+    Parameters:
+      | *[in]*  **m**     mat2x3 (left)
+      | *[in]*  **v**     vec3 (right, column vector)
+      | *[out]* **dest**  destination (result, column vector)
+
+.. 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,90 @@
+.. 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`
+
+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(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, mat2 dest)
+
+    multiply m1 and m2 to dest
+
+    .. code-block:: c
+
+       glm_mat2x4_mul(mat2x4, mat4x2, mat2);
+
+    Parameters:
+      | *[in]*  **m1**    left matrix (mat2x4)
+      | *[in]*  **m2**    right matrix (mat4x2)
+      | *[out]* **dest**  destination matrix (mat2)
+
+.. c:function:: void glm_mat2x4_mulv(mat2x4 m, vec4 v, vec2 dest)
+
+    multiply mat2x4 with vec4 (column vector) and store in dest vector
+
+    Parameters:
+      | *[in]*  **m**     mat2x4 (left)
+      | *[in]*  **v**     vec4 (right, column vector)
+      | *[out]* **dest**  destination (result, column vector)
+
+.. 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, out]* **m** matrix
-      | *[in]*      **dest**  scalar
+      | *[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(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,90 @@
+.. 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`
+
+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(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, mat3 dest)
+
+    multiply m1 and m2 to dest
+
+    .. code-block:: c
+
+       glm_mat3x2_mul(mat3x2, mat2x3, mat3);
+
+    Parameters:
+      | *[in]*  **m1**    left matrix (mat3x2)
+      | *[in]*  **m2**    right matrix (mat2x3)
+      | *[out]* **dest**  destination matrix (mat3)
+
+.. c:function:: void glm_mat3x2_mulv(mat3x2 m, vec2 v, vec3 dest)
+
+    multiply mat3x2 with vec2 (column vector) and store in dest vector
+
+    Parameters:
+      | *[in]*  **m**     mat3x2 (left)
+      | *[in]*  **v**     vec3 (right, column vector)
+      | *[out]* **dest**  destination (result, column vector)
+
+.. 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,90 @@
+.. 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`
+
+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(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, mat3 dest)
+
+    multiply m1 and m2 to dest
+
+    .. code-block:: c
+
+       glm_mat3x4_mul(mat3x4, mat4x3, mat3);
+
+    Parameters:
+      | *[in]*  **m1**    left matrix (mat3x4)
+      | *[in]*  **m2**    right matrix (mat4x3)
+      | *[out]* **dest**  destination matrix (mat3)
+
+.. c:function:: void glm_mat3x4_mulv(mat3x4 m, vec4 v, vec3 dest)
+
+    multiply mat3x4 with vec4 (column vector) and store in dest vector
+
+    Parameters:
+      | *[in]*  **m**     mat3x4 (left)
+      | *[in]*  **v**     vec4 (right, column vector)
+      | *[out]* **dest**  destination (result, column vector)
+
+.. 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_mat4_trace`
-#. :c:func:`glm_mat3_trace3`
+#. :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(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,90 @@
+.. 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`
+
+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(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, mat4 dest)
+
+    multiply m1 and m2 to dest
+
+    .. code-block:: c
+
+       glm_mat4x2_mul(mat4x2, mat2x4, mat4);
+
+    Parameters:
+      | *[in]*  **m1**    left matrix (mat4x2)
+      | *[in]*  **m2**    right matrix (mat2x4)
+      | *[out]* **dest**  destination matrix (mat4)
+
+.. c:function:: void glm_mat4x2_mulv(mat4x2 m, vec2 v, vec4 dest)
+
+    multiply mat4x2 with vec2 (column vector) and store in dest vector
+
+    Parameters:
+      | *[in]*  **m**     mat4x2 (left)
+      | *[in]*  **v**     vec2 (right, column vector)
+      | *[out]* **dest**  destination (result, column vector)
+
+.. 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,90 @@
+.. 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`
+
+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(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, mat4 dest)
+
+    multiply m1 and m2 to dest
+
+    .. code-block:: c
+
+       glm_mat4x3_mul(mat4x3, mat3x4, mat4);
+
+    Parameters:
+      | *[in]*  **m1**    left matrix (mat4x3)
+      | *[in]*  **m2**    right matrix (mat3x4)
+      | *[out]* **dest**  destination matrix (mat4)
+
+.. c:function:: void glm_mat4x3_mulv(mat4x3 m, vec3 v, vec4 dest)
+
+    multiply mat4x3 with vec3 (column vector) and store in dest vector
+
+    Parameters:
+      | *[in]*  **m**     mat4x3 (left)
+      | *[in]*  **v**     vec3 (right, column vector)
+      | *[out]* **dest**  destination (result, column vector)
+
+.. 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

@@ -2,7 +2,7 @@ How to send vector or matrix to OpenGL like API
 ==================================================
 
 *cglm*'s vector and matrix types are arrays. So you can send them directly to a
-function which accecpts pointer. But you may got warnings for matrix because it is
+function which accepts pointer. But you may got warnings for matrix because it is
 two dimensional array.
 
 Passing / Uniforming Matrix to OpenGL:
@@ -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,6 +1,6 @@
 .. default-domain:: C
 
-Options
+🛠️ Options
 ===============================================================================
 
 A few options are provided via macros.
@@ -29,12 +29,51 @@ have to compile cglm with **CGLM_ALL_UNALIGNED** macro.
  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**.
+ cglm functions will read from and write to **INVALID MEMORY LOCATIONSNs**.
 
  ALWAYS USE SAME CONFIGURATION / OPTION for **cglm** if you have multiple projects.
 
  For instance if you set CGLM_ALL_UNALIGNED in a project then set it in other projects too
 
+Clipspace Option[s]
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+By starting **v0.8.3** cglm provides options to switch between clipspace configurations.
+
+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 include all headers for you.
+
+1. **CGLM_CLIPSPACE_INCLUDE_ALL**
+2. **CGLM_FORCE_DEPTH_ZERO_TO_ONE**
+3. **CGLM_FORCE_LEFT_HANDED**
+
+
+1. **CGLM_CLIPSPACE_INCLUDE_ALL**:
+
+By defining this macro, **cglm** will include all clipspace functions for you by just using
+`#include cglm/cglm.h` or `#include cglm/struct.h` or `#include cglm/call.h`
+
+Otherwise you need to include header you want manually e.g. `#include cglm/clipspace/view_rh_zo.h`
+
+2. **CGLM_FORCE_DEPTH_ZERO_TO_ONE**
+
+This is similar to **GLM**'s **GLM_FORCE_DEPTH_ZERO_TO_ONE** option. 
+This will set clip space between 0 to 1 which makes **cglm** Vulkan, Metal friendly. 
+
+You can use functions like `glm_lookat_lh_zo()` individually. By setting **CGLM_FORCE_DEPTH_ZERO_TO_ONE**
+functions in cam.h for instance will use `_zo` versions.
+
+3. **CGLM_FORCE_LEFT_HANDED**
+
+Force **cglm** to use the left handed coordinate system by default, currently **cglm** uses right handed coordinate system as default,
+you can change this behavior with this option.
+
+**VERY VERY IMPORTANT:**
+
+Be careful if you include **cglm** in multiple projects.
+
 SSE and SSE2 Shuffle Option
 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 **_mm_shuffle_ps** generates **shufps** instruction even if registers are same.
@@ -51,6 +90,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/project.rst

@@ -27,8 +27,8 @@ Functions documentation
       the matrix should contain projection matrix.
 
     if you don't have ( and don't want to have ) an inverse matrix then use
-    glm_unproject version. You may use existing inverse of matrix in somewhere
+    :c:func:`glm_unproject` version. You may use existing inverse of matrix in somewhere
-    else, this is why glm_unprojecti exists to save save inversion cost
+    else, this is why **glm_unprojecti** exists to save inversion cost
 
     [1] space:
       - if m = invProj:     View Space
@@ -57,7 +57,7 @@ Functions documentation
     | maps the specified viewport coordinates into specified space [1]
       the matrix should contain projection matrix.
 
-    this is same as glm_unprojecti except this function get inverse matrix for
+    this is same as :c:func:`glm_unprojecti` except this function get inverse matrix for
     you.
 
     [1] space:
@@ -91,12 +91,29 @@ Functions documentation
        glm_mat4_mul(proj, view, viewProj);
        glm_mat4_mul(viewProj, model, MVP);
 
-    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]*  **pos**      object coordinates
       | *[in]*  **m**        MVP matrix
       | *[in]*  **vp**       viewport as [x, y, width, height]
       | *[out]* **dest**     projected coordinates
+
+.. c:function:: float glm_project_z(vec3 pos, mat4 m)
+
+    | map object's z coordinate to window coordinates
+
+    this is same as :c:func:`glm_project` except this function projects only Z coordinate
+    which reduces a few calculations and parameters.
+
+    Computing MVP:
+
+    .. code-block:: c
+
+       glm_mat4_mul(proj, view, viewProj);
+       glm_mat4_mul(viewProj, model, MVP);
+
+    Parameters:
+      | *[in]*  **pos**      object coordinates
+      | *[in]*  **m**        MVP matrix
+
+    Returns:
+        projected z coordinate

docs/source/quat.rst

@@ -32,6 +32,7 @@ Functions:
 #. :c:func:`glm_quat`
 #. :c:func:`glm_quatv`
 #. :c:func:`glm_quat_copy`
+#. :c:func:`glm_quat_from_vecs`
 #. :c:func:`glm_quat_norm`
 #. :c:func:`glm_quat_normalize`
 #. :c:func:`glm_quat_normalize_to`
@@ -52,6 +53,7 @@ Functions:
 #. :c:func:`glm_quat_mat3`
 #. :c:func:`glm_quat_mat3t`
 #. :c:func:`glm_quat_lerp`
+#. :c:func:`glm_quat_nlerp`
 #. :c:func:`glm_quat_slerp`
 #. :c:func:`glm_quat_look`
 #. :c:func:`glm_quat_for`
@@ -60,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
 ~~~~~~~~~~~~~~~~~~~~~~~
@@ -122,6 +125,20 @@ Functions documentation
       | *[in]*  **q**     source quaternion
       | *[out]* **dest**  destination quaternion
 
+.. c:function:: void  glm_quat_from_vecs(vec3 a, vec3 b, versor dest)
+
+    | compute unit quaternion needed to rotate a into b
+
+    References:
+      * `Finding quaternion representing the rotation from one vector to another <https://stackoverflow.com/a/11741520/183120>`_
+      * `Quaternion from two vectors <http://lolengine.net/blog/2014/02/24/quaternion-from-two-vectors-final>`_
+      * `Angle between vectors <http://www.euclideanspace.com/maths/algebra/vectors/angleBetween/minorlogic.htm>`_
+
+    Parameters:
+      | *[in]*  **a**     unit vector
+      | *[in]*  **b**     unit vector
+      | *[in]*  **dest**  unit quaternion
+
 .. c:function:: float  glm_quat_norm(versor q)
 
     | returns norm (magnitude) of quaternion
@@ -304,6 +321,25 @@ Functions documentation
       | *[in]*  **t**     interpolant (amount) clamped between 0 and 1
       | *[out]* **dest**  result quaternion
 
+.. c:function:: void glm_quat_nlerp(versor q, versor r, float  t, versor dest)
+
+    | interpolates between two quaternions
+    | taking the shortest rotation path using
+    | normalized linear interpolation (NLERP)
+
+    | This is a cheaper alternative to slerp; most games use nlerp
+    | for animations as it visually makes little difference.
+
+    References:
+      * `Understanding Slerp, Then Not Using it <http://number-none.com/product/Understanding%20Slerp,%20Then%20Not%20Using%20It>`_
+      * `Lerp, Slerp and Nlerp <https://keithmaggio.wordpress.com/2011/02/15/math-magician-lerp-slerp-and-nlerp/>`_
+
+    Parameters:
+      | *[in]*  **from**  from
+      | *[in]*  **to**    to
+      | *[in]*  **t**     interpolant (amount) clamped between 0 and 1
+      | *[out]* **dest**  result quaternion
+
 .. c:function:: void glm_quat_slerp(versor q, versor r, float  t, versor dest)
 
     | interpolates between two quaternions
@@ -374,7 +410,7 @@ Functions documentation
       | *[in]*       **q**      quaternion
       | *[in]*       **pivot**  pivot
 
-.. c:function:: void glm_quat_rotate(mat4 m, versor q, mat4 dest)
+.. c:function:: void glm_quat_rotate_atm(mat4 m, versor q, vec3 pivot)
 
     | rotate NEW transform matrix using quaternion at pivot point
     | this creates rotation matrix, it assumes you don't have a matrix
@@ -385,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(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/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/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.
+Recall that **cglm** does not alloc any memory on the heap.
-cglm functions works like memcpy; it copies data from src,
+cglm functions work like memcpy; they copy data from src,
-makes calculations then copy the result to dest.
+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.
+**vec4** and **mat4** types require 16 byte alignment.
-These types are marked with align attribute to let compiler know about this
+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.
+**UPDATE:** Starting with v0.4.5, cglm provides an option to disable the 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.
+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.
+For example, you may have called a **glm_vec4_** function for a **vec3** data type.
-It will try to write 32 byte but since **vec3** is 24 byte it should throw
+It will try to write 32 bytes, but since **vec3** is 24 bytes, it should throw
-memory access error or exit the app without saying anything.
+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.
+  | If you send GLM_VEC4_ONE or similar macros directly to a function, it may crash.
-  | Because compiler may not apply alignment as defined on **typedef** to that macro while passing it (on stack) to a function.
+  | 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**,
+For instance if you use **glm_normalize()** or **glm_vec3_normalize()** for a **vec4**,
-it will assume that passed param is **vec3** and will normalize it for **vec3**.
+it will assume that the passed param is a **vec3**, and will normalize it for **vec3**.
-Since you need to **vec4** to be normalized in your case, you will get wrong results.
+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-ext.rst

@@ -26,6 +26,7 @@ Functions:
 #. :c:func:`glm_vec2_isinf`
 #. :c:func:`glm_vec2_isvalid`
 #. :c:func:`glm_vec2_sign`
+#. :c:func:`glm_vec2_abs`
 #. :c:func:`glm_vec2_sqrt`
 
 Functions documentation
@@ -125,6 +126,14 @@ Functions documentation
       | *[in]*   **v**     vector
       | *[out]*  **dest**  sign vector (only keeps signs as -1, 0, -1)
 
+.. c:function:: void glm_vec2_abs(vec2 v, vec2 dest)
+
+    absolute value of each vector item
+
+    Parameters:
+      | *[in]*   **v**     vector
+      | *[out]*  **dest**  destination vector
+
 .. c:function:: void glm_vec2_sqrt(vec2 v, vec2 dest)
 
     square root of each vector item

docs/source/vec2.rst

@@ -10,10 +10,10 @@ Table of contents (click to go):
 
 Macros:
 
-1. GLM_vec2_ONE_INIT
+1. GLM_VEC2_ONE_INIT
-#. GLM_vec2_ZERO_INIT
+#. GLM_VEC2_ZERO_INIT
-#. GLM_vec2_ONE
+#. GLM_VEC2_ONE
-#. GLM_vec2_ZERO
+#. GLM_VEC2_ZERO
 
 Functions:
 
@@ -45,17 +45,19 @@ 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`
 
 Functions documentation
 ~~~~~~~~~~~~~~~~~~~~~~~
 
-.. c:function:: void glm_vec2(float \*v, vec2 dest)
+.. c:function:: void glm_vec2(float * v, vec2 dest)
 
     init vec2 using vec3 or vec4
 
@@ -115,7 +117,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 +315,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 +384,13 @@ Functions documentation
       | *[in]*  **to**     to value
       | *[in]*  **t**      interpolant (amount) clamped between 0 and 1
       | *[out]* **dest**   destination
+
+.. c:function:: void glm_vec2_make(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

docs/source/vec3-ext.rst

@@ -27,6 +27,7 @@ Functions:
 #. :c:func:`glm_vec3_isinf`
 #. :c:func:`glm_vec3_isvalid`
 #. :c:func:`glm_vec3_sign`
+#. :c:func:`glm_vec3_abs`
 #. :c:func:`glm_vec3_sqrt`
 
 Functions documentation
@@ -49,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)
@@ -134,6 +143,14 @@ Functions documentation
       | *[in]*   **v**     vector
       | *[out]*  **dest**  sign vector (only keeps signs as -1, 0, -1)
 
+.. c:function:: void glm_vec3_abs(vec3 v, vec3 dest)
+
+    absolute value of each vector item
+
+    Parameters:
+      | *[in]*   **v**     vector
+      | *[out]*  **dest**  destination vector
+
 .. c:function:: void glm_vec3_sqrt(vec3 v, vec3 dest)
 
     square root of each vector item

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,7 @@ Functions:
 #. :c:func:`glm_vec3_ortho`
 #. :c:func:`glm_vec3_clamp`
 #. :c:func:`glm_vec3_lerp`
+#. :c:func:`glm_vec3_make`
 
 Functions documentation
 ~~~~~~~~~~~~~~~~~~~~~~~
@@ -147,7 +148,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 +309,7 @@ Functions documentation
 
 .. c:function:: void  glm_vec3_flipsign(vec3 v)
 
-    **DEPRACATED!**
+    **DEPRECATED!**
 
     use :c:func:`glm_vec3_negate`
 
@@ -317,7 +318,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 +328,7 @@ Functions documentation
 
 .. c:function:: void  glm_vec3_inv(vec3 v)
 
-    **DEPRACATED!**
+    **DEPRECATED!**
 
     use :c:func:`glm_vec3_negate`
 
@@ -336,7 +337,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 +377,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
@@ -474,6 +475,9 @@ Functions documentation
 
     possible orthogonal/perpendicular vector
 
+    References:
+      * `On picking an orthogonal vector (and combing coconuts) <http://lolengine.net/blog/2013/09/21/picking-orthogonal-vector-combing-coconuts>`_
+
     Parameters:
       | *[in]*  **v**     vector
       | *[out]* **dest**  orthogonal/perpendicular vector
@@ -498,3 +502,13 @@ Functions documentation
       | *[in]*  **to**     to value
       | *[in]*  **t**      interpolant (amount) clamped between 0 and 1
       | *[out]* **dest**   destination
+
+.. c:function:: void glm_vec3_make(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

docs/source/vec4.rst

@@ -59,6 +59,7 @@ Functions:
 #. :c:func:`glm_vec4_clamp`
 #. :c:func:`glm_vec4_lerp`
 #. :c:func:`glm_vec4_cubic`
+#. :c:func:`glm_vec4_make`
 
 Functions documentation
 ~~~~~~~~~~~~~~~~~~~~~~~
@@ -107,6 +108,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 +131,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 +291,7 @@ Functions documentation
 
 .. c:function:: void  glm_vec4_flipsign(vec4 v)
 
-    **DEPRACATED!**
+    **DEPRECATED!**
 
     use :c:func:`glm_vec4_negate`
 
@@ -292,7 +300,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 +310,7 @@ Functions documentation
 
 .. c:function:: void  glm_vec4_inv(vec4 v)
 
-    **DEPRACATED!**
+    **DEPRECATED!**
 
     use :c:func:`glm_vec4_negate`
 
@@ -311,7 +319,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 +414,13 @@ Functions documentation
     Parameters:
       | *[in]*  **s**      parameter
       | *[out]* **dest**   destination
+
+.. c:function:: void glm_vec4_make(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

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);
  */
 
@@ -26,6 +27,14 @@
 #  include "simd/avx/affine.h"
 #endif
 
+#ifdef CGLM_NEON_FP
+#  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
  *
@@ -45,10 +54,14 @@
 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);
+#elif defined(CGLM_NEON_FP)
+  glm_mul_neon(m1, m2, dest);
 #else
   float a00 = m1[0][0], a01 = m1[0][1], a02 = m1[0][2], a03 = m1[0][3],
         a10 = m1[1][0], a11 = m1[1][1], a12 = m1[1][2], a13 = m1[1][3],
@@ -101,8 +114,12 @@ 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);
 #else
   float a00 = m1[0][0], a01 = m1[0][1], a02 = m1[0][2], a03 = m1[0][3],
         a10 = m1[1][0], a11 = m1[1][1], a12 = m1[1][2], a13 = m1[1][3],
@@ -148,8 +165,12 @@ 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);
 #else
   CGLM_ALIGN_MAT mat3 r;
   CGLM_ALIGN(8)  vec3 t;

include/cglm/affine-post.h

@@ -0,0 +1,247 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglm_affine_post_h
+#define cglm_affine_post_h
+
+/*
+ Functions:
+   CGLM_INLINE void glm_translated_to(mat4 m, vec3 v, mat4 dest);
+   CGLM_INLINE void glm_translated(mat4 m, vec3 v);
+   CGLM_INLINE void glm_translated_x(mat4 m, float to);
+   CGLM_INLINE void glm_translated_y(mat4 m, float to);
+   CGLM_INLINE void glm_translated_z(mat4 m, float to);
+   CGLM_INLINE void glm_rotated_x(mat4 m, float angle, mat4 dest);
+   CGLM_INLINE void glm_rotated_y(mat4 m, float angle, mat4 dest);
+   CGLM_INLINE void glm_rotated_z(mat4 m, float angle, mat4 dest);
+   CGLM_INLINE void glm_rotated(mat4 m, float angle, vec3 axis);
+   CGLM_INLINE void glm_rotated_at(mat4 m, vec3 pivot, float angle, vec3 axis);
+   CGLM_INLINE void glm_spinned(mat4 m, float angle, vec3 axis);
+ */
+
+#include "common.h"
+#include "util.h"
+#include "vec3.h"
+#include "vec4.h"
+#include "mat4.h"
+#include "affine-mat.h"
+
+/*!
+ * @brief translate existing transform matrix by v vector
+ *        and stores result in same matrix
+ *
+ *  this is POST transform, applies to existing transform as last transform
+ *
+ * @param[in, out]  m  affine transform
+ * @param[in]       v  translate vector [x, y, z]
+ */
+CGLM_INLINE
+void
+glm_translated(mat4 m, vec3 v) {
+  glm_vec3_add(m[3], v, m[3]);
+}
+
+/*!
+ * @brief translate existing transform matrix by v vector
+ *        and store result in dest
+ *
+ * source matrix will remain same
+ *
+ *  this is POST transform, applies to existing transform as last transform
+ *
+ * @param[in]  m    affine transform
+ * @param[in]  v    translate vector [x, y, z]
+ * @param[out] dest translated matrix
+ */
+CGLM_INLINE
+void
+glm_translated_to(mat4 m, vec3 v, mat4 dest) {
+  glm_mat4_copy(m, dest);
+  glm_translated(dest, v);
+}
+
+/*!
+ * @brief translate existing transform matrix by x factor
+ *
+ *  this is POST transform, applies to existing transform as last transform
+ *
+ * @param[in, out]  m  affine transform
+ * @param[in]       x  x factor
+ */
+CGLM_INLINE
+void
+glm_translated_x(mat4 m, float x) {
+  m[3][0] += x;
+}
+
+/*!
+ * @brief translate existing transform matrix by y factor
+ *
+ *  this is POST transform, applies to existing transform as last transform
+ *
+ * @param[in, out]  m  affine transform
+ * @param[in]       y  y factor
+ */
+CGLM_INLINE
+void
+glm_translated_y(mat4 m, float y) {
+  m[3][1] += y;
+}
+
+/*!
+ * @brief translate existing transform matrix by z factor
+ *
+ *  this is POST transform, applies to existing transform as last transform
+ *
+ * @param[in, out]  m  affine transform
+ * @param[in]       z  z factor
+ */
+CGLM_INLINE
+void
+glm_translated_z(mat4 m, float z) {
+  m[3][2] += 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 transform
+ *
+ * @param[in]   m      affine transform
+ * @param[in]   angle  angle (radians)
+ * @param[out]  dest   rotated matrix
+ */
+CGLM_INLINE
+void
+glm_rotated_x(mat4 m, float angle, mat4 dest) {
+  CGLM_ALIGN_MAT mat4 t = GLM_MAT4_IDENTITY_INIT;
+  float c, s;
+
+  c = cosf(angle);
+  s = sinf(angle);
+
+  t[1][1] =  c;
+  t[1][2] =  s;
+  t[2][1] = -s;
+  t[2][2] =  c;
+
+  glm_mul_rot(t, m, 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 transform
+ *
+ * @param[in]   m      affine transform
+ * @param[in]   angle  angle (radians)
+ * @param[out]  dest   rotated matrix
+ */
+CGLM_INLINE
+void
+glm_rotated_y(mat4 m, float angle, mat4 dest) {
+  CGLM_ALIGN_MAT mat4 t = GLM_MAT4_IDENTITY_INIT;
+  float c, s;
+
+  c = cosf(angle);
+  s = sinf(angle);
+
+  t[0][0] =  c;
+  t[0][2] = -s;
+  t[2][0] =  s;
+  t[2][2] =  c;
+
+  glm_mul_rot(t, m, 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 transform
+ *
+ * @param[in]   m      affine transform
+ * @param[in]   angle  angle (radians)
+ * @param[out]  dest   rotated matrix
+ */
+CGLM_INLINE
+void
+glm_rotated_z(mat4 m, float angle, mat4 dest) {
+  CGLM_ALIGN_MAT mat4 t = GLM_MAT4_IDENTITY_INIT;
+  float c, s;
+
+  c = cosf(angle);
+  s = sinf(angle);
+
+  t[0][0] =  c;
+  t[0][1] =  s;
+  t[1][0] = -s;
+  t[1][1] =  c;
+
+  glm_mul_rot(t, m, dest);
+}
+
+/*!
+ * @brief rotate existing transform matrix around given axis by angle
+ *
+ *  this is POST transform, applies to existing transform as last transform
+ *
+ * @param[in, out]  m      affine transform
+ * @param[in]       angle  angle (radians)
+ * @param[in]       axis   axis
+ */
+CGLM_INLINE
+void
+glm_rotated(mat4 m, float angle, vec3 axis) {
+  CGLM_ALIGN_MAT mat4 rot;
+  glm_rotate_make(rot, angle, axis);
+  glm_mul_rot(rot, m, m);
+}
+
+/*!
+ * @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 transform
+ *
+ * @param[in, out]  m      affine transform
+ * @param[in]       pivot  rotation center
+ * @param[in]       angle  angle (radians)
+ * @param[in]       axis   axis
+ */
+CGLM_INLINE
+void
+glm_rotated_at(mat4 m, vec3 pivot, float angle, vec3 axis) {
+  CGLM_ALIGN(8) vec3 pivotInv;
+
+  glm_vec3_negate_to(pivot, pivotInv);
+
+  glm_translated(m, pivot);
+  glm_rotated(m, angle, axis);
+  glm_translated(m, pivotInv);
+}
+
+/*!
+ * @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 transform
+ *
+ * @param[in, out]  m      affine transform
+ * @param[in]       angle  angle (radians)
+ * @param[in]       axis   axis
+ */
+CGLM_INLINE
+void
+glm_spinned(mat4 m, float angle, vec3 axis) {
+  CGLM_ALIGN_MAT mat4 rot;
+  glm_rotate_atm(rot, m[3], angle, axis);
+  glm_mat4_mul(rot, m, m);
+}
+
+#endif /* cglm_affine_post_h */

include/cglm/affine-pre.h

@@ -0,0 +1,304 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglm_affine_pre_h
+#define cglm_affine_pre_h
+
+/*
+ Functions:
+   CGLM_INLINE void glm_translate_to(mat4 m, vec3 v, mat4 dest);
+   CGLM_INLINE void glm_translate(mat4 m, vec3 v);
+   CGLM_INLINE void glm_translate_x(mat4 m, float to);
+   CGLM_INLINE void glm_translate_y(mat4 m, float to);
+   CGLM_INLINE void glm_translate_z(mat4 m, float to);
+   CGLM_INLINE void glm_rotate_x(mat4 m, float angle, mat4 dest);
+   CGLM_INLINE void glm_rotate_y(mat4 m, float angle, mat4 dest);
+   CGLM_INLINE void glm_rotate_z(mat4 m, float angle, mat4 dest);
+   CGLM_INLINE void glm_rotate(mat4 m, float angle, vec3 axis);
+   CGLM_INLINE void glm_rotate_at(mat4 m, vec3 pivot, float angle, vec3 axis);
+   CGLM_INLINE void glm_rotate_atm(mat4 m, vec3 pivot, float angle, vec3 axis);
+   CGLM_INLINE void glm_spin(mat4 m, float angle, vec3 axis);
+ */
+
+#include "common.h"
+#include "util.h"
+#include "vec3.h"
+#include "vec4.h"
+#include "mat4.h"
+#include "affine-mat.h"
+
+/*!
+ * @brief translate existing transform matrix by v vector
+ *        and stores result in same matrix
+ *
+ * @param[in, out]  m  affine transform
+ * @param[in]       v  translate vector [x, y, z]
+ */
+CGLM_INLINE
+void
+glm_translate(mat4 m, vec3 v) {
+#if defined(CGLM_SIMD)
+  glmm_128 m0, m1, m2, m3;
+
+  m0 = glmm_load(m[0]);
+  m1 = glmm_load(m[1]);
+  m2 = glmm_load(m[2]);
+  m3 = glmm_load(m[3]);
+
+  glmm_store(m[3],
+             glmm_fmadd(m0, glmm_set1(v[0]),
+                        glmm_fmadd(m1, glmm_set1(v[1]),
+                                   glmm_fmadd(m2, glmm_set1(v[2]), m3))));
+#else
+  glm_vec4_muladds(m[0], v[0], m[3]);
+  glm_vec4_muladds(m[1], v[1], m[3]);
+  glm_vec4_muladds(m[2], v[2], m[3]);
+#endif
+}
+
+/*!
+ * @brief translate existing transform matrix by v vector
+ *        and store result in dest
+ *
+ * source matrix will remain same
+ *
+ * @param[in]  m    affine transform
+ * @param[in]  v    translate vector [x, y, z]
+ * @param[out] dest translated matrix
+ */
+CGLM_INLINE
+void
+glm_translate_to(mat4 m, vec3 v, mat4 dest) {
+  glm_mat4_copy(m, dest);
+  glm_translate(dest, v);
+}
+
+/*!
+ * @brief translate existing transform matrix by x factor
+ *
+ * @param[in, out]  m  affine transform
+ * @param[in]       x  x factor
+ */
+CGLM_INLINE
+void
+glm_translate_x(mat4 m, float x) {
+#if defined(CGLM_SIMD)
+  glmm_store(m[3], glmm_fmadd(glmm_load(m[0]), glmm_set1(x), glmm_load(m[3])));
+#else
+  vec4 v1;
+  glm_vec4_scale(m[0], x, v1);
+  glm_vec4_add(v1, m[3], m[3]);
+#endif
+}
+
+/*!
+ * @brief translate existing transform matrix by y factor
+ *
+ * @param[in, out]  m  affine transform
+ * @param[in]       y  y factor
+ */
+CGLM_INLINE
+void
+glm_translate_y(mat4 m, float y) {
+#if defined(CGLM_SIMD)
+  glmm_store(m[3], glmm_fmadd(glmm_load(m[1]), glmm_set1(y), glmm_load(m[3])));
+#else
+  vec4 v1;
+  glm_vec4_scale(m[1], y, v1);
+  glm_vec4_add(v1, m[3], m[3]);
+#endif
+}
+
+/*!
+ * @brief translate existing transform matrix by z factor
+ *
+ * @param[in, out]  m  affine transform
+ * @param[in]       z  z factor
+ */
+CGLM_INLINE
+void
+glm_translate_z(mat4 m, float z) {
+#if defined(CGLM_SIMD)
+  glmm_store(m[3], glmm_fmadd(glmm_load(m[2]), glmm_set1(z), glmm_load(m[3])));
+#else
+  vec4 v1;
+  glm_vec4_scale(m[2], z, v1);
+  glm_vec4_add(v1, m[3], m[3]);
+#endif
+}
+
+/*!
+ * @brief rotate existing transform matrix around X axis by angle
+ *        and store result in dest
+ *
+ * @param[in]   m      affine transform
+ * @param[in]   angle  angle (radians)
+ * @param[out]  dest   rotated matrix
+ */
+CGLM_INLINE
+void
+glm_rotate_x(mat4 m, float angle, mat4 dest) {
+  CGLM_ALIGN_MAT mat4 t = GLM_MAT4_IDENTITY_INIT;
+  float c, s;
+
+  c = cosf(angle);
+  s = sinf(angle);
+
+  t[1][1] =  c;
+  t[1][2] =  s;
+  t[2][1] = -s;
+  t[2][2] =  c;
+
+  glm_mul_rot(m, t, dest);
+}
+
+/*!
+ * @brief rotate existing transform matrix around Y axis by angle
+ *        and store result in dest
+ *
+ * @param[in]   m      affine transform
+ * @param[in]   angle  angle (radians)
+ * @param[out]  dest   rotated matrix
+ */
+CGLM_INLINE
+void
+glm_rotate_y(mat4 m, float angle, mat4 dest) {
+  CGLM_ALIGN_MAT mat4 t = GLM_MAT4_IDENTITY_INIT;
+  float c, s;
+
+  c = cosf(angle);
+  s = sinf(angle);
+
+  t[0][0] =  c;
+  t[0][2] = -s;
+  t[2][0] =  s;
+  t[2][2] =  c;
+
+  glm_mul_rot(m, t, dest);
+}
+
+/*!
+ * @brief rotate existing transform matrix around Z axis by angle
+ *        and store result in dest
+ *
+ * @param[in]   m      affine transform
+ * @param[in]   angle  angle (radians)
+ * @param[out]  dest   rotated matrix
+ */
+CGLM_INLINE
+void
+glm_rotate_z(mat4 m, float angle, mat4 dest) {
+  CGLM_ALIGN_MAT mat4 t = GLM_MAT4_IDENTITY_INIT;
+  float c, s;
+
+  c = cosf(angle);
+  s = sinf(angle);
+
+  t[0][0] =  c;
+  t[0][1] =  s;
+  t[1][0] = -s;
+  t[1][1] =  c;
+
+  glm_mul_rot(m, t, dest);
+}
+
+/*!
+ * @brief 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.
+ *
+ * @param[in, out]  m      affine transform
+ * @param[in]       angle  angle (radians)
+ * @param[in]       axis   axis
+ */
+CGLM_INLINE
+void
+glm_rotate(mat4 m, float angle, vec3 axis) {
+  CGLM_ALIGN_MAT mat4 rot;
+  glm_rotate_make(rot, angle, axis);
+  glm_mul_rot(m, rot, m);
+}
+
+/*!
+ * @brief rotate existing transform
+ *        around given axis by angle at given pivot point (rotation center)
+ *
+ * @param[in, out]  m      affine transform
+ * @param[in]       pivot  rotation center
+ * @param[in]       angle  angle (radians)
+ * @param[in]       axis   axis
+ */
+CGLM_INLINE
+void
+glm_rotate_at(mat4 m, vec3 pivot, float angle, vec3 axis) {
+  CGLM_ALIGN(8) vec3 pivotInv;
+
+  glm_vec3_negate_to(pivot, pivotInv);
+
+  glm_translate(m, pivot);
+  glm_rotate(m, angle, axis);
+  glm_translate(m, pivotInv);
+}
+
+/*!
+ * @brief 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.
+ *
+ * @param[out] m      affine transform
+ * @param[in]  pivot  rotation center
+ * @param[in]  angle  angle (radians)
+ * @param[in]  axis   axis
+ */
+CGLM_INLINE
+void
+glm_rotate_atm(mat4 m, vec3 pivot, float angle, vec3 axis) {
+  CGLM_ALIGN(8) vec3 pivotInv;
+
+  glm_vec3_negate_to(pivot, pivotInv);
+
+  glm_translate_make(m, pivot);
+  glm_rotate(m, angle, axis);
+  glm_translate(m, pivotInv);
+}
+
+/*!
+ * @brief rotate existing transform matrix 
+ *        around given axis by angle around self (doesn't affected by position)
+ *
+ * @param[in, out]  m      affine transform
+ * @param[in]       angle  angle (radians)
+ * @param[in]       axis   axis
+ */
+CGLM_INLINE
+void
+glm_spin(mat4 m, float angle, vec3 axis) {
+  CGLM_ALIGN_MAT mat4 rot;
+  glm_rotate_atm(rot, m[3], angle, axis);
+  glm_mat4_mul(m, rot, m);
+}
+
+#endif /* cglm_affine_pre_h */

include/cglm/affine.h

@@ -24,6 +24,7 @@
    CGLM_INLINE void glm_rotate(mat4 m, float angle, vec3 axis);
    CGLM_INLINE void glm_rotate_at(mat4 m, vec3 pivot, float angle, vec3 axis);
    CGLM_INLINE void glm_rotate_atm(mat4 m, vec3 pivot, float angle, vec3 axis);
+   CGLM_INLINE void glm_spin(mat4 m, float angle, vec3 axis);
    CGLM_INLINE void glm_decompose_scalev(mat4 m, vec3 s);
    CGLM_INLINE bool glm_uniscaled(mat4 m);
    CGLM_INLINE void glm_decompose_rs(mat4 m, mat4 r, vec3 s);
@@ -40,126 +41,10 @@
 #include "mat4.h"
 #include "affine-mat.h"
 
-/*!
- * @brief translate existing transform matrix by v vector
- *        and stores result in same matrix
- *
- * @param[in, out]  m  affine transfrom
- * @param[in]       v  translate vector [x, y, z]
- */
-CGLM_INLINE
-void
-glm_translate(mat4 m, vec3 v) {
-#if defined( __SSE__ ) || defined( __SSE2__ )
-  glmm_store(m[3],
-             _mm_add_ps(_mm_add_ps(_mm_mul_ps(glmm_load(m[0]),
-                                              _mm_set1_ps(v[0])),
-                                   _mm_mul_ps(glmm_load(m[1]),
-                                              _mm_set1_ps(v[1]))),
-                        _mm_add_ps(_mm_mul_ps(glmm_load(m[2]),
-                                              _mm_set1_ps(v[2])),
-                                   glmm_load(m[3]))))
-  ;
-#else
-  vec4 v1, v2, v3;
-
-  glm_vec4_scale(m[0], v[0], v1);
-  glm_vec4_scale(m[1], v[1], v2);
-  glm_vec4_scale(m[2], v[2], v3);
-
-  glm_vec4_add(v1, m[3], m[3]);
-  glm_vec4_add(v2, m[3], m[3]);
-  glm_vec4_add(v3, m[3], m[3]);
-#endif
-}
-
-/*!
- * @brief translate existing transform matrix by v vector
- *        and store result in dest
- *
- * source matrix will remain same
- *
- * @param[in]  m    affine transfrom
- * @param[in]  v    translate vector [x, y, z]
- * @param[out] dest translated matrix
- */
-CGLM_INLINE
-void
-glm_translate_to(mat4 m, vec3 v, mat4 dest) {
-  glm_mat4_copy(m, dest);
-  glm_translate(dest, v);
-}
-
-/*!
- * @brief translate existing transform matrix by x factor
- *
- * @param[in, out]  m  affine transfrom
- * @param[in]       x  x factor
- */
-CGLM_INLINE
-void
-glm_translate_x(mat4 m, float x) {
-#if defined( __SSE__ ) || defined( __SSE2__ )
-  glmm_store(m[3],
-             _mm_add_ps(_mm_mul_ps(glmm_load(m[0]),
-                                   _mm_set1_ps(x)),
-                        glmm_load(m[3])))
-  ;
-#else
-  vec4 v1;
-  glm_vec4_scale(m[0], x, v1);
-  glm_vec4_add(v1, m[3], m[3]);
-#endif
-}
-
-/*!
- * @brief translate existing transform matrix by y factor
- *
- * @param[in, out]  m  affine transfrom
- * @param[in]       y  y factor
- */
-CGLM_INLINE
-void
-glm_translate_y(mat4 m, float y) {
-#if defined( __SSE__ ) || defined( __SSE2__ )
-  glmm_store(m[3],
-             _mm_add_ps(_mm_mul_ps(glmm_load(m[1]),
-                                   _mm_set1_ps(y)),
-                        glmm_load(m[3])))
-  ;
-#else
-  vec4 v1;
-  glm_vec4_scale(m[1], y, v1);
-  glm_vec4_add(v1, m[3], m[3]);
-#endif
-}
-
-/*!
- * @brief translate existing transform matrix by z factor
- *
- * @param[in, out]  m  affine transfrom
- * @param[in]       z  z factor
- */
-CGLM_INLINE
-void
-glm_translate_z(mat4 m, float z) {
-#if defined( __SSE__ ) || defined( __SSE2__ )
-  glmm_store(m[3],
-             _mm_add_ps(_mm_mul_ps(glmm_load(m[2]),
-                                   _mm_set1_ps(z)),
-                        glmm_load(m[3])))
-  ;
-#else
-  vec4 v1;
-  glm_vec4_scale(m[2], z, v1);
-  glm_vec4_add(v1, m[3], m[3]);
-#endif
-}
-
 /*!
  * @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
@@ -173,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
  */
@@ -190,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
@@ -206,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
@@ -219,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
@@ -229,87 +114,12 @@ glm_scale_uni(mat4 m, float s) {
   glm_scale_to(m, v, m);
 }
 
-/*!
- * @brief rotate existing transform matrix around X axis by angle
- *        and store result in dest
- *
- * @param[in]   m      affine transfrom
- * @param[in]   angle  angle (radians)
- * @param[out]  dest   rotated matrix
- */
-CGLM_INLINE
-void
-glm_rotate_x(mat4 m, float angle, mat4 dest) {
-  CGLM_ALIGN_MAT mat4 t = GLM_MAT4_IDENTITY_INIT;
-  float c, s;
-
-  c = cosf(angle);
-  s = sinf(angle);
-
-  t[1][1] =  c;
-  t[1][2] =  s;
-  t[2][1] = -s;
-  t[2][2] =  c;
-
-  glm_mul_rot(m, t, dest);
-}
-
-/*!
- * @brief rotate existing transform matrix around Y axis by angle
- *        and store result in dest
- *
- * @param[in]   m      affine transfrom
- * @param[in]   angle  angle (radians)
- * @param[out]  dest   rotated matrix
- */
-CGLM_INLINE
-void
-glm_rotate_y(mat4 m, float angle, mat4 dest) {
-  CGLM_ALIGN_MAT mat4 t = GLM_MAT4_IDENTITY_INIT;
-  float c, s;
-
-  c = cosf(angle);
-  s = sinf(angle);
-
-  t[0][0] =  c;
-  t[0][2] = -s;
-  t[2][0] =  s;
-  t[2][2] =  c;
-
-  glm_mul_rot(m, t, dest);
-}
-
-/*!
- * @brief rotate existing transform matrix around Z axis by angle
- *        and store result in dest
- *
- * @param[in]   m      affine transfrom
- * @param[in]   angle  angle (radians)
- * @param[out]  dest   rotated matrix
- */
-CGLM_INLINE
-void
-glm_rotate_z(mat4 m, float angle, mat4 dest) {
-  CGLM_ALIGN_MAT mat4 t = GLM_MAT4_IDENTITY_INIT;
-  float c, s;
-
-  c = cosf(angle);
-  s = sinf(angle);
-
-  t[0][0] =  c;
-  t[0][1] =  s;
-  t[1][0] = -s;
-  t[1][1] =  c;
-
-  glm_mul_rot(m, t, dest);
-}
-
 /*!
  * @brief creates NEW rotation matrix by angle and axis
  *
  * 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
  */
@@ -337,67 +147,6 @@ glm_rotate_make(mat4 m, float angle, vec3 axis) {
   m[3][3] = 1.0f;
 }
 
-/*!
- * @brief rotate existing transform matrix around given axis by angle
- *
- * @param[in, out]  m      affine transfrom
- * @param[in]       angle  angle (radians)
- * @param[in]       axis   axis
- */
-CGLM_INLINE
-void
-glm_rotate(mat4 m, float angle, vec3 axis) {
-  CGLM_ALIGN_MAT mat4 rot;
-  glm_rotate_make(rot, angle, axis);
-  glm_mul_rot(m, rot, m);
-}
-
-/*!
- * @brief rotate existing transform
- *        around given axis by angle at given pivot point (rotation center)
- *
- * @param[in, out]  m      affine transfrom
- * @param[in]       pivot  rotation center
- * @param[in]       angle  angle (radians)
- * @param[in]       axis   axis
- */
-CGLM_INLINE
-void
-glm_rotate_at(mat4 m, vec3 pivot, float angle, vec3 axis) {
-  CGLM_ALIGN(8) vec3 pivotInv;
-
-  glm_vec3_negate_to(pivot, pivotInv);
-
-  glm_translate(m, pivot);
-  glm_rotate(m, angle, axis);
-  glm_translate(m, pivotInv);
-}
-
-/*!
- * @brief 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.
- *
- * @param[out] m      affine transfrom
- * @param[in]  pivot  rotation center
- * @param[in]  angle  angle (radians)
- * @param[in]  axis   axis
- */
-CGLM_INLINE
-void
-glm_rotate_atm(mat4 m, vec3 pivot, float angle, vec3 axis) {
-  CGLM_ALIGN(8) vec3 pivotInv;
-
-  glm_vec3_negate_to(pivot, pivotInv);
-
-  glm_translate_make(m, pivot);
-  glm_rotate(m, angle, axis);
-  glm_translate(m, pivotInv);
-}
-
 /*!
  * @brief decompose scale vector
  *
@@ -455,7 +204,7 @@ glm_decompose_rs(mat4 m, mat4 r, vec3 s) {
   glm_vec4_scale(r[1], 1.0f/s[1], r[1]);
   glm_vec4_scale(r[2], 1.0f/s[2], r[2]);
 
-  /* Note from Apple Open Source (asume that the matrix is orthonormal):
+  /* Note from Apple Open Source (assume that the matrix is orthonormal):
      check for a coordinate system flip.  If the determinant
      is -1, then negate the matrix and the scaling factors. */
   glm_vec3_cross(m[0], m[1], v);
@@ -471,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]
@@ -483,4 +232,7 @@ glm_decompose(mat4 m, vec4 t, mat4 r, vec3 s) {
   glm_decompose_rs(m, r, s);
 }
 
+#include "affine-pre.h"
+#include "affine-post.h"
+
 #endif /* cglm_affine_h */

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
@@ -228,6 +228,8 @@ glm_aabb_aabb(vec3 box[2], vec3 other[2]) {
  * https://github.com/erich666/GraphicsGems/blob/master/gems/BoxSphere.c
  * Solid Box - Solid Sphere test.
  *
+ * Sphere Representation in cglm: [center.x, center.y, center.z, radii]
+ *
  * @param[in]   box    solid bounding box
  * @param[in]   s      solid sphere
  */
@@ -237,13 +239,13 @@ glm_aabb_sphere(vec3 box[2], vec4 s) {
   float dmin;
   int   a, b, c;
 
-  a = s[0] >= box[0][0];
+  a = (s[0] < box[0][0]) + (s[0] > box[1][0]);
-  b = s[1] >= box[0][1];
+  b = (s[1] < box[0][1]) + (s[1] > box[1][1]);
-  c = s[2] >= box[0][2];
+  c = (s[2] < box[0][2]) + (s[2] > box[1][2]);
 
-  dmin  = glm_pow2(s[0] - box[a][0])
+  dmin  = glm_pow2((s[0] - box[!(a - 1)][0]) * (a != 0))
-        + glm_pow2(s[1] - box[b][1])
+        + glm_pow2((s[1] - box[!(b - 1)][1]) * (b != 0))
-        + glm_pow2(s[2] - box[c][2]);
+        + glm_pow2((s[2] - box[!(c - 1)][2]) * (c != 0));
 
   return dmin <= glm_pow2(s[3]);
 }

include/cglm/call.h

@@ -15,15 +15,25 @@ extern "C" {
 #include "call/vec2.h"
 #include "call/vec3.h"
 #include "call/vec4.h"
+#include "call/ivec2.h"
+#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/affine.h

@@ -81,6 +81,10 @@ CGLM_EXPORT
 void
 glmc_rotate_atm(mat4 m, vec3 pivot, float angle, vec3 axis);
 
+CGLM_EXPORT
+void
+glmc_spin(mat4 m, float angle, vec3 axis);
+
 CGLM_EXPORT
 void
 glmc_decompose_scalev(mat4 m, vec3 s);
@@ -97,6 +101,52 @@ CGLM_EXPORT
 void
 glmc_decompose(mat4 m, vec4 t, mat4 r, vec3 s);
 
+/* affine-post */
+
+CGLM_EXPORT
+void
+glmc_translated(mat4 m, vec3 v);
+
+CGLM_EXPORT
+void
+glmc_translated_to(mat4 m, vec3 v, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_translated_x(mat4 m, float x);
+
+CGLM_EXPORT
+void
+glmc_translated_y(mat4 m, float y);
+
+CGLM_EXPORT
+void
+glmc_translated_z(mat4 m, float z);
+
+CGLM_EXPORT
+void
+glmc_rotated_x(mat4 m, float angle, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_rotated_y(mat4 m, float angle, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_rotated_z(mat4 m, float angle, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_rotated(mat4 m, float angle, vec3 axis);
+
+CGLM_EXPORT
+void
+glmc_rotated_at(mat4 m, vec3 pivot, float angle, vec3 axis);
+
+CGLM_EXPORT
+void
+glmc_spinned(mat4 m, float angle, vec3 axis);
+
 /* affine-mat */
 
 CGLM_EXPORT

include/cglm/call/cam.h

@@ -15,22 +15,16 @@ extern "C" {
 
 CGLM_EXPORT
 void
-glmc_frustum(float left,
+glmc_frustum(float left,   float right,
-             float right,
+             float bottom, float top,
-             float bottom,
+             float nearZ,  float farZ,
-             float top,
-             float nearVal,
-             float farVal,
              mat4  dest);
 
 CGLM_EXPORT
 void
-glmc_ortho(float left,
+glmc_ortho(float left,   float right,
-           float right,
+           float bottom, float top,
-           float bottom,
+           float nearZ,  float farZ,
-           float top,
-           float nearVal,
-           float farVal,
            mat4  dest);
 
 CGLM_EXPORT
@@ -55,11 +49,7 @@ glmc_ortho_default_s(float aspect, float size, mat4 dest);
 
 CGLM_EXPORT
 void
-glmc_perspective(float fovy,
+glmc_perspective(float fovy, float aspect, float nearZ, float farZ, mat4 dest);
-                 float aspect,
-                 float nearVal,
-                 float farVal,
-                 mat4 dest);
 
 CGLM_EXPORT
 void
@@ -88,8 +78,8 @@ glmc_look_anyup(vec3 eye, vec3 dir, mat4 dest);
 CGLM_EXPORT
 void
 glmc_persp_decomp(mat4 proj,
-                  float * __restrict nearVal,
+                  float * __restrict nearZ,
-                  float * __restrict farVal,
+                  float * __restrict farZ,
                   float * __restrict top,
                   float * __restrict bottom,
                   float * __restrict left,
@@ -114,16 +104,16 @@ glmc_persp_decomp_y(mat4 proj,
 CGLM_EXPORT
 void
 glmc_persp_decomp_z(mat4 proj,
-                    float * __restrict nearVal,
+                    float * __restrict nearZ,
-                    float * __restrict farVal);
+                    float * __restrict farZ);
 
 CGLM_EXPORT
 void
-glmc_persp_decomp_far(mat4 proj, float * __restrict farVal);
+glmc_persp_decomp_far(mat4 proj, float * __restrict farZ);
 
 CGLM_EXPORT
 void
-glmc_persp_decomp_near(mat4 proj, float * __restrict nearVal);
+glmc_persp_decomp_near(mat4 proj, float * __restrict nearZ);
 
 CGLM_EXPORT
 float

include/cglm/call/clipspace/ortho_lh_no.h

@@ -0,0 +1,46 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_ortho_lh_no_h
+#define cglmc_ortho_lh_no_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_ortho_lh_no(float left,    float right,
+                 float bottom,  float top,
+                 float nearZ,   float farZ,
+                 mat4  dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_aabb_lh_no(vec3 box[2], mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_aabb_p_lh_no(vec3 box[2], float padding, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_aabb_pz_lh_no(vec3 box[2], float padding, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_default_lh_no(float aspect, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_default_s_lh_no(float aspect, float size, mat4 dest);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_ortho_lh_no_h */

include/cglm/call/clipspace/ortho_lh_zo.h

@@ -0,0 +1,46 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_ortho_lh_zo_h
+#define cglmc_ortho_lh_zo_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_ortho_lh_zo(float left,    float right,
+                 float bottom,  float top,
+                 float nearZ,   float farZ,
+                 mat4  dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_aabb_lh_zo(vec3 box[2], mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_aabb_p_lh_zo(vec3 box[2], float padding, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_aabb_pz_lh_zo(vec3 box[2], float padding, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_default_lh_zo(float aspect, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_default_s_lh_zo(float aspect, float size, mat4 dest);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_ortho_lh_zo_h */

include/cglm/call/clipspace/ortho_rh_no.h

@@ -0,0 +1,46 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_ortho_rh_no_h
+#define cglmc_ortho_rh_no_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_ortho_rh_no(float left,    float right,
+                 float bottom,  float top,
+                 float nearZ,   float farZ,
+                 mat4  dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_aabb_rh_no(vec3 box[2], mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_aabb_p_rh_no(vec3 box[2], float padding, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_aabb_pz_rh_no(vec3 box[2], float padding, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_default_rh_no(float aspect, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_default_s_rh_no(float aspect, float size, mat4 dest);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_ortho_rh_no_h */

include/cglm/call/clipspace/ortho_rh_zo.h

@@ -0,0 +1,46 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_ortho_rh_zo_h
+#define cglmc_ortho_rh_zo_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_ortho_rh_zo(float left,    float right,
+                 float bottom,  float top,
+                 float nearZ,   float farZ,
+                 mat4  dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_aabb_rh_zo(vec3 box[2], mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_aabb_p_rh_zo(vec3 box[2], float padding, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_aabb_pz_rh_zo(vec3 box[2], float padding, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_default_rh_zo(float aspect, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_ortho_default_s_rh_zo(float aspect, float size, mat4 dest);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_ortho_rh_zo_h */

include/cglm/call/clipspace/persp_lh_no.h

@@ -0,0 +1,87 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_persp_lh_no_h
+#define cglmc_persp_lh_no_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_frustum_lh_no(float left,    float right,
+                   float bottom,  float top,
+                   float nearZ,   float farZ,
+                   mat4  dest);
+
+CGLM_EXPORT
+void
+glmc_perspective_lh_no(float fovy,
+                       float aspect,
+                       float nearVal,
+                       float farVal,
+                       mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_persp_move_far_lh_no(mat4 proj, float deltaFar);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_lh_no(mat4 proj,
+                        float * __restrict nearZ, float * __restrict farZ,
+                        float * __restrict top,   float * __restrict bottom,
+                        float * __restrict left,  float * __restrict right);
+
+CGLM_EXPORT
+void
+glmc_persp_decompv_lh_no(mat4 proj, float dest[6]);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_x_lh_no(mat4 proj,
+                          float * __restrict left,
+                          float * __restrict right);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_y_lh_no(mat4 proj,
+                          float * __restrict top,
+                          float * __restrict bottom);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_z_lh_no(mat4 proj,
+                          float * __restrict nearZ,
+                          float * __restrict farZ);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_far_lh_no(mat4 proj, float * __restrict farZ);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_near_lh_no(mat4 proj, float * __restrict nearZ);
+
+CGLM_EXPORT
+void
+glmc_persp_sizes_lh_no(mat4 proj, float fovy, vec4 dest);
+
+CGLM_EXPORT
+float
+glmc_persp_fovy_lh_no(mat4 proj);
+
+CGLM_EXPORT
+float
+glmc_persp_aspect_lh_no(mat4 proj);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_persp_lh_no_h */

include/cglm/call/clipspace/persp_lh_zo.h

@@ -0,0 +1,87 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_persp_lh_zo_h
+#define cglmc_persp_lh_zo_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_frustum_lh_zo(float left,    float right,
+                   float bottom,  float top,
+                   float nearZ,   float farZ,
+                   mat4  dest);
+
+CGLM_EXPORT
+void
+glmc_perspective_lh_zo(float fovy,
+                       float aspect,
+                       float nearVal,
+                       float farVal,
+                       mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_persp_move_far_lh_zo(mat4 proj, float deltaFar);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_lh_zo(mat4 proj,
+                        float * __restrict nearZ, float * __restrict farZ,
+                        float * __restrict top,   float * __restrict bottom,
+                        float * __restrict left,  float * __restrict right);
+
+CGLM_EXPORT
+void
+glmc_persp_decompv_lh_zo(mat4 proj, float dest[6]);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_x_lh_zo(mat4 proj,
+                          float * __restrict left,
+                          float * __restrict right);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_y_lh_zo(mat4 proj,
+                          float * __restrict top,
+                          float * __restrict bottom);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_z_lh_zo(mat4 proj,
+                          float * __restrict nearZ,
+                          float * __restrict farZ);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_far_lh_zo(mat4 proj, float * __restrict farZ);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_near_lh_zo(mat4 proj, float * __restrict nearZ);
+
+CGLM_EXPORT
+void
+glmc_persp_sizes_lh_zo(mat4 proj, float fovy, vec4 dest);
+
+CGLM_EXPORT
+float
+glmc_persp_fovy_lh_zo(mat4 proj);
+
+CGLM_EXPORT
+float
+glmc_persp_aspect_lh_zo(mat4 proj);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_persp_lh_zo_h */

include/cglm/call/clipspace/persp_rh_no.h

@@ -0,0 +1,87 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_persp_rh_no_h
+#define cglmc_persp_rh_no_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_frustum_rh_no(float left,    float right,
+                   float bottom,  float top,
+                   float nearZ,   float farZ,
+                   mat4  dest);
+
+CGLM_EXPORT
+void
+glmc_perspective_rh_no(float fovy,
+                       float aspect,
+                       float nearVal,
+                       float farVal,
+                       mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_persp_move_far_rh_no(mat4 proj, float deltaFar);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_rh_no(mat4 proj,
+                        float * __restrict nearZ, float * __restrict farZ,
+                        float * __restrict top,   float * __restrict bottom,
+                        float * __restrict left,  float * __restrict right);
+
+CGLM_EXPORT
+void
+glmc_persp_decompv_rh_no(mat4 proj, float dest[6]);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_x_rh_no(mat4 proj,
+                          float * __restrict left,
+                          float * __restrict right);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_y_rh_no(mat4 proj,
+                          float * __restrict top,
+                          float * __restrict bottom);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_z_rh_no(mat4 proj,
+                          float * __restrict nearZ,
+                          float * __restrict farZ);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_far_rh_no(mat4 proj, float * __restrict farZ);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_near_rh_no(mat4 proj, float * __restrict nearZ);
+
+CGLM_EXPORT
+void
+glmc_persp_sizes_rh_no(mat4 proj, float fovy, vec4 dest);
+
+CGLM_EXPORT
+float
+glmc_persp_fovy_rh_no(mat4 proj);
+
+CGLM_EXPORT
+float
+glmc_persp_aspect_rh_no(mat4 proj);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_persp_rh_no_h */

include/cglm/call/clipspace/persp_rh_zo.h

@@ -0,0 +1,87 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_persp_rh_zo_h
+#define cglmc_persp_rh_zo_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_frustum_rh_zo(float left,    float right,
+                   float bottom,  float top,
+                   float nearZ,   float farZ,
+                   mat4  dest);
+
+CGLM_EXPORT
+void
+glmc_perspective_rh_zo(float fovy,
+                       float aspect,
+                       float nearVal,
+                       float farVal,
+                       mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_persp_move_far_rh_zo(mat4 proj, float deltaFar);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_rh_zo(mat4 proj,
+                        float * __restrict nearZ, float * __restrict farZ,
+                        float * __restrict top,   float * __restrict bottom,
+                        float * __restrict left,  float * __restrict right);
+
+CGLM_EXPORT
+void
+glmc_persp_decompv_rh_zo(mat4 proj, float dest[6]);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_x_rh_zo(mat4 proj,
+                          float * __restrict left,
+                          float * __restrict right);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_y_rh_zo(mat4 proj,
+                          float * __restrict top,
+                          float * __restrict bottom);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_z_rh_zo(mat4 proj,
+                          float * __restrict nearZ,
+                          float * __restrict farZ);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_far_rh_zo(mat4 proj, float * __restrict farZ);
+
+CGLM_EXPORT
+void
+glmc_persp_decomp_near_rh_zo(mat4 proj, float * __restrict nearZ);
+
+CGLM_EXPORT
+void
+glmc_persp_sizes_rh_zo(mat4 proj, float fovy, vec4 dest);
+
+CGLM_EXPORT
+float
+glmc_persp_fovy_rh_zo(mat4 proj);
+
+CGLM_EXPORT
+float
+glmc_persp_aspect_rh_zo(mat4 proj);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_persp_rh_zo_h */

include/cglm/call/clipspace/project_no.h

@@ -0,0 +1,31 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_project_no_h
+#define cglmc_project_no_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_unprojecti_no(vec3 pos, mat4 invMat, vec4 vp, vec3 dest);
+
+CGLM_EXPORT
+void
+glmc_project_no(vec3 pos, mat4 m, vec4 vp, vec3 dest);
+
+CGLM_EXPORT
+float
+glmc_project_z_no(vec3 pos, mat4 m);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_project_no_h */

include/cglm/call/clipspace/project_zo.h

@@ -0,0 +1,31 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_project_zo_h
+#define cglmc_project_zo_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_unprojecti_zo(vec3 pos, mat4 invMat, vec4 vp, vec3 dest);
+
+CGLM_EXPORT
+void
+glmc_project_zo(vec3 pos, mat4 m, vec4 vp, vec3 dest);
+
+CGLM_EXPORT
+float
+glmc_project_z_zo(vec3 pos, mat4 m);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_project_zo_h */

include/cglm/call/clipspace/view_lh_no.h

@@ -0,0 +1,31 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_view_lh_no_h
+#define cglmc_view_lh_no_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_lookat_lh_no(vec3 eye, vec3 center, vec3 up, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_look_lh_no(vec3 eye, vec3 dir, vec3 up, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_look_anyup_lh_no(vec3 eye, vec3 dir, mat4 dest);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_view_lh_no_h */

include/cglm/call/clipspace/view_lh_zo.h

@@ -0,0 +1,31 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_view_lh_zo_h
+#define cglmc_view_lh_zo_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_lookat_lh_zo(vec3 eye, vec3 center, vec3 up, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_look_lh_zo(vec3 eye, vec3 dir, vec3 up, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_look_anyup_lh_zo(vec3 eye, vec3 dir, mat4 dest);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_view_lh_zo_h */

include/cglm/call/clipspace/view_rh_no.h

@@ -0,0 +1,31 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_view_rh_no_h
+#define cglmc_view_rh_no_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_lookat_rh_no(vec3 eye, vec3 center, vec3 up, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_look_rh_no(vec3 eye, vec3 dir, vec3 up, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_look_anyup_rh_no(vec3 eye, vec3 dir, mat4 dest);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_view_rh_no_h */

include/cglm/call/clipspace/view_rh_zo.h

@@ -0,0 +1,31 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_view_rh_zo_h
+#define cglmc_view_rh_zo_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_lookat_rh_zo(vec3 eye, vec3 center, vec3 up, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_look_rh_zo(vec3 eye, vec3 dir, vec3 up, mat4 dest);
+
+CGLM_EXPORT
+void
+glmc_look_anyup_rh_zo(vec3 eye, vec3 dir, mat4 dest);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_view_rh_zo_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

@@ -0,0 +1,179 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_ivec2_h
+#define cglmc_ivec2_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_ivec2(int * __restrict v, ivec2 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec2_copy(ivec2 a, ivec2 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec2_zero(ivec2 v);
+
+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);
+
+CGLM_EXPORT
+void
+glmc_ivec2_adds(ivec2 v, int s, ivec2 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec2_sub(ivec2 a, ivec2 b, ivec2 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec2_subs(ivec2 v, int s, ivec2 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec2_mul(ivec2 a, ivec2 b, ivec2 dest);
+
+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);
+
+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);
+
+CGLM_EXPORT
+void
+glmc_ivec2_minv(ivec2 a, ivec2 b, ivec2 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec2_clamp(ivec2 v, int minVal, int maxVal);
+
+CGLM_EXPORT
+void
+glmc_ivec2_abs(ivec2 v, ivec2 dest);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_ivec2_h */

include/cglm/call/ivec3.h

@@ -0,0 +1,183 @@
+/*
+ * Copyright (c);, Recep Aslantas.
+ *
+ * MIT License (MIT);, http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_ivec3_h
+#define cglmc_ivec3_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_ivec3(ivec4 v4, ivec3 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec3_copy(ivec3 a, ivec3 dest);
+
+CGLM_EXPORT
+void 
+glmc_ivec3_zero(ivec3 v);
+
+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);
+
+CGLM_EXPORT
+void
+glmc_ivec3_adds(ivec3 v, int s, ivec3 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec3_sub(ivec3 a, ivec3 b, ivec3 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec3_subs(ivec3 v, int s, ivec3 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec3_mul(ivec3 a, ivec3 b, ivec3 dest);
+
+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);
+
+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);
+
+CGLM_EXPORT
+void
+glmc_ivec3_minv(ivec3 a, ivec3 b, ivec3 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec3_clamp(ivec3 v, int minVal, int maxVal);
+
+CGLM_EXPORT
+void
+glmc_ivec3_abs(ivec3 v, ivec3 dest);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_ivec3_h */

include/cglm/call/ivec4.h

@@ -0,0 +1,147 @@
+/*
+ * Copyright (c), Recep Aslantas.
+ *
+ * MIT License (MIT), http://opensource.org/licenses/MIT
+ * Full license can be found in the LICENSE file
+ */
+
+#ifndef cglmc_ivec4_h
+#define cglmc_ivec4_h
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#include "../cglm.h"
+
+CGLM_EXPORT
+void
+glmc_ivec4(ivec3 v3, int last, ivec4 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec4_copy(ivec4 a, ivec4 dest);
+
+CGLM_EXPORT
+void 
+glmc_ivec4_zero(ivec4 v);
+
+CGLM_EXPORT
+void
+glmc_ivec4_one(ivec4 v);
+
+CGLM_EXPORT
+void
+glmc_ivec4_add(ivec4 a, ivec4 b, ivec4 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec4_adds(ivec4 v, int s, ivec4 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec4_sub(ivec4 a, ivec4 b, ivec4 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec4_subs(ivec4 v, int s, ivec4 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec4_mul(ivec4 a, ivec4 b, ivec4 dest);
+
+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);
+
+CGLM_EXPORT
+float
+glmc_ivec4_distance(ivec4 a, ivec4 b);
+
+CGLM_EXPORT
+void
+glmc_ivec4_maxv(ivec4 a, ivec4 b, ivec4 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec4_minv(ivec4 a, ivec4 b, ivec4 dest);
+
+CGLM_EXPORT
+void
+glmc_ivec4_clamp(ivec4 v, int minVal, int maxVal);
+
+CGLM_EXPORT
+void
+glmc_ivec4_abs(ivec4 v, ivec4 dest);
+
+#ifdef __cplusplus
+}
+#endif
+#endif /* cglmc_ivec4_h */

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(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(float * __restrict src, mat2x3 dest);
+
+CGLM_EXPORT
+void
+glmc_mat2x3_mul(mat2x3 m1, mat3x2 m2, mat2 dest);
+
+CGLM_EXPORT
+void
+glmc_mat2x3_mulv(mat2x3 m, vec3 v, vec2 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(float * __restrict src, mat2x4 dest);
+
+CGLM_EXPORT
+void
+glmc_mat2x4_mul(mat2x4 m1, mat4x2 m2, mat2 dest);
+
+CGLM_EXPORT
+void
+glmc_mat2x4_mulv(mat2x4 m, vec4 v, vec2 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(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(float * __restrict src, mat3x2 dest);
+
+CGLM_EXPORT
+void
+glmc_mat3x2_mul(mat3x2 m1, mat2x3 m2, mat3 dest);
+
+CGLM_EXPORT
+void
+glmc_mat3x2_mulv(mat3x2 m, vec2 v, vec3 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 */