Hacktoberfest 2026：メンテナが10月に向けて印を付けた、オープンで初心者向けの issue。 Hacktoberfest の issue を見る

Slice.to_global puts the second-to-last node row on the domain's upper faces

オープン 初心者向け
  • 何か問題がありますか？
#123 コメント 0 件 リアクション 0 件 担当者 0 名 GitHub で見る

@chraibi がすでに取り組んでいます。

2026年10月6日 から。

  • #124 @chraibi による — オープン

評価

難易度
2/5
見積もり時間
1〜3時間
初心者へのやさしさ
82/100
issue の種類
バグ
明瞭さ
明確に書かれている
活発さ
活発
技術スタック
python
領域
data

調査の方向性

slcf/slice.py の「We ignore border points」ブロック（1.11.7 では 614–638 行目）を読み、Slice.to_global() 用に提供された再現を実行してください。完了とは、ノードベースのスライスが全球最大の面およびコーナーにおける FDS が最後に書き込んだ値を保持し、セル中心のスライスは変更されないことを意味します。

索引モデルが issue の本文から書いたものです。

説明

Slice.to_global() (slcf/slice.py, the "We ignore border points ..." block;
1.11.7 lines 614-638, unchanged on the default branch at c5aac2a) writes the wrong values
into the node row at the global maximum of a node-based slice. For each axis it
first drops the mesh's last node (reduced_shape[axis + 1] -= 1) and only then
saves the "border points" from the shortened array, so the node appended at the
global maximum is the original second-to-last node. Cell-centred slices skip
this block and are not affected. Same code block as #121 (dropped upper-face nodes), distinct defect.

Reproduction (one mesh [0, 1] x [0, 1], node spacing 0.5, horizontal
node-based slice, FDS data K = x + 10 y; the mesh and subslice are fdsreader's
own classes, only the data is set):

import numpy as np
from fdsreader.fds_classes.mesh import Mesh
from fdsreader.slcf.slice import Slice, SubSlice
from fdsreader.utils.extent import Extent

class OneMeshSlice:
    to_global = Slice.to_global
    subslices = Slice.subslices
    orientation, id, n_t, times, cell_centered = 3, "s", 1, np.array([0.0]), False

    def __init__(self):
        nodes = {"x": np.array([0.0, 0.5, 1.0]), "y": np.array([0.0, 0.5, 1.0]),
                 "z": np.array([0.0, 1.0, 2.0])}
        mesh = Mesh(nodes, {"x": (0.0, 1.0), "y": (0.0, 1.0), "z": (0.0, 2.0)}, "M")
        sub = SubSlice(self, "", None, Extent(0.0, 1.0, 0.0, 1.0, 1.0, 1.0), mesh)
        X, Y = np.meshgrid(nodes["x"], nodes["y"], indexing="ij")
        sub._data = (X + 10 * Y)[None].astype(np.float32)  # as FDS wrote it, (t, x, y)
        self._subslices = {"M": sub}
        self.extent = Extent(0.0, 1.0, 0.0, 1.0, 1.0, 1.0)

print(np.asarray(OneMeshSlice().to_global())[0])

FDS data (rows over x): to_global() returns:
[[0. 5. 10. ] [[0. 5. 5. ]
[0.5 5.5 10.5] [0.5 5.5 5.5]
[1. 6. 11. ]] [0.5 5.5 5.5]]

The y = 1 entries repeat y = 0.5, the x = 1 row repeats x = 0.5, and the corner
(1, 1) is 5.5 (the node (0.5, 0.5)) instead of 11. Reproduced on 1.11.5,
1.11.6, 1.11.7, 1.11.9, 1.12.0 and 1.12.1. The same setup with a cell-centred
slice returns the data unchanged.

Expected: the node row at the global maximum is the last node row FDS wrote.

On a real deck (L-shaped corridor, 8 meshes at 0.25 m, SOOT EXTINCTION
COEFFICIENT at z = 2.0 m, FDS 6.10.1) 153 nodes on the upper faces y = 28 m
and x = 45 m are affected; the values differ from the written last node by up
to 0.33 1/m.

Suggested fix: save the border slice (temp_data, and temp_mask when
masked) from the unreduced array, i.e. before slc_data is shortened, and
append it after the repeat as now.

Found through fdsvismap, which reads slices with to_global().

主要言語
Python
スター
75
フォーク
28
平均マージ
2分
マージ済み PR（30日）
2

環境構築

はじめの一歩

  1. issue を最後まで読み、次にプロジェクトのコントリビューションガイドを読みます。
  2. 着手することを issue にコメントします — 二人が同じ作業をするのを防げます。
  3. リポジトリをフォークし、ブランチを切って変更します。
  4. issue 番号を参照したプルリクエストを送ります。

FireDynamics/fdsreader のほかの issue

FireDynamics/fdsreader の issue をすべて見る

似ている issue

Python の issue をもっと見る

新しい issue をメールで受け取る

初心者向けの GitHub issue を短くまとめたダイジェスト。