numerical

May 10, 2023

Changes in this version

+ type FiniteVector interface

+ At func(i int) float64

+ Len func() int

+ WithDim func(i int, value float64) T

+ type GridSearch3D struct

+ Recursions int

+ XStops int

+ YStops int

+ ZStops int

+ func (g *GridSearch3D) Maximize(min, max Vec3, f func(Vec3) float64) (Vec3, float64)

+ func (g *GridSearch3D) Minimize(min, max Vec3, f func(Vec3) float64) (Vec3, float64)

+ type RecursiveLineSearch struct

+ LineSearch LineSearch

+ func (r *RecursiveLineSearch[T]) Maximize(min, max T, f func(T) float64) (solution T, fVal float64)

+ func (r *RecursiveLineSearch[T]) Minimize(min, max T, f func(T) float64) (solution T, fVal float64)

type Vec

+ func (v Vec) At(i int) float64

+ func (v Vec) Len() int

+ func (v Vec) WithDim(i int, x float64) Vec

type Vec2

+ func (v Vec2) At(i int) float64

+ func (v Vec2) Len() int

+ func (v Vec2) WithDim(i int, x float64) Vec2

type Vec3

+ func (v Vec3) At(i int) float64

+ func (v Vec3) Len() int

+ func (v Vec3) WithDim(i int, x float64) Vec3

type Vec4

+ func (v Vec4) At(i int) float64

+ func (v Vec4) Len() int

+ func (v Vec4) WithDim(i int, x float64) Vec4

Aug 3, 2022

Changes in this version

+ type BiCGSTAB struct

+ B Vec

+ Op func(v Vec) Vec

+ func NewBiCGSTAB(op func(v Vec) Vec, b, initGuess Vec) *BiCGSTAB

+ func (b *BiCGSTAB) Iter() Vec

+ type BiCGSTABSolver struct

+ MAETolerance float64

+ MSETolerance float64

+ MaxIters int

+ func (b *BiCGSTABSolver) SolveLinearSystem(op func(v Vec) Vec, x, initGuess Vec) Vec

+ type GridSearch2D struct

+ Recursions int

+ XStops int

+ YStops int

+ func (g *GridSearch2D) Maximize(min, max Vec2, f func(Vec2) float64) (Vec2, float64)

+ func (g *GridSearch2D) Minimize(min, max Vec2, f func(Vec2) float64) (Vec2, float64)

+ type LargeLinearSolver interface

+ SolveLinearSystem func(op func(v Vec) Vec, b, initGuess Vec) Vec

+ type LineSearch struct

+ Recursions int

+ Stops int

+ func (l *LineSearch) Maximize(min, max float64, f func(float64) float64) (x, fVal float64)

+ func (l *LineSearch) Minimize(min, max float64, f func(float64) float64) (x, fVal float64)

type SparseCholesky

+ func (s *SparseCholesky) ApplyInverseVec2(x []Vec2) []Vec2

+ func (s *SparseCholesky) ApplyVec2(x []Vec2) []Vec2

type SparseMatrix

+ func (s *SparseMatrix) Apply(x Vec) Vec

+ func (s *SparseMatrix) ApplyVec2(x []Vec2) []Vec2

+ func (s *SparseMatrix) Transpose() *SparseMatrix

type Vec

+ func (v Vec) Zeros() Vec

type Vec2

+ func (v Vec2) DistSquared(v1 Vec2) float64

+ func (v Vec2) Max(v1 Vec2) Vec2

+ func (v Vec2) Min(v1 Vec2) Vec2

+ func (v Vec2) Zeros() Vec2

type Vec3

+ func (v Vec3) DistSquared(v1 Vec3) float64

+ func (v Vec3) Max(v1 Vec3) Vec3

+ func (v Vec3) Min(v1 Vec3) Vec3

+ func (v Vec3) Zeros() Vec3

type Vec4

+ func (v Vec4) DistSquared(v1 Vec4) float64

+ func (v Vec4) Max(v1 Vec4) Vec4

+ func (v Vec4) Min(v1 Vec4) Vec4

+ func (v Vec4) Zeros() Vec4

+ type Vector interface

+ Add func(T) T

+ Dist func(T) float64

+ DistSquared func(T) float64

+ Max func(T) T

+ Min func(T) T

+ Norm func() float64

+ Scale func(float64) T

+ Sub func(T) T

+ Zeros func() T

Jul 9, 2022

Changes in this version

+ type Matrix2 [4]float64

+ func NewMatrix2Columns(c1, c2 Vec2) *Matrix2

+ func NewMatrix2Rotation(theta float64) *Matrix2

+ func (m *Matrix2) Add(m1 *Matrix2) *Matrix2

+ func (m *Matrix2) Det() float64

+ func (m *Matrix2) Eigenvalues() [2]complex128

+ func (m *Matrix2) Inverse() *Matrix2

+ func (m *Matrix2) InvertInPlace()

+ func (m *Matrix2) InvertInPlaceDet(det float64)

+ func (m *Matrix2) Mul(m1 *Matrix2) *Matrix2

+ func (m *Matrix2) MulColumn(c Vec2) Vec2

+ func (m *Matrix2) MulColumnInv(c Vec2, det float64) Vec2

+ func (m *Matrix2) SVD(u, s, v *Matrix2)

+ func (m *Matrix2) Scale(s float64)

+ func (m *Matrix2) Transpose() *Matrix2

+ type Matrix3 [9]float64

+ func NewMatrix3Columns(c1, c2, c3 Vec3) *Matrix3

+ func NewMatrix3Rotation(axis Vec3, angle float64) *Matrix3

+ func (m *Matrix3) Add(m1 *Matrix3) *Matrix3

+ func (m *Matrix3) Cols() [3]Vec3

+ func (m *Matrix3) Det() float64

+ func (m *Matrix3) Eigenvalues() [3]complex128

+ func (m *Matrix3) Inverse() *Matrix3

+ func (m *Matrix3) InvertInPlace()

+ func (m *Matrix3) InvertInPlaceDet(det float64)

+ func (m *Matrix3) Mul(m1 *Matrix3) *Matrix3

+ func (m *Matrix3) MulColumn(c Vec3) Vec3

+ func (m *Matrix3) MulColumnInv(c Vec3, det float64) Vec3

+ func (m *Matrix3) Rows() [3]Vec3

+ func (m *Matrix3) SVD(u, s, v *Matrix3)

+ func (m *Matrix3) Scale(s float64)

+ func (m *Matrix3) Transpose() *Matrix3

+ type Matrix4 [16]float64

+ func NewMatrix4Columns(c1, c2, c3, c4 Vec4) *Matrix4

+ func NewMatrix4Identity() *Matrix4

+ func (m *Matrix4) Add(m1 *Matrix4) *Matrix4

+ func (m *Matrix4) CharPoly() Polynomial

+ func (m *Matrix4) Cols() [4]Vec4

+ func (m *Matrix4) Det() float64

+ func (m *Matrix4) Mul(m1 *Matrix4) *Matrix4

+ func (m *Matrix4) MulColumn(v Vec4) Vec4

+ func (m *Matrix4) Rows() [4]Vec4

+ func (m *Matrix4) SVD(u, s, v *Matrix4)

+ func (m *Matrix4) Scale(s float64) *Matrix4

+ func (m *Matrix4) Sub(m1 *Matrix4) *Matrix4

+ func (m *Matrix4) Transpose() *Matrix4

+ type Vec2 [2]float64

+ func NewVec2RandomNormal() Vec2

+ func NewVec2RandomUnit() Vec2

+ func (v Vec2) Add(v1 Vec2) Vec2

+ func (v Vec2) Dist(v1 Vec2) float64

+ func (v Vec2) Dot(v1 Vec2) float64

+ func (v Vec2) Norm() float64

+ func (v Vec2) Normalize() Vec2

+ func (v Vec2) ProjectOut(v1 Vec2) Vec2

+ func (v Vec2) Scale(f float64) Vec2

+ func (v Vec2) Sub(v1 Vec2) Vec2

+ func (v Vec2) Sum() float64

type Vec3

+ func LeastSquares3(a []Vec3, b []float64, epsilon float64) Vec3

+ func LeastSquaresReg3(a []Vec3, b []float64, lambda, epsilon float64) Vec3

+ func NewVec3RandomUnit() Vec3

+ func (v Vec3) Cross(v1 Vec3) Vec3

+ func (v Vec3) Dot(v1 Vec3) float64

+ func (v Vec3) Norm() float64

+ func (v Vec3) Normalize() Vec3

+ func (v Vec3) OrthoBasis() (Vec3, Vec3)

+ func (v Vec3) ProjectOut(v1 Vec3) Vec3

+ type Vec4 [4]float64

+ func NewVec4RandomNormal() Vec4

+ func NewVec4RandomUnit() Vec4

+ func (v Vec4) Add(v1 Vec4) Vec4

+ func (v Vec4) Dist(v1 Vec4) float64

+ func (v Vec4) Dot(v1 Vec4) float64

+ func (v Vec4) Norm() float64

+ func (v Vec4) Normalize() Vec4

+ func (v Vec4) OrthoBasis() (Vec4, Vec4, Vec4)

+ func (v Vec4) ProjectOut(v1 Vec4) Vec4

+ func (v Vec4) Scale(f float64) Vec4

+ func (v Vec4) Sub(v1 Vec4) Vec4

+ func (v Vec4) Sum() float64

Jan 2, 2022

Changes in this version

+ type KMeans struct

+ Centers []Vec

+ Data []Vec

+ func NewKMeans(data []Vec, numCenters int) *KMeans

+ func (k *KMeans) Assign(vecs []Vec) []int

+ func (k *KMeans) Iterate() float64

type Vec

+ func (v Vec) Dist(v1 Vec) float64

+ func (v Vec) DistSquared(v1 Vec) float64

Nov 25, 2021

Changes in this version

type Polynomial

+ func (p Polynomial) Scale(c float64) Polynomial

Jul 28, 2021

Changes in this version

+ const DefaultGSSIters

+ func GSS(min, max float64, iters int, f func(float64) float64) float64

+ type Polynomial []float64

+ func (p Polynomial) Add(p1 Polynomial) Polynomial

+ func (p Polynomial) Derivative() Polynomial

+ func (p Polynomial) Eval(x float64) float64

+ func (p Polynomial) IterRealRoots(f func(x float64) bool)

+ func (p Polynomial) Mul(p1 Polynomial) Polynomial

+ func (p Polynomial) RealRoots() []float64

+ func (p Polynomial) String() string

+ type SparseCholesky struct

+ func NewSparseCholesky(mat *SparseMatrix) *SparseCholesky

+ func (a *SparseCholesky) ApplyInverseVec3(x []Vec3) []Vec3

+ func (a *SparseCholesky) ApplyVec3(x []Vec3) []Vec3

+ type SparseMatrix struct

+ func NewSparseMatrix(size int) *SparseMatrix

+ func (a *SparseMatrix) ApplyVec3(x []Vec3) []Vec3

+ func (a *SparseMatrix) Iterate(row int, f func(col int, x float64))

+ func (a *SparseMatrix) Permute(perm []int) *SparseMatrix

+ func (a *SparseMatrix) RCM() []int

+ func (a *SparseMatrix) Set(row, col int, x float64)

+ type Vec []float64

+ func (v Vec) Add(v1 Vec) Vec

+ func (v Vec) Dot(v1 Vec) float64

+ func (v Vec) Norm() float64

+ func (v Vec) NormSquared() float64

+ func (v Vec) Normalize() Vec

+ func (v Vec) ProjectOut(v1 Vec) Vec

+ func (v Vec) Scale(s float64) Vec

+ func (v Vec) Sub(v1 Vec) Vec

+ type Vec3 [3]float64

+ func NewVec3RandomNormal() Vec3

+ func (v Vec3) Add(v1 Vec3) Vec3

+ func (v Vec3) Dist(v1 Vec3) float64

+ func (v Vec3) Scale(f float64) Vec3

+ func (v Vec3) Sub(v1 Vec3) Vec3

+ func (v Vec3) Sum() float64