Erfcx.java

  1. /*
  2.  * Licensed to the Apache Software Foundation (ASF) under one or more
  3.  * contributor license agreements.  See the NOTICE file distributed with
  4.  * this work for additional information regarding copyright ownership.
  5.  * The ASF licenses this file to You under the Apache License, Version 2.0
  6.  * (the "License"); you may not use this file except in compliance with
  7.  * the License.  You may obtain a copy of the License at
  8.  *
  9.  *      http://www.apache.org/licenses/LICENSE-2.0
  10.  *
  11.  * Unless required by applicable law or agreed to in writing, software
  12.  * distributed under the License is distributed on an "AS IS" BASIS,
  13.  * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
  14.  * See the License for the specific language governing permissions and
  15.  * limitations under the License.
  16.  */
  17. package org.apache.commons.numbers.gamma;

  18. /**
  19.  * Scaled complementary error function.
  20.  *
  21.  * <p>\[ \operatorname{erfcx}(z) = \operatorname{erfc}(z)\ e^{z^2} \]
  22.  *
  23.  * <p>For large z the value is approximately:
  24.  *
  25.  * <p>\[ \operatorname{erfcx}(z) = \frac{1}{z \sqrt{\pi}} \]
  26.  *
  27.  * @see Erfc
  28.  * @since 1.1
  29.  */
  30. public final class Erfcx {
  31.     /** Private constructor. */
  32.     private Erfcx() {
  33.         // intentionally empty.
  34.     }

  35.     /**
  36.      * Returns the scaled complementary error function.
  37.      *
  38.      * <p>Special cases:
  39.      * <ul>
  40.      * <li>If the argument is 0, then the result is 1.
  41.      * <li>If the argument is +infinity, then the result is 0.
  42.      * <li>If the argument is negative and {@code exp(x*x)} is infinite, then the result is +infinity.
  43.      * <li>If the argument is nan, then the result is nan.
  44.      * </ul>
  45.      *
  46.      * @param x Value.
  47.      * @return the scaled complementary error function.
  48.      */
  49.     public static double value(double x) {
  50.         return BoostErf.erfcx(x);
  51.     }
  52. }