超材料电磁分析
理论与物理
什么是超材料
超材料是指具有自然界所没有的特性的材料吗？能实现隐形斗篷之类的吗？
它是人工实现负折射率和负磁导率的周期性结构体。粗略地说，通过将远小于波长的金属图案（如SRR或细线阵列）规则排列，可以产生自然界材料无法实现的电磁响应。
诶，那隐形斗篷真的能造出来吗？
完全的隐形斗篷还只是理论上的事情。Pendry（2006）的变换光学理论虽然优美，但要实现宽带、低损耗是极其困难的。不过，已经实用化的应用倒是不少：
- 天线小型化：利用负磁导率缩短天线的电长度（例：手机基站）
- EMC对策的EBG结构：利用电磁带隙抑制印刷电路板的杂散辐射
- FSS（频率选择表面）：用于雷达罩、隐身技术的频率滤波
- 超透镜：超越衍射极限的成像（已在毫米波频段验证）
嘿，原来已经有实用级别的东西了。但是用CAE怎么分析呢？感觉和普通的FEM不一样…
问到了核心。超材料分析有3个特点：
- 用单元胞的周期边界条件，仅用一个胞元就能模拟无限周期结构
- 用Drude-Lorentz模型描述频率相关的等效介电常数和磁导率
- 用S参数反演法从仿真结果反算等效物性值
我们按顺序来看。
Drude-Lorentz模型
Drude-Lorentz模型是什么？光听名字就觉得很难…
这是一种将超材料的电磁响应表示为频率函数的色散模型。有效介电常数可以用Lorentz型表示：
同样地，对于像SRR（开口谐振环）这样的磁谐振结构，有效磁导率也是Lorentz型：
$\omega_{pe}$和$\omega_{0e}$这些表示什么？
| 参数 | 含义 | 物理背景 |
|---|---|---|
| $\omega_{pe}$ | 电等离子体频率 | 由细线阵列的几何形状决定 |
| $\omega_{0e}$ | 电谐振频率 | Drude型中$\omega_{0e}=0$ |
| $\gamma_e$ | 电损耗系数 | 取决于金属电导率和表面粗糙度 |
| $\omega_{0m}$ | 磁谐振频率 | 由SRR的间隙电容和自感决定 |
| $F$ | 填充率 | SRR的面积占有率 |
| $\gamma_m$ | 磁损耗系数 | 金属电阻和辐射损耗的总和 |
原来如此，意思是在谐振频率附近介电常数或磁导率会变成负的吗？
没错。在$\omega_{0m} < \omega < \omega_{pm}$（磁等离子体频率）的频带内，$\mu_{\text{eff}} < 0$。如果同时存在$\varepsilon_{\text{eff}} < 0$的频带与之重叠，那么折射率$n = \sqrt{\varepsilon_{\text{eff}} \mu_{\text{eff}}}$就会变成负的——这就是Veselago（1968）预言的“左手系超材料”。
超材料中的Maxwell方程
在超材料里面，Maxwell方程的形式会改变吗？
形式本身不变。但决定性区别在于本构关系（$\mathbf{D} = \varepsilon(\omega)\mathbf{E}$、$\mathbf{B} = \mu(\omega)\mathbf{H}$）具有频率色散性。在时域FDTD中，通过ADE（辅助微分方程）法或RC电路等效模型来纳入色散性。
这里$k_0 = \omega/c$是自由空间波数。关键在于$\varepsilon_r(\omega)$和$\mu_r(\omega)$是遵循Drude-Lorentz模型的复数。实部为负的区域，波不会传播而成为倏逝波，或者作为左手系的后向波传播。
Bloch-Floquet周期边界条件
超材料是假设无限周期结构排列的吧？但计算机无法模拟无限个吧？
问得好。这时就用Bloch-Floquet周期边界条件。只建模一个单元胞，在相对面上给电场施加相位差：
$\mathbf{a}$是晶格矢量，$\mathbf{k}$是Bloch波数矢量。让这个$\mathbf{k}$沿着布里渊区的不可约区域（如$\Gamma \to X \to M \to \Gamma$）扫描，就能得到色散关系（能带图）。存在带隙的频率带正是EBG（电磁带隙）。
布里渊区在固体物理里听过。电磁学里也是同样的概念啊！
正是如此。光子晶体和超材料在数学上与固体物理的能带理论是同一框架。区别在于，光子晶体的周期与波长同量级，而超材料的周期远小于波长（$a \ll \lambda$），因此可以作为等效的连续介质用ε、μ来描述。
S参数反演法（NRW法）
从仿真得到的S参数，怎么求出等效的ε、μ呢？
使用Nicolson-Ross-Weir（NRW）法的扩展版。从超材料平板的$S_{11}$（反射）和$S_{21}$（透射）反算阻抗$Z$和折射率$n$：
求出$n$和$Z$后，等效物性值就很简单了：
太巧妙了！但感觉好像有陷阱…
敏锐。最大的陷阱是分支选择的模糊性。$e^{jnk_0 d}$的反函数会出现log函数，所以$n$的虚部存在$2m\pi/(k_0 d)$的任意整数$m$的不确定性。当平板厚度$d$超过波长的一半时，选择正确的分支就变得不简单了。对策我们会在故障排除章节详细讨论。
超材料的诞生——“负折射率”从理论到实证花了40年
同时具有负介电常数和负磁导率的“超材料”概念由Veselago（1968）从理论上预言，但实证要等到2000年的Smith等人。实验中组合金属环和细线阵列，在微波频段观测到了负折射。Pendry等人的“完美透镜”理论（2000）和变换光学隐身理论（2006）拥有超过千次的引用数，为电磁学历史翻开了新的一页。
各项的物理意义
- 等离子体频率 $\omega_{pe}$：在细线阵列中由$\omega_{pe}^2 = 2\pi c^2 / (a^2 \ln(a/r))$决定。$a$是晶格常数，$r$是细线半径。细线越细，等离子体频率越低，因此即使在GHz频段也能实现“金属性”（$\varepsilon < 0$）响应。
- 磁谐振频率 $\omega_{0m}$：对于SRR，可以用LC电路模型估算为$\omega_{0m} = 1/\sqrt{LC}$。由间隙电容$C$和环的自感$L$决定，因此几何形状的优化直接关系到谐振频率的设计。
- 损耗项 $\gamma$：由于其实部与$\varepsilon < 0$的频带重叠，因此决定性支配着超材料的性能（品质因数：$\text{FOM} = |\text{Re}(n)| / \text{Im}(n)$）。设计上选择铜或铝的Q值高的频段很重要。
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 角频率 $\omega$ | rad/s | $\omega = 2\pi f$。GHz频段下$\omega \sim 10^{10}$ rad/s |
| 波数 $k_0$ | rad/m | $k_0 = \omega/c$。10 GHz下$k_0 \approx 209$ rad/m |
| 晶格常数 $a$ | m | 超材料中$a \ll \lambda$。典型情况$a < \lambda/10$ |
| 电导率 $\sigma$ | S/m | 铜：5.96×10⁷、铝：3.77×10⁷ |
| S参数 | 无量纲（复数） | dB表示：$|S_{21}|_{\text{dB}} = 20\log_{10}|S_{21}|$ |
数值解法与实现
FEM边单元法进行单元胞分析
超材料的FEM分析，用结构分析一样的节点单元可以吗？
绝对不行。高频电磁场分析如果不使用边单元（Nedelec单元/Whitney单元），会产生大量伪模式（非物理的虚假解）。边单元能自动保证电场切向分量的连续性，排除不满足$\nabla \cdot \mathbf{D} = 0$的伪解。
原来如此，和结构分析那样在节点上插值温度或位移完全不同呢。
是的。边单元的自由度位于单元的边上，直接插值矢量场的切向分量。弱形式如下：
$\mathbf{N}_i$是边单元的基函数。将其离散化后得到广义特征值问题：
$[S]$相当于旋度-旋度刚度矩阵，$[T]$相当于质量矩阵。将周期边界条件结合进去，并以Bloch波数$\mathbf{k}$为参数进行扫描，就能得到能带图。
FDTD法进行超材料分析
FDTD也能分析超材料吗？色散材料在时域里好像很麻烦…
CST Studio Suite的时域求解器正是基于FDTD的。色散材料用ADE（辅助微分方程）法处理。对于Drude-Lorentz模型，引入辅助变量$\mathbf{P}$（极化电流）：
这个辅助微分方程在Yee网格上显式地进行时间积分。FDTD的优点在于一次计算就能得到宽频带的S参数。入射高斯脉冲并获取时间响应，进行FFT就能立刻得到S11、S21的频率特性。
FEM和FDTD，该用哪个好呢？
| 比较项目 | FEM（频域） | FDTD（时域） | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 能带图计算 | 直接得到（特征值问题） | 间接（需要多次计算） | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 宽频带S参数 | 需要频率扫描 | 一次计算即可获取 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 色散材料 | 每个频率代入ε(ω) | 需要ADE法的辅助变量 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 曲面形状精度 | 高（四面体网格） | 低（阶梯近似） | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 内存效率 | 稀疏矩阵，效率高 | 整个空间用网格填充，较大<
ベンダーの系譜と製品統合の経緯各ソフトの成り立ちって、結構ドラマチックだったりしますか？ Ansys HFSS次はAnsys HFSSの話ですね。どんな内容ですか？ CST Studio SuiteCST Studioって、具体的にはどういうことですか？ Computer Simulation Technology (ドイツ) が開発。2016年にDassault Systèmesが買収しSIMULIAに統合。 現在の所属: Dassault Systèmes SIMULIA COMSOL Multiphysics「COMSOL Multiphysics」について教えてください！ 1986年スウェーデンで設立。MATLAB連携のFEMLABとして開始、後にCOMSOLに改名。マルチフィジックスに強み。 現在の所属: COMSOL AB 待って待って、が開発したってことは、つまりこういうケースでも使えますか？ ファイル形式と相互運用性異なるソルバー間でモデルを変換する際は、要素タイプの対応関係、材料モデルの互換性、荷重・境界条件の表現差異に注意が必要になるんだ。特に高次要素や特殊要素（コヒーシブ要素、ユーザー定義要素等）はソルバー間で直接変換できない場合が多い。 なるほど…フォーマットって一見シンプルだけど、実はすごく奥が深いんですね。 実務上の注意点教科書には載ってない「現場の知恵」みたいなものってありますか？ メッシュ収束性の確認、境界条件の妥当性検証、材料パラメータの感度分析がすごく大事なんだ。
メタマテリアル電磁解析の全体像がつかめました！ 明日から実務で意識してみます。 うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。 Coffee Break よもやま話
メタマテリアルの誕生——「負の屈折率」は理論から実証まで40年かかった負の誘電率・透磁率を同時に持つ「メタマテリアル」の概念はVeselago（1968）が理論的に予言したが、実証は2000年のSmith et al.まで待つことになった。実験では金属リングとワイヤのアレイを組み合わせ、マイクロ波帯で負の屈折を観測した。「ものが消える」完全マント（Cloaking）装置の可能性は世界的な注目を集め、Pendryらの理論（2006）は千件超の被引用数を誇る。CAE解析では周期境界条件とSパラメータ反転法によりメタマテリアルの等価定数を抽出する手法が標準だ。 各項の物理的意味
仮定条件と適用限界
数値解法と実装数値手法の詳細具体的にはどんなアルゴリズムでメタマテリアル電磁解析を解くんですか？ 離散化の定式化形状関数 $N_i$ を用いて未知量を近似: $$ u^h(\mathbf{x}) = \sum_{i=1}^{n} N_i(\mathbf{x}) \, u_i $$
これを数式で表すとこうなるよ。 $$ K_e = \int_{\Omega_e} B^T \, D \, B \, d\Omega \approx \sum_{g=1}^{n_g} w_g \, B^T(\xi_g) \, D \, B(\xi_g) \, |J(\xi_g)| $$ 基礎方程式の離散形これを数式で表すとこうなるよ。 $$ n = \sqrt{\varepsilon_r\mu_r} $$
$$ \varepsilon_{eff}(\omega) = 1 - \frac{\omega_p^2}{\omega^2-\omega_0^2+j\gamma\omega} $$
うーん、式だけだとピンとこないです… 何を表してるんですか？ 連続体の支配方程式を離散化すると、以下の代数方程式系が得られる: $$ [K]\{u\} = \{F\} $$
ここで $[K]$ は全体剛性マトリクス（または同等のシステムマトリクス）、$\{u\}$ は未知節点変数ベクトル、$\{F\}$ は外力ベクトルなんだ。 あっ、そういうことか！ 連続体の支配方程式をってそういう仕組みだったんですね。 要素技術「要素技術」って聞いたことはあるんですけど、ちゃんと理解できてないかもしれません…
積分スキーム積分スキームって、具体的にはどういうことですか？ ここまで聞いて、要素タイプがなぜ重要か、やっと腹落ちしました！ 収束性と安定性収束しなくなったら、まず何をチェックすればいいですか？
収束速度: 二次要素で $O(h^2)$ のオーダーで誤差が減少（滑らかな解の場合） なるほど…メッシュを細分化って一見シンプルだけど、実はすごく奥が深いんですね。 ソルバー設定の推奨事項辺要素（Nedelec要素）電磁場解析に特化した要素。接線成分の連続性を自動的に保証し、スプリアスモードを排除。3D高周波解析の標準。 節点要素スカラーポテンシャル定式化に使用。静磁場のスカラーポテンシャル法や静電場解析で有効。 FEM vs BEM（境界要素法）FEM: 非線形材料・非均質媒質に対応。BEM: 無限領域（開領域問題）を自然に扱える。ハイブリッドFEM-BEMも有効。 非線形収束（磁気飽和）B-Hカーブの非線形性をニュートン・ラフソン法で処理。残差基準: $||R||/||R_0|| < 10^{-4}$が一般的。 周波数領域解析時間高調波仮定により定常問題に帰着。複素数演算が必要だが、広帯域特性は時間領域解析で取得。 時間領域の時間刻み最高周波数成分の1/20以下の時間刻みが必要。暗黙的時間積分ではより大きな刻みも可能だが精度に注意。 周波数領域と時間領域の使い分け周波数領域解析は「ラジオの特定の周波数に合わせる」ようなもの——1つの周波数での応答を効率的に計算できる。時間領域解析は「全チャンネルを同時に録画する」ようなもの——あらゆる周波数成分を含む過渡現象を再現できるが計算コストが高い。 実践ガイド実践ガイド先生、「実践ガイド」について教えてください！ メタマテリアル電磁解析の実務的な解析フローと注意点を解説する。 解析フロー最初の一歩から教えてください！ 何から始めればいいですか？ 2. 求解 (Solving)
メッシュ生成のベストプラクティスメッシュの良し悪しってどうやって判断するんですか？ 要素品質指標メッシュ密度の決定メッシュ密度の決定って、具体的にはどういうことですか？
境界条件の設定指針境界条件って、ここを間違えると全部ダメになるって聞いたんですけど…
あっ、そういうことか！ 過拘束に注意ってそういう仕組みだったんですね。 商用ツール別の実装手順いろんなソフトがあるんですよね？ それぞれの特徴を教えてください！
Ansys HFSS次はAnsys HFSSの話ですね。どんな内容ですか？ CST Studio SuiteCST Studioって、具体的にはどういうことですか？ Computer Simulation Technology (ドイツ) が開発。2016年にDassault Systèmesが買収しSIMULIAに統合。 現在の所属: Dassault Systèmes SIMULIA 先生の説明分かりやすい！ ツール名のモヤモヤが晴れました。 よくある失敗と対策初心者がやりがちな失敗パターンってありますか？ 事前に知っておきたいです！
品質保証チェックリスト教科書には載ってない「現場の知恵」みたいなものってありますか？
メタマテリアル電磁解析の全体像がつかめました！ 明日から実務で意識してみます。 うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。 Coffee Break よもやま話
「平らなレンズ」——メタマテリアルが変えたミリ波イメージング負の屈折率メタマテリアルは焦点距離ゼロの「完全レンズ」を理論上実現できる。実用上の壁は「損失」で、金属共鳴構造の虚部成分がGHzから先の周波数で増大し、解像度向上の恩恵を打ち消す。それでもmmW帯（77 GHz, 150 GHz）では損失が管理可能な水準となり、自動車レーダや空港セキュリティスキャナへの応用が試みられている。CAEでは複素誘電率・透磁率のフィッティングモデル（Drude-Lorentz）を用いた時間領域解析が設計の基本となる。 解析フローのたとえモータの電磁界解析は「ギターの調律」に近い感覚です。弦の太さ（コイル巻数）とブリッジの位置（磁石配置）を調整して、最も美しい音色（効率の良いトルク特性）を引き出す。1つのパラメータを変えると全体のバランスが変わる——だからパラメトリックスタディが重要なんです。 初心者が陥りやすい落とし穴「空気領域？ なんで空気をメッシュで切るの？」——初めて電磁界解析に触れた人がほぼ全員抱く疑問です。答えは「磁力線は鉄心の外にも広がるから」。解析領域を鉄心ぎりぎりにすると、行き場を失った磁束が壁に「ぶつかって」反射し、実際にはありえない磁束集中が起きます。部屋が狭すぎてボールが壁に跳ね返りまくる状態を想像してみてください。 境界条件の考え方遠方の境界条件って地味ですが超重要です。「ここから先は無限に広がる空間」ということを数値的に表現する必要がある。設定を間違えると、まるで「見えない壁」があるかのように磁束が跳ね返されてしまいます。 ソフトウェア比較商用ツール比較いろんなソフトがあるんですよね？ それぞれの特徴を教えてください！ メタマテリアル電磁解析に対応する主要な商用CAEツールの機能比較と、各製品の歴史的背景を詳述する。 対応ツール一覧で、メタマテリアル電磁解析をやるにはどんなソフトが使えるんですか？
Ansys HFSS次はAnsys HFSSの話ですね。どんな内容ですか？ CST Studio SuiteCST Studioって、具体的にはどういうことですか？ Computer Simulation Technology (ドイツ) が開発。2016年にDassault Systèmesが買収しSIMULIAに統合。 現在の所属: Dassault Systèmes SIMULIA COMSOL Multiphysics「COMSOL Multiphysics」について教えてください！ 1986年スウェーデンで設立。MATLAB連携のFEMLABとして開始、後にCOMSOLに改名。マルチフィジックスに強み。 現在の所属: COMSOL AB 機能比較マトリクス変換時のリスク変換時のリスクって、具体的にはどういうことですか？
あっ、そういうことか！ 異なるツール間でのモってそういう仕組みだったんですね。 ライセンス形態「ライセンス形態」って聞いたことはあるんですけど、ちゃんと理解できてないかもしれません…
選定の指針結局どれを選べばいいか、判断基準を教えてもらえますか？ メタマテリアル電磁解析のツール選定においては以下を考慮:
メタマテリアル電磁解析の全体像がつかめました！ 明日から実務で意識してみます。 うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。 Coffee Break よもやま話
メタマテリアル解析ツール——CST Studio Suite vs COMSOLメタマテリアル解析にはCST Studio Suite（時間領域FDTD・周波数領域FEM）とCOMSOL Wave Optics/RF Moduleが多く使われる。CSTの「Unit Cell」シミュレーション機能は周期境界条件と平面波入射の設定が容易で、Sパラメータ抽出からバンド図計算までワンストップで対応する。COMSOLは多物理（熱・力学）連成解析が得意で、温度変化による誘電定数変動が特性に与える影響を詳細評価できる。Ansys LSTCのLS-DYNAはEMモジュールでメタマテリアルの動的変形時の特性変化も解析可能だ。 選定で最も重要な3つの問い
先端技術先端トピックと研究動向メタマテリアル電磁解析の分野って、これからどう進化していくんですか？ メタマテリアル電磁解析における最新の研究動向と先進的手法を見ていこう。 最新の数値手法次は最新の数値手法の話ですね。どんな内容ですか？ うーん、式だけだとピンとこないです… 何を表してるんですか？ 高性能計算 (HPC) への対応
|
| 並列化手法 | 概要 | 適用ソルバー |
|---|---|---|
| MPI (領域分割) | 分散メモリ型。大規模問題の標準 | 全主要ソルバー |
| OpenMP | 共有メモリ型。ノード内並列 | 多くのソルバー |
| GPU (CUDA/OpenCL) | GPGPU活用。特に陽解法で有効 | LS-DYNA, Fluent等 |
| ハイブリッド MPI+OpenMP | ノード間+ノード内並列 | 大規模HPC環境 |
トラブルシューティング
トラブルシューティング
常见错误与对策
先生もメタマテリアル電磁解析で徹夜デバッグしたことありますか？（笑）
1. 収束失敗
収束失敗って、具体的にはどういうことですか？
症状: ソルバーが指定反復回数内に収束せず異常終了
考えられる原因:
- メッシュ品質の不足（過度に歪んだ要素）
- 材料パラメータの不適切な設定
- 不適切な初期条件
- 非線形性が強すぎる（荷重ステップの不足）
つまり収束失敗のところで手を抜くと、後で痛い目を見るってことですね。肝に銘じます！
2. 非物理的な結果
次は非物理的な結果の話ですね。どんな内容ですか？
症状: 応力/変位/温度等が物理的に非現実的な値
考えられる原因:
- 境界条件の誤設定
- 単位系の混在（SI単位と工学単位の混同）
- 不適切な要素タイプの選択
- 応力特異点の存在
対策:
- 反力の合計を確認（力の釣り合い）
- 単位系の一貫性を確認
- 要素タイプの適切性を再検討
- 特異点除去またはサブモデリング
先輩が「収束失敗だけはちゃんとやれ」って言ってた意味が分かりました。
3. 計算時間の超過
計算時間の超過って、具体的にはどういうことですか？
症状: 計算が想定時間の何倍もかかる
対策:
- メッシュの粗密分布の最適化
- 対称性の活用（1/2, 1/4モデル）
- ソルバー設定の最適化（反復法、前処理の選択）
- 並列計算の活用
4. メモリ不足
「メモリ不足」について教えてください！
症状: Out of Memory エラー
先輩が「収束失敗だけはちゃんとやれ」って言ってた意味が分かりました。
対策:
- アウトオブコア解法の使用
- メッシュ規模の削減
- 64bit版ソルバーの使用確認
- メモリ割り当ての増加
おお〜、収束失敗の話、めちゃくちゃ面白いです！ もっと聞かせてください。
Nastran代表的エラー
代表的エラーって、具体的にはどういうことですか？
- FATAL 2012: 特異剛性マトリクス → 拘束条件の見直し
- USER WARNING 5291: 要素品質不良 → メッシュ修正
- SYSTEM FATAL 3008: メモリ不足 → MEM設定の調整
Abaqus代表的エラー
「代表的エラー」について教えてください！
- Excessive distortion: 要素の過大変形 → NLGEOM確認、メッシュ改善
- Zero pivot: 拘束不足 → 境界条件追加
- Time increment too small: 収束失敗 → ステップ設定見直し
なるほど。じゃあツール名ができていれば、まずは大丈夫ってことですか？
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——メタマテリアル電磁解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
なった
詳しく
報告