Skip to content
Snippets Groups Projects

Draft: Fix sparse lbm

Open Markus Holzer requested to merge holzer/lbmpy:UpdateAPISparse into master
2 unresolved threads

This MR updates the sparse LBM by adding correct boundary handling and renewal of the API.

Merge request reports

Loading
Loading

Activity

Filter activity
  • Approvals
  • Assignees & reviewers
  • Comments (from bots)
  • Comments (from users)
  • Commits & branches
  • Edits
  • Labels
  • Lock status
  • Mentions
  • Merge request status
  • Tracking
153 161
154 162 class SparseLbBoundaryMapper:
155 163 NEIGHBOR_IDX_NAME = 'nidx{}'
156 DIR_SYMBOL = TypedSymbol("dir", np.uint32)
164 DIR_SYMBOL = TypedSymbol("dir", np.int64)
  • 203 if fa.field == f_out:
    204 f_dir = 'out'
    205 elif fa.field == f_in:
    206 f_dir = 'in'
    207 else:
    208 continue
    209 i = 1 if f_dir == "in" else 2
    210 accessor_subs[fa] = pdf_field_sparse.absolute_access((index_field(i),), ())
    211
    212 self.boundary_eqs = assignments.new_with_substitutions(accessor_subs)
    213 self.boundary_eqs_orig = assignments
    189 214 self.method = method
    190 215 self.index_field_dtype = index_field_dtype
    191 216 self.neighbor_offsets = neighbor_offsets
    192 217 self.index_field = index_field
    218 self.timesptep_indexing = indexing
  • Looks okay overall, but I am not familiar enough with the sparse LB to review this in detail.

  • Michael Kuron removed review request for @kuron

    removed review request for @kuron

  • Markus Holzer added 53 commits

    added 53 commits

    Compare with previous version

  • Please register or sign in to reply
    Loading