Skip to content

feat: add comprehensive linear algebra operations with BLAS backend #13

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
CodeWithKyrian merged 1 commit into from
Apr 17, 2026
Merged

feat: add comprehensive linear algebra operations with BLAS backend #13

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
49 changes: 9 additions & 40 deletions .github/workflows/release.yml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -43,48 +43,17 @@ jobs:
shell: bash
run: ./scripts/build.sh release

- name: Create distribution package (Unix)
if: matrix.archive_ext == 'tar.gz'
shell: bash
run: |
DIST_NAME="dist-${{ matrix.platform }}"
mkdir -p "${DIST_NAME}"

# Copy all files except excluded directories
find . -maxdepth 1 \
-not -path "./.git" \
-not -path "." \
-not -path "./${DIST_NAME}" \
-not -path "./rust" \
-not -path "./docs" \
-exec cp -r {} "${DIST_NAME}/" \;

# Create tarball
tar -czf "${DIST_NAME}.tar.gz" "${DIST_NAME}"

echo "Created ${DIST_NAME}.tar.gz"
ls -lh "${DIST_NAME}.tar.gz"

- name: Create distribution package (Windows)
if: matrix.archive_ext == 'zip'
shell: pwsh
- name: Install zip (Windows)
if: runner.os == 'Windows'
shell: powershell
run: |
$distName = "dist-${{ matrix.platform }}"
New-Item -ItemType Directory -Name $distName -Force

# Copy all files except excluded directories
$exclude = @('.git', $distName, 'rust', 'docs')
Get-ChildItem -Path . -Force | Where-Object {
$_.Name -notin $exclude
} | ForEach-Object {
Copy-Item -Path $_.FullName -Destination $distName -Recurse -Force
if (!(Get-Command zip -ErrorAction SilentlyContinue)) {
choco install zip --no-progress -y
}

# Create zip
Compress-Archive -Path $distName -DestinationPath "${distName}.zip" -Force

Write-Host "Created ${distName}.zip"
Get-ChildItem "${distName}.zip"

- name: Create distribution package
shell: bash
run: ./scripts/package.sh '${{ matrix.platform }}'

- name: Upload to GitHub Release
if: github.event_name == 'release'
Expand Down
36 changes: 18 additions & 18 deletions SPEC.md
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -334,29 +334,29 @@ Sort kind selection is enum-based via `SortKind`:
**Priority**: MEDIUM (requires BLAS)

**Requirements**:
- [ ] 9.2.1: `$array->svd($full_matrices = true)` - Singular value decomposition
- [ ] 9.2.2: `$array->qr()` - QR decomposition
- [ ] 9.2.3: `$array->cholesky()` - Cholesky decomposition
- [ ] 9.2.4: `$array->eig()` - Eigenvalue decomposition
- [ ] 9.2.5: `$array->lu()` - LU decomposition
- [x] 9.2.1: `$array->svd($full_matrices = true)` - Singular value decomposition
- [x] 9.2.2: `$array->qr()` - QR decomposition
- [x] 9.2.3: `$array->cholesky()` - Cholesky decomposition
- [ ] 9.2.4: `$array->eig()` - Eigenvalue decomposition (blocked on complex number support)
- [ ] 9.2.5: `$array->lu()` - LU decomposition (not exposed; used internally by solve/inv/det)

### 9.3 Matrix Properties (REQ-9.3)
**Priority**: MEDIUM (requires BLAS)

**Requirements**:
- [ ] 9.3.1: `$array->det()` - Determinant
- [ ] 9.3.2: `$array->inv()` - Matrix inverse
- [ ] 9.3.3: `$array->pinv()` - Pseudo-inverse
- [x] 9.3.1: `$array->det()` - Determinant
- [x] 9.3.2: `$array->inv()` - Matrix inverse
- [x] 9.3.3: `$array->pinv()` - Pseudo-inverse
- [x] 9.3.4: `$array->norm($ord = null, $axis = null, $keepdims = false)` - Matrix/vector norm (supports 1, 2, inf, -inf, fro)
- [ ] 9.3.5: `$array->cond($p = null)` - Condition number
- [ ] 9.3.6: `$array->rank()` - Matrix rank
- [x] 9.3.5: `$array->cond($p = null)` - Condition number
- [x] 9.3.6: `$array->rank()` - Matrix rank

### 9.4 Solving Linear Systems (REQ-9.4)
**Priority**: MEDIUM (requires BLAS)

**Requirements**:
- [ ] 9.4.1: `NDArray::solve($a, $b)` - Solve linear equations
- [ ] 9.4.2: `NDArray::lstsq($a, $b)` - Least squares solution
- [x] 9.4.1: `$a->solve($b)` - Solve linear equations
- [x] 9.4.2: `$a->lstsq($b)` - Least squares solution

## 10. Iteration and Application

Expand Down Expand Up @@ -419,11 +419,11 @@ Sort kind selection is enum-based via `SortKind`:
**Priority**: MEDIUM

**Requirements**:
- [ ] 12.2.1: Matrix multiplication (GEMM)
- [ ] 12.2.2: Dot products (DOT)
- [ ] 12.2.3: Matrix-vector operations (GEMV)
- [ ] 12.2.4: Linear system solving
- [ ] 12.2.5: Eigenvalue/SVD decompositions
- [x] 12.2.1: Matrix multiplication (GEMM)
- [x] 12.2.2: Dot products (DOT)
- [x] 12.2.3: Matrix-vector operations (GEMV)
- [x] 12.2.4: Linear system solving
- [x] 12.2.5: SVD decompositions (eigenvalue blocked on complex number support)

## 13. Memory Management

Expand Down Expand Up @@ -617,7 +617,7 @@ Sort kind selection is enum-based via `SortKind`:
| Broadcasting | Automatic | CRITICAL | Done |
| Matrix mult | `@` operator | HIGH | Planned |
| Reductions | `sum()`, `mean()` | HIGH | Planned |
| SVD | `np.linalg.svd()` | MEDIUM | Planned |
| SVD | `np.linalg.svd()` | MEDIUM | Done |
| FFT | `np.fft.fft()` | LOW | Future |

## Appendix B: Glossary
Expand Down
Loading