[論文レビュー] The Cosmic Linear Anisotropy Solving System (CLASS) I: Overview

Julien Lesgourgues|arXiv (Cornell University)|Apr 14, 2011

Cosmology and Gravitation Theories参考文献 9被引用数 489

ひとこと要約

CLASSは、宇宙論的摂動理論のための新しい、使いやすく、非常に高精度なボルツマンコードであり、宇宙背景放射（CMB）の非等方性および物質パワー スペクトルを、よりモジュラーで、エラー処理が強く、柔軟性に富んだ設計で計算することを目的としている。構造化されたエラー報告とモジュラー設計により、高い精度を達成しており、CAMB や CMBFAST よりも高速であるため、現代の宇宙論におけるパラメータ推定に最適である。

ABSTRACT

The Cosmic Linear Anisotropy Solving System (CLASS) is a new accurate Boltzmann code, designed to offer a more user-friendly and flexible coding environment to cosmologists. CLASS is very structured, easy to modify, and offers a rigorous way to control the accuracy of output quantities. It is also incidentally a bit faster than other codes. In this overview, we present the general principles of CLASS and its basic structure. We insist on the friendliness and flexibility aspects, while accuracy, physical approximations and performances are discussed in a series of companion papers.

研究の動機と目的

• 既存のコード（CAMB や CMBFAST など）よりも使いやすさ、柔軟性、精度制御の面で優れた新しいボルツマンコードの開発。
• CMBFAST や CMBEASY のようなレガシーコードの限界に対処すること。これらのコードはもはやメンテナンスされておらず、レコビリアリゼーションのアルゴリズムも古くなっている。
• 特にプラランク後のデータに対して、絶対的な精度を検証できる独立した、自己完結的なコードベースを提供すること。
• 1回の呼び出しあたりの計算時間を短縮することで、パラメータ空間の探索を効率的に行えるようにすること。
• 質量のあるニュートリノ、非コールドダークマター、テンソルモード、レンズ効果を含む複雑な宇宙論的モデルを、拡張可能なモジュラー設計によりサポートすること。

提案手法

• 明確な責任分担を持つ、マクロベースのエラー処理システム（class_call、class_test）を用いた、モジュラーなC言語アーキテクチャの実装。エラー報告が堅牢かつトレース可能である。
• CMB非等方性の計算に、CMBFAST と同様のラインオブサイト統合法を採用するが、数値的安定性と精度制御を向上させた。
• 人間が読みやすい .ini ファイルと .pre ファイルを介して、宇宙論的パラメータと精度設定を入力可能。初期値をデフォルトとして設定し、即時の起動を可能にした。
• 複数のゲージ選択（v1.0では同期ゲージが利用可能。ニュートンゲージは開発中）をサポート。アディアバティックモードおよび等曲率モードの柔軟な初期条件処理を可能にした。
• OpenMP を介した並列化を、1つのコンパイルフラグで有効化可能。マルチコアシステムでのパフォーマンス向上を実現した。
• 非線形モジュールに正規化アルゴリズムを埋め込み、将来的に HALOFIT などの非線形近似の実装を可能にする。

実験結果

リサーチクエスチョン

• RQ1 CAMB や CMBFAST よりも著しく使いやすく、保守性に優れた新しいボルツマンコードをどのように設計できるか？
• RQ2 宇宙論的シミュレーションコードにおいて、高い精度とトレース可能なエラー報告を保証するための、どのようなアーキテクチャ的・プログラミング的実践が有効か？
• RQ3 既存のコードと同等またはそれ以上の性能を達成しつつ、より複雑な宇宙論的モデルをサポートできる新しいコードは可能か？
• RQ4 精度設定を一元化・キャリブレーションすることで、さまざまな宇宙論的状況において信頼性のある精度制御を実現できるか？
• RQ5 モジュラーで拡張可能なコードベースは、将来の拡張（例：宇宙線のパワー スペクトル、改善された非線形構造形成モデル）をどの程度サポートできるか？

主な発見

• CLASS v1.0 は、CMB TT、TE、EE、BB、およびレンズ効果を考慮したスペクトル、物質パワー スペクトル、複数の種類の物質に対する伝達関数の計算をサポートしている。
• コードには質量のあるニュートリノ、調整可能な位相空間分布を持つ非コールドダークマターの残余物、オプションでテンソルモードと宇宙定数が含まれている。
• エラーメッセージはマクロを介して自動生成され、関数呼び出しの間でチェーン化され、詳細な行番号付きのデバッグ情報を提供する。
• 主なコンパイラと互換性があり、OpenMP を介した並列実行をサポートしており、マルチコアシステムでのパフォーマンス向上を実現している。
• モジュラー設計により、ユーザーは再コンパイルなしに、入力ファイルで特定のモジュール（例：等曲率モード、テンソルモード）を有効または無効にできる。
• CosmoMC や MultiNest、CosmoPMC といったパラメータ抽出ツールとの統合を想定しており、宇宙論における効率的なベイズ推論を可能にしている。