patrick-kidger / patrick-kidger/optimistix

Least squares performance

Open
#88 2 comments 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

question
Dominant language
Python
Stars
623
Forks
54
PR merge metrics
No merged PRs in 30d

Description

Hi Patrick, I see you are using general solvers to solve the least squares problem, by contructing an objective out of residuals. I have found this approach is less performant than using Levenberg Marquardt. Do you get good performance?

Contributor guide

Open the contributing guide

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

Research direction

Start in optimistix/_least_squares.py at line 43, then compare the current residual-objective approach with the proposed Levenberg–Marquardt approach. Before coding, establish the expected performance comparison; done should be a decided, documented approach or a concrete LM implementation with evidence that it improves least-squares performance.

Written by the indexing model from the issue text.

Assessment

Tech stack
python
Domain
performance
Issue type
Feature
Difficulty
5/5
Estimated time
Over a week
Activity status
Stale
Clarity
Needs clarification
Newbie friendliness
25/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.