Parallelization of a concentrated load using GridapDistributed
まだ誰も着手していません。
評価
- 難易度
- 4/5
- 見積もり時間
- 3〜5日
- 初心者へのやさしさ
- 25/100
- issue の種類
- バグ
- 明瞭さ
- 説明が足りない
- 活発さ
- 停滞
- 技術スタック
- julia
調査の方向性
分散 local_views(model) の例と issue で報告された _point_to_cell!(cache, p) のアサーションから始め、その後この経路をシリアルな DiracDelta 実装と比較してください。点がアクティブなローカルセル内で見つからない理由を特定し、集中荷重に対してサポートされている並列アプローチを文書化または実装し、完了を示す再現可能な線形弾性の例を用意してください。
索引モデルが issue の本文から書いたものです。
説明
I'm trying to parallelize a linear elastic problem solver with a concentrated load applied on the middle point of ine boundary. In series code the load is described using DiracDelta function as following:
function DiracDelta(model::DiscreteModel{D}, pvec::Vector{Point{D,T}}) where {D,T}
trian = Triangulation(model)
cell_to_pindices = _cell_to_pindices(pvec,trian)
cell_ids = collect(keys(cell_to_pindices))
cell_points = collect(values(cell_to_pindices))
points = map(i->pvec[cell_points[i]], 1:length(cell_ids))
weights_x_cell = Fill.(one(T),length.(cell_points))
pquad = map(i -> GenericQuadrature(points[i],weights_x_cell[i]), 1:length(cell_ids))
trianv = view(trian,cell_ids)
pmeas = Measure(CellQuadrature(pquad,points,weights_x_cell,trianv,PhysicalDomain(),PhysicalDomain()))
GenericDiracDelta{0,D,NotGridEntity}(trianv,pmeas)
end
Since in parallel mode, the model is a GenericDistributedDiscreteModel() which is not supported in original definition, I tried to write a parallel version of it:
tmp = map(local_views(model)) do model
p = Point(1.0,0.5)
trian = Triangulation(model)
cache = _point_to_cell_cache(KDTreeSearch(),trian)
cell = _point_to_cell!(cache, p)
...
end
Yet the _point_to_cell!(cache,p) throughs an error: ERROR: AssertionError: Point (1.0, 0.5) is not inside any active cell
I'm quite sure that point p is in the grid and the serial code runs well. Any idea how this problem results from? Or is there any other parallel application of a concentrated load on a point or a line?
- 主要言語
- Julia
- スター
- 150
- フォーク
- 28
- PR マージ指標
- 30日以内にマージされた PR はありません
環境構築
このプロジェクトには開発コンテナ、Dockerfile、コントリビューションガイドがありません。まず README を読み、一般的な手順ははじめてのコントリビューションガイドを参照してください。
はじめの一歩
- issue を最後まで読み、次にプロジェクトのコントリビューションガイドを読みます。
- 着手することを issue にコメントします — 二人が同じ作業をするのを防げます。
- リポジトリをフォークし、ブランチを切って変更します。
- issue 番号を参照したプルリクエストを送ります。
gridap/GridapDistributed.jl のほかの issue
-
A generic method of the abstract type DistributedDiscreteModel{Dc,Dp} directly accesses to a field name .face_gids再び着手できるかも @amartinhuertas が 121 日前に担当しましたが、オープン中のプルリクエストはありません。 オープン
gridap/GridapDistributed.jl#211 · 担当者 1 名 ·
-
Automatically add ghost cells to triangulation whenever a CellField is created out of a DistributedTriangulation再び着手できるかも @JordiManyer が 154 日前に担当しましたが、オープン中のプルリクエストはありません。 オープン
gridap/GridapDistributed.jl#209 · 担当者 2 名 ·
-
難易度 4/5 3〜5日 初心者へのやさしさ 35/100
gridap/GridapDistributed.jl#189 · コメント 1 件 ·
-
Issue with DG-SIP and FullyAssembledRows再び着手できるかも @amartinhuertas が 359 日前に担当しましたが、オープン中のプルリクエストはありません。 オープン
gridap/GridapDistributed.jl#186 · コメント 20 件 · 担当者 1 名 ·
-
Potential issue when interpolating DistributedCelllFields that are only defined on own cells再び着手できるかも @JordiManyer が 373 日前に担当しましたが、オープン中のプルリクエストはありません。 オープン
gridap/GridapDistributed.jl#184 · コメント 1 件 · 担当者 2 名 ·
gridap/GridapDistributed.jl の issue をすべて見る
似ている issue
-
難易度 2/5 1〜3時間 初心者へのやさしさ 72/100
-
documentation invalid
難易度 1/5 1〜3時間 初心者へのやさしさ 88/100
JuliaPhysics/BeamletOptics.jl#190 ·
メンテナーはふだん 2 日以内に返信
-
mul! with a Symmetric or Hermitian factor throws for a matrix view whose columns are not contiguousオープン
難易度 2/5 1〜3時間 初心者へのやさしさ 72/100
JuliaLang/LinearAlgebra.jl#1829 ·
メンテナーはふだん 1 日以内に返信
-
`EnsembleSummary` `ci_type = :SEM` ribbon is `sqrt(N)` too wide for scalar states対応中かも @isaacsas が 1 日前に担当しました。 オープン
難易度 2/5 1〜3時間 初心者へのやさしさ 88/100
SciML/SciMLBase.jl#1654 ·
メンテナーはふだん 3 日以内に返信
-
難易度 2/5 1〜3時間 初心者へのやさしさ 78/100
MagneticResonanceImaging/MRIReco.jl#320 ·
メンテナーはふだん 1 日以内に返信