inline double PSParallelCompact::normal_distribution(double density) {

  {- -------------------------------------------
  (1) 正規分布における確率を計算するだけ (以下の計算式は教科書通りの正規分布の式).

      (なお, この正規分布の期待値(_dwl_mean), 及び標準偏差(_dwl_std_dev)は
       それぞれ ParallelOldDeadWoodLimiterMean オプション
       及び ParallelOldDeadWoodLimiterStdDev オプションで変更できる.
       See: PSParallelCompact::initialize_dead_wood_limiter())
      ---------------------------------------- -}

      assert(_dwl_initialized, "uninitialized");
      const double squared_term = (density - _dwl_mean) / _dwl_std_dev;
      return _dwl_first_term * exp(-0.5 * squared_term * squared_term);
    }