Macaulay2 / Macaulay2/M2

crash in rawLinAlgSolve over rationals

Open
#2,831 1 comment 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

Engine Linear Algebra
Dominant language
Macaulay2
Stars
435
Forks
297
Avg merge
4d 20h
Merged PRs (30d)
11

Description

```m2
i1 : solve(0 * id_(QQ^1), map(QQ^1, QQ^1, 0))
-- SIGSEGV
-* stack trace, pid: 376853
0# stack_trace(std::ostream&, bool) at /tmp/macaulay2-20230205-22019-789qlq/M2-release-1.21/M2/Macaulay2/bin/main.cpp:127
1# segv_handler at /tmp/macaulay2-20230205-22019-789qlq/M2-release-1.21/M2/Macaulay2/bin/main.cpp:244
2# 0x00007F1FF5C3EA00 in /lib64/libc.so.6
3# MutableMat >::solveLinear(MutableMatrix const*) const at /tmp/macaulay2-20230205-22019-789qlq/M2-release-1.21/M2/Macaulay2/e/mutablemat-imp.hpp:91
4# rawLinAlgSolve at /tmp/macaulay2-20230205-22019-789qlq/M2-release-1.21/M2/Macaulay2/e/interface/mutable-matrix.cpp:844
5# interface_rawLinAlgSolve at /tmp/macaulay2-20230205-22019-789qlq/M2-release-1.21/M2/Macaulay2/d/interface.dd:3112
```

Contributor guide

No contributing guide indexed for this repository

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

Reproduce the exact solve command in Macaulay2 and inspect MutableMat>::solveLinear in M2/Macaulay2/e/mutablemat-imp.hpp. Trace how rawLinAlgSolve in M2/Macaulay2/e/interface/mutable-matrix.cpp reaches that method, then verify the rational zero-matrix case completes without a SIGSEGV.

Written by the indexing model from the issue text.

Assessment

Tech stack
cpp
Domain
backend, tooling
Issue type
Bug
Difficulty
3/5
Estimated time
1-2 days
Activity status
Active
Clarity
Clearly specified
Newbie friendliness
70/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.