pimath 0.0.131 → 0.0.133

package/lib/maths/coefficients/nthRoot.ts

@@ -1,148 +0,0 @@
-/**
- * NthRoot is something like "a+b\sqrt{3}
- */
-export class NthRoot {
-    private _radical: number;
-    private _nth: number;
-    private _coefficient: number;
-    private _isValid:boolean;
-
-    constructor(...values: number[]) {
-        this._radical = 1;
-        this._coefficient = 1;
-        this._nth = 2;
-        this._isValid = true;
-
-        if(values !== undefined){
-            this.parse(values[0], values[1], values[2])
-        }
-    }
-
-    // ------------------------------------------
-    // Getter and setter
-    // ------------------------------------------
-    get radical(): number {
-        return this._radical;
-    }
-
-    set radical(value: number) {
-        this._radical = value;
-    }
-
-    get nth(): number {
-        return this._nth;
-    }
-
-    set nth(value: number) {
-        if (Number.isSafeInteger(value) && value >= 2) {
-            this._nth = value;
-        } else {
-            // Error setting the nth root.
-            console.log('Error setting the nth root');
-            this._nth = 2;
-        }
-    }
-
-    get coefficient(): number {
-        return this._coefficient;
-    }
-
-    set coefficient(value: number) {
-        this._coefficient = value;
-    }
-
-    get tex(): string {
-        let C: string;
-
-        if (this._coefficient === 1) {
-            C = '';
-        } else if (this._coefficient === -1) {
-            C = '-';
-        } else {
-            C = this._coefficient.toString();
-        }
-
-        if (this._radical === 1) {
-            return `${this._coefficient}`;
-        } else {
-            if (this._nth === 2) {
-                return `${C}\\sqrt{${this._radical}}`
-            } else {
-                return `${C}\\sqrt[${this._nth}]{${this._radical}}`
-            }
-        }
-    }
-
-    get display(): string {
-        let C: string;
-
-        if (this._coefficient === 1) {
-            C = '';
-        } else if (this._coefficient === -1) {
-            C = '-';
-        } else {
-            C = this._coefficient.toString();
-        }
-
-        if (this._radical === 1) {
-            return `${this._coefficient}`;
-        } else {
-            if (this._nth === 2) {
-                return `${C}sqrt{${this._radical}}`
-            } else {
-                return `${C}root(${this._nth}){${this._radical}}`
-            }
-        }
-    }
-
-    get value(): number {
-        return this._coefficient * Math.pow(this._radical, 1 / this._nth);
-    }
-
-    // ------------------------------------------
-    // Creation / parsing functions
-    // ------------------------------------------
-    parse = (radical: number, nthroot?: number, coefficient?: number): NthRoot => {
-        this._coefficient = (coefficient === undefined) ? 1 : coefficient;
-        this._nth = (nthroot === undefined) ? 2 : nthroot;
-        this._radical = (radical === undefined) ? 1 : radical;
-
-        if(this._nth%2===0 && this._radical<0){
-            this._isValid = false;
-        }
-        return this;
-    };
-
-    // ------------------------------------------
-    // Mathematical operations
-    // ------------------------------------------
-    reduce = (): NthRoot => {
-        // Max value to test.
-        let V = Math.floor(Math.pow(this._radical, 1 / this._nth));
-        while (V > 1) {
-            if (this._radical % Math.pow(V, this._nth) === 0) {
-                // It's dividable by V^n
-                this._coefficient *= V;
-                this._radical = this._radical / Math.pow(V, this._nth);
-
-                // Redifine the new testing value (this is optimization)
-                V = Math.floor(Math.pow(this._radical, 1 / this._nth));
-                continue;
-            }
-            V--;
-        }
-        return this;
-    };
-
-    multiply = (N: NthRoot): NthRoot => {
-        this._radical *= N.radical;
-        return this.reduce();
-    };
-
-    // ------------------------------------------
-    // Help functions
-    // ------------------------------------------
-    hasRadical = ():boolean => {
-        return !(this._radical===1 || this._radical===0 || this._isValid===false)
-    };
-}

package/lib/maths/geometry/circle.ts

@@ -1,379 +0,0 @@
-import {Point} from "./point";
-import {Line, LinePropriety} from "./line";
-import {Vector} from "./vector";
-import {Triangle} from "./triangle";
-import {Numeric} from "../numeric.ts";
-import {Fraction} from "../coefficients/fraction";
-import {Equation} from "../algebra/equation";
-import {Polynom} from "../algebra/polynom";
-import {Monom} from "../algebra/monom";
-
-export class Circle {
-    private _center: Point;
-    private _squareRadius: Fraction;
-    private _cartesian: Equation;
-    private _exists: boolean;
-
-    constructor(...values: unknown[]) {
-        this._exists = false
-
-        if (values !== undefined) {
-            this.parse(...values)
-        }
-    }
-
-    get center(): Point {
-        return this._center;
-    }
-
-    get squareRadius(): Fraction {
-        return this._squareRadius
-    }
-
-    get cartesian(): Equation {
-        return this._cartesian
-    }
-
-    get exists(): boolean {
-        return this._exists;
-    }
-
-    get radius(): { tex: string, display: string, value: number } {
-        if (this._squareRadius.isSquare()) {
-            return {
-                tex: this._squareRadius.clone().sqrt().tex,
-                display: this._squareRadius.clone().sqrt().display,
-                value: this._squareRadius.clone().sqrt().value
-            }
-        } else {
-            return {
-                tex: `\\sqrt{${this._squareRadius.tex}}`,
-                display: `sqrt(${this._squareRadius.display})`,
-                value: this._squareRadius.clone().sqrt().value
-            }
-        }
-        return this._squareRadius
-    }
-
-    get tex(): string {
-        if (this._exists) {
-            let cx, cy
-            if (this._center.x.isZero()) {
-                cx = 'x^2'
-            } else {
-                cx = `\\left(x${this._center.x.isNegative() ? '+' : '-'}${this._center.x.clone().abs().tex}\\right)^2`
-            }
-            if (this._center.y.isZero()) {
-                cy = 'y^2'
-            } else {
-                cy = `\\left(y${this._center.y.isNegative() ? '+' : '-'}${this._center.y.clone().abs().tex}\\right)^2`
-            }
-            return `${cx}+${cy}=${this._squareRadius.tex}`
-        } else {
-            return `\\text{le cercle n'existe pas.}`
-        }
-    }
-
-    get developed(): string {
-        return this._cartesian.tex
-    }
-
-    get display(): string {
-        if (this._exists) {
-            let cx, cy
-            if (this._center.x.isZero()) {
-                cx = 'x^2'
-            } else {
-                cx = `(x${this._center.x.isNegative() ? '+' : '-'}${this._center.x.clone().abs().tex})^2`
-            }
-            if (this._center.y.isZero()) {
-                cy = 'y^2'
-            } else {
-                cy = `(y${this._center.y.isNegative() ? '+' : '-'}${this._center.y.clone().abs().tex})^2`
-            }
-            return `${cx}+${cy}=${this._squareRadius.display}`
-        } else {
-            return `\\text{le cercle n'existe pas.}`
-        }
-    }
-
-    /**
-     * Get the relative position between circle and line. It corresponds to the number of intersection.
-     * @param {Line} L
-     * @returns {number}
-     */
-    relativePosition = (L: Line): number => {
-        let distance = L.distanceTo(this.center), radius = Math.sqrt(this._squareRadius.value)
-
-        if (distance.value - radius > 0.0000000001) {
-            return 0 // external
-        } else if (Math.abs(distance.value - radius) < 0.0000000001) {
-            return 1 // tangent
-        } else {
-            return 2 // external
-        }
-    }
-
-    lineIntersection = (L: Line): Point[] => {
-        let intersectionPoints: Point[] = [], solX: Fraction
-
-        if (this._cartesian === null) {
-            return []
-        }
-        const equX = this._cartesian.clone(), lineX = L.equation.clone().isolate('x'),
-            lineY = L.equation.clone().isolate('y')
-
-        if (lineX instanceof Equation && lineY instanceof Equation) {
-            equX.replaceBy('y', lineY.right).simplify()
-            equX.solve()
-
-            for (let x of equX.solutions) {
-                if (x.exact === false && isNaN(x.value)) {
-                    continue
-                }
-
-                solX = new Fraction(x.exact === false ? x.value : x.exact)
-                intersectionPoints.push(new Point(solX.clone(), lineY.right.evaluate(solX)))
-            }
-        }
-
-        return intersectionPoints
-    }
-
-    tangents = (P: Point | Fraction): Line[] => {
-        if (P instanceof Fraction) {
-            return this._tangentsWithSlope(P)
-        } else if (this.isPointOnCircle(P)) {
-            return this._tangentsThroughOnePointOnTheCircle(P)
-        } else if (this.center.distanceTo(P).value > this.radius.value) {
-            //TODO:  Must check it's outside the circle
-            return this._tangentsThroughOnePointOutsideTheCircle(P)
-        } else {
-            console.log('No tangents as the point is inside !')
-        }
-        return []
-    }
-
-    isPointOnCircle = (P: Point): Boolean => {
-        return this._cartesian.test({x: P.x, y: P.y})
-    }
-
-    getPointsOnCircle = (numberIsInteger?: boolean): Point[] => {
-        if (numberIsInteger === undefined) {
-            numberIsInteger = false
-        }
-
-        // It means searching for pythagorician triples that make a perfect square.
-        // (x-4)^2 + (y+3)^2 = 15
-
-        let triplets = Numeric.pythagoricianTripletsWithTarget(this._squareRadius.value, true)
-
-        let points: Point[] = [], pt
-        triplets.forEach(triplet => {
-            // Allow positive / negative values
-            // x-a = t  => x = a + t
-            // x-a = -t => x = a - t
-
-            for (let k of [[1, 1], [-1, 1], [-1, -1], [1, -1]]) {
-                pt = new Point(
-                    this.center.x.clone().add(k[0] * triplet[0]),
-                    this.center.y.clone().add(k[1] * triplet[1])
-                )
-                // Check if the point is not already in points.
-                if (!pt.isInListOfPoints(points)) {
-                    points.push(pt)
-                }
-            }
-        })
-        return points
-    }
-
-    clone(): Circle {
-        this._center = this._center.clone()
-        this._squareRadius = this._squareRadius.clone()
-        this._calculateCartesian()
-        return this
-    }
-
-    private _tangentsThroughOnePointOnTheCircle = (P: Point): Line[] => {
-        let CT = new Vector(this._center, P)
-        return [new Line(P, CT, LinePropriety.Perpendicular)]
-    }
-
-    private _tangentsThroughOnePointOutsideTheCircle = (P: Point): Line[] => {
-        // y = mx + h
-        // px, py => h = -m px + py => mx - y -m.px + py = 0 =>
-        // Centre: cx, cy, radius: r
-        // (m.cx - cy -m.px + py)^2 = r^2  * (m^2  + 1)
-        // (m(cx-py) - (cy - py))^2 = r^2  * (m^2  + 1)
-
-        let cx_px = this.center.x.clone().subtract(P.x), cy_py = this.center.y.clone().subtract(P.y),
-            polyLeft = new Polynom('x'), polyRight = new Polynom('x^2+1')
-
-        polyLeft.multiply(cx_px).subtract(cy_py).pow(2)
-        polyRight.multiply(this.squareRadius)
-
-        let equ = new Equation(polyLeft, polyRight)
-        equ.moveLeft().simplify().solve()
-
-        return equ.solutions.map(sol => {
-            //  h = -m px + py
-            let h, equ = new Equation('y', 'x')
-
-            if (sol.exact instanceof Fraction) {
-                h = P.x.clone().opposed().multiply(sol.exact).add(P.y)
-                equ.right.multiply(sol.exact).add(h)
-            } else {
-                h = P.x.clone().opposed().multiply(sol.value).add(P.y)
-                equ.right.multiply(sol.value).add(h)
-            }
-
-            return new Line(equ)
-        })
-
-    }
-
-    private _tangentsWithSlope = (slope: Fraction): Line[] => {
-        // d(C;t)=r => ac1+bc2 + x = +- sqrt(a^2 + b^2)*r
-        // x = -ac1-bc2  +-  sqrt(a^2 + b^2)*r
-        // y = a/bx + h => ax-by + H = 0
-
-        const a = slope.numerator, b = -slope.denominator, c1 = this._center.x.clone(), c2 = this._center.y.clone(),
-            r = this._squareRadius
-
-        let sq = this._squareRadius.clone().multiply(slope.numerator ** 2 + slope.denominator ** 2),
-            x1 = c1.clone().multiply(a).opposed().subtract(c2.clone().multiply(b)).add(sq.clone().sqrt()),
-            x2 = c1.clone().multiply(a).opposed().subtract(c2.clone().multiply(b)).subtract(sq.clone().sqrt())
-
-        return [new Line(a, b, x1), new Line(a, b, x2)]
-    }
-
-    private _reset(): Circle {
-        this._center = null
-        this._squareRadius = null
-        this._cartesian = null
-        this._exists = false
-
-        return this
-    }
-
-    private parse(...values: unknown[]): Circle {
-        // Data can be given in these formats:
-        // one value, a string -> make it an Equation
-        // one value, an Equation
-        // one value, a circle -> clone it
-        // two values: two points (center and pointThrough)
-        // two values: point and Fraction (center and radius)
-        // three values: Point, Fraction, Boolean (center, square radius, true)
-
-        this._reset()
-
-        if (typeof values[0] === 'string') {
-            this._parseEquation(new Equation(values[0]))
-        } else if (values[0] instanceof Equation) {
-            this._parseEquation(values[0])
-        } else if (values[0] instanceof Circle) {
-            this._parseCopyCircle(values[0])
-        } else if (values[0] instanceof Point && values.length > 1) {
-            if (values[1] instanceof Point) {
-                if (values[2] instanceof Point) {
-                    this._parseThroughtThreePoints(values[0], values[1], values[2])
-                } else {
-                    this._parseCenterAndPointThrough(values[0], values[1])
-                }
-            } else if (values[1] instanceof Fraction || typeof values[1] === 'number') {
-                this._parseCenterAndRadius(values[0], values[1], (typeof values[2] === "boolean") ? values[2] : false)
-            }
-        }
-
-        // Calculate once the different values.
-        if (this._exists) {
-            this._calculateCartesian()
-
-            // If the square radius is zero or positive, the circle exists.
-            if (this._squareRadius !== undefined && this._squareRadius.isNegative()) {
-                this._exists = false
-            }
-        }
-
-        return this
-    }
-
-    private _calculateCartesian() {
-        this._cartesian = (new Equation(new Polynom(`(x-(${this._center.x.display}))^2+(y-(${this._center.y.display}))^2`), new Polynom(`${this._squareRadius.display}`))).moveLeft()
-    }
-
-    private _parseCopyCircle(circle: Circle): Circle {
-        this._center = circle.center.clone()
-        this._squareRadius = circle.squareRadius.clone()
-        this._calculateCartesian()
-        this._exists = circle.exists
-        return this
-    }
-
-    private _parseCenterAndRadius(center: Point, radius: Fraction | number, square?: boolean): Circle {
-        this._center = center.clone()
-
-        if (square) {
-            this._squareRadius = (new Fraction(radius))
-        } else {
-            this._squareRadius = new Fraction(radius).pow(2)
-        }
-
-        this._exists = true
-        return this
-    }
-
-    private _parseCenterAndPointThrough(center: Point, pointThrough: Point): Circle {
-        this._center = center.clone()
-        this._squareRadius = new Vector(this._center, pointThrough).normSquare
-        this._exists = true
-        return this
-    }
-
-    private _parseEquation(equ: Equation): Circle {
-        this._exists = false
-
-        // Move everything to the left.
-        equ.moveLeft()
-
-        if (equ.degree('x').value === 2 && equ.degree('y').value === 2) {
-            // Both must be of degree 2.
-            let x2 = equ.left.monomByDegree(2, 'x'), y2 = equ.left.monomByDegree(2, 'y'), x1: Monom, y1: Monom, c: Monom
-
-            // Both square monoms must have the same coefficient.
-            if (x2.coefficient.isEqual(y2.coefficient)) {
-                equ.divide(x2.coefficient)
-
-                x1 = equ.left.monomByDegree(1, 'x')
-                y1 = equ.left.monomByDegree(1, 'y')
-
-                c = equ.left.monomByDegree(0)
-
-                this._center = new Point(x1.coefficient.clone().divide(2).opposed(), y1.coefficient.clone().divide(2).opposed())
-
-                this._squareRadius = c.coefficient.clone().opposed()
-                    .add(this._center.x.clone().pow(2))
-                    .add(this._center.y.clone().pow(2))
-
-                this._calculateCartesian()
-                this._exists = true
-            } else {
-                // The circle is not a valid circle
-                this._center = null
-                this._squareRadius = null
-                this._exists = false
-            }
-        }
-        return this
-    }
-
-    private _parseThroughtThreePoints(A: Point, B: Point, C: Point): Circle {
-        let T = new Triangle(A, B, C), mAB = T.remarquables.mediators.AB.clone(),
-            mAC = T.remarquables.mediators.AC.clone()
-        this.parse(mAB.intersection(mAC).point, A)
-
-        return this
-    }
-
-}