Coloring for symmetric banded matrices
@amontoison 已经在做这个了。
开始于 2026年1月30日。
评估
这个 Issue 还没有评估数据。
描述
The code is quite short but it took me a long a time to understand how we can do the decompression directly from the peridic coloring.
We only need to solve bidiagonal systems for the decompression of each 2-colored tree.
It also allows to not require a buffer like in the acyclic coloring.
using BandedMatrices
T = Float64
n = 10
kd = 3
k = 2 * kd + 1
A = brand(T,n,n,kd,kd)
A = A + A'
color = Int[mod(j,kd+1)+1 for j=1:n]
B = zeros(T, n, kd+1)
for j = 1:n
cj = color[j]
B[:,cj] .+= A[:,j]
end
C = copy(A)
fill!(C, 0)
for i = 1:n
ci = color[i]
C[i,i] = B[i, ci]
end
for j = 1:kd
for i = j+1:kd+1
posi = i
posj = j
c = color[posi]
C[posi,posj] = B[posj,c]
C[posj,posi] = C[posi,posj]
flag = true
while flag
if posi > posj
posi, posj = posj, posi
end
if posi + kd + 1 <= n
c = color[posi + kd + 1]
C[posi+kd+1,posj] = B[posj,c] - C[posj,posi]
C[posj,posi+kd+1] = C[posi+kd+1,posj]
posi = posi + kd + 1
else
flag = false
end
end
end
end
- 主要语言
- Julia
- 星标
- 46
- 派生
- 14
- 平均合并
- 6 小时 30 分钟
- 30 天内合并 PR
- 2
环境准备
这个项目没有提供开发容器、Dockerfile 或贡献指南,环境需要你自己搭建:先看它的 README,通用步骤见我们的新手贡献指南。
从这里开始
- 先读完整个 Issue,再读项目的贡献指南。
- 在 Issue 下留言说明你要接手 —— 这能避免两个人做同样的事。
- Fork 仓库,在一个分支上完成修改。
- 提交 Pull Request,并在描述里引用这个 Issue 编号。
JuliaDiff/SparseMatrixColorings.jl 的其他 Issue
-
难度 3/5 1-2 天 新手友好度 35/100
JuliaDiff/SparseMatrixColorings.jl#342 · 1 条评论 ·
-
Contributors未关闭wontfix
难度 5/5 一周以上 新手友好度 20/100
JuliaDiff/SparseMatrixColorings.jl#331 · 4 条评论 ·
-
Compiled (trimmed) binaries可能已有人在做 @amontoison 于 43 天前认领。 未关闭
难度 5/5 一周以上 新手友好度 35/100
JuliaDiff/SparseMatrixColorings.jl#329 · 2 条评论 · 2 个 reaction ·
-
难度 2/5 1-3 小时 新手友好度 50/100
JuliaDiff/SparseMatrixColorings.jl#301 · 1 条评论 ·
-
Decompress into a single triangle可能已有人在做 @amontoison 于 291 天前认领。 未关闭epic feature
难度 5/5 一周以上 新手友好度 20/100
JuliaDiff/SparseMatrixColorings.jl#280 · 4 条评论 · 1 个 reaction ·
查看 JuliaDiff/SparseMatrixColorings.jl 的全部 Issue
相似的 Issue
-
难度 2/5 1-3 小时 新手友好度 78/100
WaterLily-jl/WaterLily.jl#367 ·
维护者通常 1 天内回复
-
R_to_gamma and R_to_v are wrong for |Z| != 1 (e.g. heavy ions)可能已有人在做 @DavidSagan 今天认领。 未关闭
难度 2/5 1-3 小时 新手友好度 88/100
bmad-sim/BeamTracking.jl#366 ·
维护者通常 1 天内回复
-
TSDAE's docstring implies a theta argument it does not take可能已有人在做 @singhharsh1708 于 1 天前认领。 未关闭documentation
难度 1/5 1 小时以内 新手友好度 74/100
SciML/PETScDiffEq.jl#282 ·
维护者通常 1 天内回复
-
enhancement good first issue help wanted
难度 2/5 半天 新手友好度 82/100
-
SavedSubsystem with an integer save_idxs maps the index the wrong way round可能已有人在做 @ChrisRackauckas 今天认领。 未关闭
难度 2/5 1-3 小时 新手友好度 78/100
SciML/SciMLBase.jl#1643 ·
维护者通常 2 天内回复