Samsara\Fermat\Provider\SequenceProvider::nthBernoulliNumber()

Complexity

Conditions	10
Paths	21

Size

Total Lines	89
Code Lines	58

Duplication

Lines	0
Ratio	0 %

Code Coverage

Tests	50
CRAP Score	10.0056

Importance

Changes	2
Bugs	0	Features	1

<?php
    public static function nthBernoulliNumber($n, ?int $scale = null): DecimalInterface
    {

        $scale = $scale ?? 5;

        $internalScale = (int)ceil($scale*(log10($scale)+1));

        $n = Numbers::makeOrDont(Numbers::IMMUTABLE, $n, $internalScale)->setMode(CalcMode::Precision);

        if (!$n->isWhole()) {
            throw new IntegrityConstraint(
                'Only integers may be indexes for Bernoulli numbers',
                'Ensure only integers are provided as indexes',
                'An attempt was made to get a Bernoulli number with a non-integer index'
            );
        }

        if ($n->isLessThan(0)) {
            throw new IntegrityConstraint(
                'Index must be non-negative',
                'Provide only non-negative indexes for Bernoulli number generation',
                'An attempt was made to get a Bernoulli number with a negative index'
            );
        }

        if ($n->isEqual(0)) {
            return Numbers::makeOne($scale);
        }

        if ($n->isEqual(1)) {
            return Numbers::make(Numbers::IMMUTABLE, '-0.5', $scale);
        }

        if ($n->modulo(2)->isEqual(1)) {
            return Numbers::makeZero($scale);
        }

        $tau = Numbers::makeTau($internalScale)->setMode(CalcMode::Precision);

        $d = Numbers::make(Numbers::IMMUTABLE, 4, $internalScale)->setMode(CalcMode::Precision)->add(
            $n->factorial()->ln($internalScale)->subtract(
                $n->multiply($tau->log10($internalScale))
            )->truncate()

The method truncate() does not exist on Samsara\Fermat\Types\Base\Interfaces\Numbers\FractionInterface.

The method truncate() does not exist on Samsara\Fermat\Types\Base\Interfaces\Numbers\NumberInterface. It seems like you code against a sub-type of Samsara\Fermat\Types\Base\Interfaces\Numbers\NumberInterface such as Samsara\Fermat\Types\Base\Interfaces\Numbers\DecimalInterface or Samsara\Fermat\Types\Decimal.

        )->add(
            $n->numberOfIntDigits()
        );
        $s = $d->multiply(
            Numbers::makeNaturalLog10($internalScale)
        )->multiply(
            '0.5'
        )->divide($n, $internalScale)->exp($internalScale)->truncate()->add(1);

The method exp() does not exist on Samsara\Fermat\Types\Base\Interfaces\Numbers\FractionInterface.

The method exp() does not exist on Samsara\Fermat\Types\Base\Interfaces\Numbers\NumberInterface. It seems like you code against a sub-type of Samsara\Fermat\Types\Base\Interfaces\Numbers\NumberInterface such as Samsara\Fermat\Types\Base\Interfaces\Numbers\DecimalInterface or Samsara\Fermat\Types\Decimal.

        $internalScale = ($d->isGreaterThan($internalScale)) ? $d->asInt() : $internalScale;

The method asInt() does not exist on Samsara\Fermat\Values\ImmutableFraction.

The method asInt() does not exist on Samsara\Fermat\Values\MutableFraction.

The method asInt() does not exist on Samsara\Fermat\Types\Fraction.

The method isGreaterThan() does not exist on Samsara\Fermat\Types\Base\Interfaces\Numbers\NumberInterface. It seems like you code against a sub-type of said class. However, the method does not exist in Samsara\Fermat\Types\Base\Number. Are you sure you never get one of those?

The method asInt() does not exist on Samsara\Fermat\Types\Base\Interfaces\Numbers\FractionInterface.

The method asInt() does not exist on Samsara\Fermat\Types\Base\Interfaces\Numbers\NumberInterface. It seems like you code against a sub-type of Samsara\Fermat\Types\Base\Interfaces\Numbers\NumberInterface such as Samsara\Fermat\Types\Base\Interfaces\Numbers\DecimalInterface or Samsara\Fermat\Types\Decimal.

        $s = $s->truncateToScale($internalScale);

The method truncateToScale() does not exist on Samsara\Fermat\Types\Base\Interfaces\Numbers\NumberInterface. It seems like you code against a sub-type of Samsara\Fermat\Types\Base\Interfaces\Numbers\NumberInterface such as Samsara\Fermat\Types\Base\Interfaces\Numbers\DecimalInterface or Samsara\Fermat\Types\Decimal.

The method truncateToScale() does not exist on Samsara\Fermat\Types\Base\Interfaces\Numbers\FractionInterface.

        $n = $n->truncateToScale($internalScale);
        $tau = Numbers::make2Pi($internalScale)->setMode(CalcMode::Precision);
        $p = Numbers::makeOne($internalScale)->setMode(CalcMode::Precision);
        $t1 = Numbers::makeOne($internalScale)->setMode(CalcMode::Precision);
        $t2 = Numbers::makeOne($internalScale)->setMode(CalcMode::Precision);

        while ($p->isLessThanOrEqualTo($s)) {
            $p = self::_nextprime($p);
            $pn = $p->pow($n);
            $pn1 = $pn->subtract(1);
            $t1 = $pn->multiply($t1);
            $t2 = $pn1->multiply($t2);
        }

        $z = $t1->divide($t2, $internalScale);
        $oz = Numbers::makeZero($internalScale)->setMode(CalcMode::Precision);

        while (!$oz->isEqual($z)) {
            $oz = $z;
            $p = self::_nextprime($p);
            $pn = $p->pow($n);
            $pn1 = $z->divide($pn, $internalScale);
            $z = $z->add($pn1);
        }

        $f = $n->factorial();
        $taun = $tau->pow($n);

        $z = $z->multiply(2)->multiply($f)->divide($taun, $internalScale);

        if ($n->modulo(4)->isEqual(0)) {
            $z = $z->multiply(-1);
        }

        return $z->round($scale);

The method round() does not exist on Samsara\Fermat\Types\Base\Interfaces\Numbers\FractionInterface.

The method round() does not exist on Samsara\Fermat\Types\Base\Interfaces\Numbers\NumberInterface. It seems like you code against a sub-type of Samsara\Fermat\Types\Base\Interfaces\Numbers\NumberInterface such as Samsara\Fermat\Types\Base\Interfaces\Numbers\DecimalInterface or Samsara\Fermat\Types\Decimal.

The method round() does not exist on Samsara\Fermat\Types\Fraction.

}