Questions tagged [linear-algebra]
The linear-algebra tag has no usage guidance.
11
questions
0
votes
1
answer
24
views
Implementing Preconditioned conjugate gradient
I have implemented the Preconditioned Conjugate Gradient (PCG) method for solving a system of linear equations in Python and I would appreciate it if someone could verify its correctness since I am ...
2
votes
1
answer
156
views
Python integer optimization with PuLP involving polynomials that I think can be faster
I have a problem where we are given multivariate polynomials in variables like y_i and z_j and the goal is to express them as sums of products of terms of the form y_i-z_j with nonnegative integer ...
0
votes
0
answers
18
views
Efficiently computing a batch of results given a batch assignment vector and series of corresponding matrices
I have a 1D tensor of tokens that belong to different batches. The batch sizes here are uneven. Each batch needs to be multiplied with a corresponding weight matrix. My current approach is using a ...
4
votes
2
answers
442
views
Generic matrix library in Java
I have this repository. The most important source files follow.
com.github.coderodde.math.linear.matrix.AbstractMatrix.java:
...
2
votes
0
answers
53
views
Implementing PLU Decomposition in Rust
I am using abstalg, because I wanted to learn more about how to do this with generic fields, (So you can implement LU decomposition not just for floats). I used ...
2
votes
1
answer
75
views
An implementation of Givens method for diagonalising symmetric matrices
So I tried to create a code which aims to diagonalise symmetric matrices via Givens algorithm. The program only works for symmetric matrices, for we know a symmetric matrix is always diagonalisable, ...
2
votes
1
answer
298
views
Matrix-vector multiplication with block-column distribution and MPI_Reduce_scatter
I'm studying MPI and I'm trying to solve the exercise (number 3.5) taken from Pacheco's book Introduction to parallel programming. Essentially the exercise is:
Write an MPI program that computes the ...
4
votes
1
answer
165
views
Yet another mathematical vector class
This is an implementation of a mathematical vector.
I am currently working on a linear algebra project and before proceeding further with other classes (matrix, polar_vector etc.) I would like to have ...
2
votes
2
answers
216
views
Coordinate Descent Non-negative Least Squares optimization
In my fast implementation of Non-negative Matrix Factorization (Rcpp Machine Learning Library, RcppML), about 40% of the runtime is spent solving Non-negative Least Squares (NNLS) systems (the rest of ...
5
votes
1
answer
191
views
Why is this tensor contraction so slow and how can I make it faster?
I am trying to calculate a specific complicated tensor contraction (of complex valued tensors) which appears as part of some other calculation and is called multiple times in a row, about 100 times or ...
1
vote
0
answers
103
views
Solving a linear system with Gaussian Elimination
I decided to implement a solver for linear systems of equations based on the gaussian elimination and reduction to upper triangular form. The gaussian elimination is quite simple to implement but the ...