`Matrix(sol)` can return other types than `Matrix`
Nobody has claimed this yet.
Assessment
- Difficulty
- 3/5
- Estimated time
- 1-2 days
- Newbie friendliness
- 30/100
Research direction
Start with src/vector_of_array.jl around lines 90-95 and reproduce Matrix(sol) for a DE solution whose states are CuVectors. Compare the resulting type with Matrix(Matrix(sol)), then inspect nearby array-construction behavior. Done means the intended Matrix conversion behavior is settled and covered by a regression test.
Written by the indexing model from the issue text.
Description
Calling Matrix(sol) on an *DE solution returns a CuMatrix when the states are CuVectors.
This is because Matrix(sol) is implemented using a generic hcat (see snippet below). This is a very obvious and useful thing to do, but would it make more sense as an implementation of convert(AbstractMatrix, sol) rather than Matrix(sol)?
I know constructors are not guaranteed to return the invoked type, but it's confusing that you have to call Matrix(Matrix(sol)) to actually get a Matrix.
- Dominant language
- Julia
- Stars
- 233
- Forks
- 76
- Avg merge
- 2h 17m
- Merged PRs (30d)
- 9
Contributor guide
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
More from SciML/RecursiveArrayTools.jl
-
Difficulty 1/5 Under an hour Newbie friendliness 68/100
SciML/RecursiveArrayTools.jl#255 ·
-
JET: similar_type(::Any) no-matching-method in RecursiveArrayToolsRaggedArrays (pre-existing QA red) Open
Difficulty 4/5 3-5 days Newbie friendliness 48/100
SciML/RecursiveArrayTools.jl#620 ·
-
Difficulty 4/5 3-5 days Newbie friendliness 45/100
SciML/RecursiveArrayTools.jl#619 ·
-
bug
Difficulty 4/5 3-5 days Newbie friendliness 45/100
SciML/RecursiveArrayTools.jl#570 ·
-
Difficulty 4/5 3-5 days Newbie friendliness 28/100
SciML/RecursiveArrayTools.jl#522 ·
All issues in SciML/RecursiveArrayTools.jl
Similar issues
-
Difficulty 2/5 1-3 hours Newbie friendliness 82/100
-
interpolate! from a Float32 LatitudeLongitudeGrid fails on Metal: Float64 division in find_λ_range Open
Difficulty 2/5 1-3 hours Newbie friendliness 88/100
CliMA/Oceananigans.jl#6056 ·
-
Difficulty 1/5 1-3 hours Newbie friendliness 88/100
JuliaPluto/PlutoPlotly.jl#72 ·
-
Difficulty 2/5 1-3 hours Newbie friendliness 65/100
epiforecasts/BVDOutbreakSize#811 ·
-
bug
Difficulty 2/5 1-3 hours Newbie friendliness 85/100
SciML/NonlinearSolve.jl#1310 · 1 comment ·