@ccpc/math 0.1.0 → 0.1.3

package/types/algorithm/intersect/surfaces_x_complex.d.ts

@@ -0,0 +1,29 @@
+import { Surface } from '../../geometry/surface';
+import { Vec3 } from '../../base/vec3';
+import { Tol } from '../../base/tol';
+import { ISurfacesXInfo } from './x_info';
+import { Curve3 } from '../../geometry/curve3d';
+export declare class SurfacesXComplex {
+    private _surface1;
+    private _surface2;
+    private _surfacePatchs1;
+    private _surfacePatchs2;
+    private _singularCurvePatchs1;
+    private _singularCurvePatchs2;
+    private _tol;
+    constructor(surf1: Surface, surf2: Surface, tol?: Tol);
+    allIntersects(): ISurfacesXInfo[];
+    singleIntersectCurve(refPoint: Vec3, refDir?: Vec3, useHighPrecision?: boolean, convertToNurbs?: boolean): Curve3 | undefined;
+    private _splitSurface;
+    private _getNextPatchPairForOnePtCurve;
+    private _updateConnectPatchPairs;
+    private _getNextPatchPair;
+    private _spliceIntersectCurves;
+    /**
+     * 在intercurve首部或者尾部拼接另一条交线
+     * @ curve 用于拼接的曲线 // 目前仅支持直线
+     * @ tol 容差
+     */
+    private _spliceCurve;
+    private _getOverlapSurfaceIntersection;
+}

package/types/algorithm/intersect/surfaces_x_special.d.ts

@@ -0,0 +1,16 @@
+import { Tol } from '../../base/tol';
+import { Surface } from '../../geometry/surface';
+import { ISurfacesXInfo } from './x_info';
+/**
+ * （内部使用）几何法和解析法计算得到交线，限于一些特殊情况求交
+ * @param surface 曲面
+ * @param surface 曲面
+ * @returns 交线可能不止一条，故返回交线的数组
+ */
+export declare class SurfacesXSpecial {
+    static execute(surface1: Surface, surface2: Surface, tol?: Tol): ISurfacesXInfo[] | undefined;
+    private static _planePlane;
+    private static _planeCylinder;
+    private static _planeEllipseCylinder;
+    private static _cylinderCylinder;
+}

package/types/algorithm/intersect/surfaces_x_util.d.ts

@@ -0,0 +1,41 @@
+import { ISurfacesXInfo, ISurfaceSurfaceIntersectPointInfo } from './x_info';
+import { Surface } from '../../geometry/surface';
+import { Interval } from '../../base/interval';
+import { Tol } from '../../base/tol';
+import { Vec3 } from '../../base/vec3';
+import { types } from '../../type_define/i_types';
+export declare class SurfacesXUtil {
+    private _surface;
+    private _surf1BoundaryUVs;
+    private _surf2BoundaryUVs;
+    private _useHighPrecision;
+    private _tol;
+    constructor(surf1: Surface, surf2: Surface, tol?: Tol, surf1RangeUV?: Interval[], surf2RangeUV?: Interval[]);
+    /**
+     * 标准的surface与surface通用求交的函数，会返回所有的交（交线、交点、重合面（未实现））
+     * 数值方法（迭代法）计算得到交线，通用曲面求交
+     * @returns 交线可能不止一条，故返回交线的数组
+     */
+    calAllIntersects(rangeU1?: Interval, rangeV1?: Interval, rangeU2?: Interval, rangeV2?: Interval, useHighPrecision?: boolean): ISurfacesXInfo[];
+    calSingleIntersect(referPoint: Vec3, referDir?: Vec3, useHighPrecision?: boolean): ISurfacesXInfo | undefined;
+    findSurfaceIntersectPoint(referPoint: Vec3): ISurfaceSurfaceIntersectPointInfo | undefined;
+    calSingleSelfIntersect(refUV1: types.IXY, refUV2: types.IXY, referDir?: Vec3, useHighPrecision?: boolean): ISurfacesXInfo | undefined;
+    private _initBoudaryUVs;
+    private _initSurfacePatch;
+    private _shrinkSurfacePatchRangesByAnotherSurfaceBox;
+    private _isSurfaceRangeUVFinite;
+    private _shrinkSurfacePatchRanges;
+    private _surfacePatchsIntersectPoints;
+    private _combineSurfacePatchPairs;
+    private _FilterSurfacePatchPairs;
+    private _refreshSurfacePatchs;
+    private _calcSurfacePatchPairIntersectPoint;
+    private _removeInitialIntersectPts;
+    private _calTwoSurfaceSingleIntersect;
+    private _estimateNewIntersectPt;
+    private _calcNextInterscetPointInfo;
+    private _estimateSurfaceUV;
+    private _estimateCurveT;
+    private _adjustEstimatedPtByBoundary;
+    private _refineNewIntersectPtIteratively;
+}

package/types/algorithm/intersect/x_info.d.ts

@@ -0,0 +1,65 @@
+import { Curve3 } from '../../geometry/curve3d';
+import { Vec3 } from '../../base/vec3';
+import { Interval } from '../../base/interval';
+import { Vec } from '../../base/vec';
+import { Surface } from '../../geometry/surface';
+import { types } from '../../type_define/i_types';
+import { Vec2 } from '../../base/vec2';
+/**
+ * 曲线与曲线的交点信息
+ */
+export interface ICurvesXInfo<PointType extends Vec> {
+    /** 交点 */
+    point: PointType;
+    /** 交点在第一条曲线上的参数 */
+    param1: number;
+    /** 交点在第二条曲线上的参数 */
+    param2: number;
+    /** 是否重合 */
+    isOverlap: boolean;
+    /** 若是重合，重合段在第一条曲线上的参数 */
+    overlap1?: Interval;
+    /** 若是重合，重合段在第二条曲线上的参数 */
+    overlap2?: Interval;
+    /** 若是重合，两重合段同向时为 true */
+    overlapSameDirection?: boolean;
+}
+export interface ICurvesXInfo2d extends ICurvesXInfo<Vec2> {
+}
+export interface ICurvesXInfo3d extends ICurvesXInfo<Vec3> {
+}
+export interface ICurveSurfXPointInfo {
+    point: Vec3;
+    curveT: number;
+    uvPara: types.IXY;
+}
+/**
+ * 曲线与曲面的求交信息
+ */
+export interface ICvSurfXInfo {
+    /** 交点 */
+    point: Vec3;
+    /** 交点在曲线上的参数 */
+    curveT: number;
+    /** 交点在曲面上的参数uv */
+    surfaceUV: types.IXY;
+    /** 若是重合，返回重合段的交线的参数域 */
+    overlapRange?: Interval;
+}
+export interface ISurfaceSurfaceIntersectPointInfo {
+    point: Vec3;
+    uvPara1: types.IXY;
+    uvPara2: types.IXY;
+    isSingularity?: boolean;
+}
+/**
+ * 曲面与曲面的求交信息
+ */
+export interface ISurfacesXInfo {
+    /** 交线 */
+    curve?: Curve3;
+    /** 交点（仅一点相交） */
+    point?: Vec3;
+    /** 交面（面面重合） */
+    surface?: Surface;
+}

package/types/algorithm/loop_property/loop-area.d.ts

@@ -0,0 +1,49 @@
+import { types } from '../../type_define/i_types';
+import { Curve2 } from '../../geometry/curve2';
+import { Vec3 } from '../../base/vec3';
+import { DiscreteParam } from '../../base/discrete_param';
+import { Curve3 } from '../../geometry/curve3d';
+/**
+ *
+ * 计算二维区域面积
+ */
+export declare class LoopArea {
+    /**
+     * 计算有序点所围成的二维区域面积
+     * http://en.wikipedia.org/wiki/Shoelace_formula
+     * @param pts
+     */
+    static areaOfPoints(pts: types.IXY[]): number;
+    /**
+     * 计算二维曲线所围成的区域面积
+     * https://en.wikipedia.org/wiki/Green%27s_theorem
+     * @param loop
+     */
+    static areaOfLoop(loop: Curve2[]): number;
+    /**
+     * 计算有序点列在投影方向围成的区域面积
+     * @param pts 三维点序列
+     * @param refDir 参考方向向量
+     */
+    static areaOfPoint3ds(points: types.IXYZ[], refDir: Vec3): number;
+    /**
+     * 计算有序点列的区域面积向量。向量方向为多边形最大投影面积方向，长度为该投影面积
+     * @param points
+     */
+    static areaVectorOfPoint3ds(points: types.IXYZ[]): Vec3;
+    /**
+     * 用离散的方法计算有序曲线组的区域面积向量。向量方向为多边形最大投影面积方向，长度为该投影面积
+     * @param curves
+     * @param discreteParam
+     */
+    static areaVectorOfCurve3ds(curves: Curve3[], discreteParam?: DiscreteParam): Vec3;
+    private static _areaOfCurve2d;
+    /**
+     * 返回三角形的面积
+     * 注意：p1 p2 会被用作临时变量而发生变化
+     * @param p1
+     * @param p2
+     * @param p3
+     */
+    private static _areaOfTriangle;
+}

package/types/algorithm/loop_property/loop-centroid.d.ts

@@ -0,0 +1,40 @@
+import { types } from '../../type_define/i_types';
+import { Vec2 } from '../../base/vec2';
+import { Curve2 } from '../../geometry/curve2';
+import { Vec3 } from '../../base/vec3';
+import { DiscreteParam } from '../../base/discrete_param';
+import { Curve3 } from '../../geometry/curve3d';
+/**
+ *
+ * 计算二维区域形心
+ */
+export declare class LoopCentroid {
+    /**
+     * 计算有序点所围成的二维区域的形心
+     * https://en.wikipedia.org/wiki/Centroid#CITEREFBourke1997
+     * @param pts
+     * @param tol
+     */
+    static centroidOfPoints(pts: types.IXY[], tol?: number): Vec2;
+    /**
+     * 计算二维曲线所围成的区域形心
+     * @param loop
+     */
+    static centroidOfLoop(loop: Curve2[], areaEps?: number): Vec2;
+    static centroidOfPoint3ds(pts: types.IXYZ[], areaEps?: number): Vec3;
+    static centroidInfoOfPoint3ds(pts: types.IXYZ[], areaEps?: number): {
+        centroid: Vec3;
+        area: number;
+        areaNormal: Vec3;
+    };
+    static centroidOfCurve3ds(curves: Curve3[], discreteParam?: DiscreteParam, areaEps?: number): Vec3;
+    static centroidInfoOfCurve3ds(curves: Curve3[], discreteParam?: DiscreteParam, areaEps?: number): {
+        centroid: Vec3;
+        area: number;
+        areaNormal: Vec3;
+    };
+    private static _integralForCurve2d;
+    private static _integralForLine2d;
+    private static _integralForPoints;
+    private static _integralForArc2d;
+}

package/types/algorithm/merge_geometry/halfplane.d.ts

@@ -0,0 +1,24 @@
+import { Ln2 } from '../../geometry/ln2';
+import { Vec2 } from '../../base/vec2';
+import { types } from '../../type_define/i_types';
+export declare class HalfPlane {
+    static createByLine2d(line: Ln2): HalfPlane;
+    w: types.IXY;
+    b: number;
+    constructor(w: types.IXY, b: types.IXY | number);
+    get normal(): Vec2;
+    distance(p: types.IXY): number;
+    clone(): HalfPlane;
+    toLine2d(): Ln2;
+    offset(c: number): HalfPlane;
+    intersect(h: HalfPlane): types.IXY;
+    parallel(h: HalfPlane, tol?: number): boolean;
+    dump(): {
+        w: {
+            x: number;
+            y: number;
+        };
+        b: number;
+    };
+    load(data: any): HalfPlane;
+}

package/types/algorithm/merge_geometry/merge_curve.d.ts

@@ -0,0 +1,14 @@
+import { Curve2 } from '../../geometry/curve2';
+import { Ln2 } from '../../geometry/ln2';
+import { IMergeCurveInfo } from '../geometry_merge';
+export declare class MergeCurve {
+    /**
+     * 仅支持 arc 和 line
+     * 默认容差为1e-6
+     */
+    static mergeCurve2d(curves: Curve2[], distol?: number): Curve2[];
+    static mergeCurve2dEx(curves: Curve2[], distol?: number): IMergeCurveInfo;
+    static mergeCurveLine2d(curves: Curve2[], distol: number): Curve2[];
+    static mergeCurveLine2dEx(lines: Ln2[], distol: number): IMergeCurveInfo;
+    static mergeArc2d(curves: Curve2[], eps?: number): IMergeCurveInfo;
+}

package/types/algorithm/merge_geometry/merge_point.d.ts

@@ -0,0 +1,18 @@
+export declare class Disjoint {
+    set: number[];
+    constructor();
+    clear(): void;
+    find(a: number): number;
+    merge(a: number, b: number): void;
+}
+export declare class MergePoint {
+    private _tmp;
+    private _disjoint;
+    constructor();
+    clear(): void;
+    merge(points: ([number, number] | [number, number, number])[], tol?: number): {
+        index: number[];
+        points: ([number, number] | [number, number, number])[];
+    };
+    private _mergeex;
+}

package/types/algorithm/mesh/clip_mesh.d.ts

@@ -0,0 +1,27 @@
+import { types } from '../../type_define/i_types';
+import { Plane } from '../../geometry/plane';
+import { Polygon } from '../../topology/polygon';
+export interface IMeshClipResult {
+    clipped: boolean;
+    aboveMesh: types.IRenderMesh;
+    downMesh?: types.IRenderMesh;
+}
+/**
+ * 平面切Mesh
+ */
+export declare class ClipMesh {
+    /**
+     * 平面切mesh
+     * @param rn mesh
+     * @param plane 平面
+     * @param keepDoubleSide 保留平面两侧的mesh(默认只保留上方 keepDoubleSide = false)
+     * @param deleteCoplanar 删除共面的三角面片(default deleteCoplanar = true)
+     */
+    static clipFaceMesh(rn: types.IRenderMesh, plane: Plane, keepDoubleSide?: boolean, deleteCoplanar?: boolean): IMeshClipResult;
+    static meshPlaneIntersect(rn: types.IRenderMesh, plane: Plane): Polygon;
+    private static _calIntersectInfo;
+    private static _getVertexInfo;
+    private static _cloneRenderMesh;
+    private static _calNormal;
+    private static _reorderPts;
+}

package/types/algorithm/mesh/extrude_clip.d.ts

@@ -0,0 +1,236 @@
+import { Loop, Ln2, Ln3, Matrix4, Box2, Box3, Vec3 } from '../..';
+import { types } from '../../type_define/i_types';
+import { Coord3 } from '../../base/coord3';
+declare enum PositionType {
+    inner = -1,
+    on = 0,
+    outer = 1
+}
+export interface ITypedTriangleVertex {
+    type: PositionType;
+    pos: types.IXYZ;
+}
+interface IVertex {
+    v: types.IXYZ;
+    n: types.IXYZ;
+    uv: types.IXY;
+    uv1?: types.IXY;
+}
+/**
+ * 1. 草图不能包含曲线，如果存在曲线，需事先做离散处理
+ * 2. 草图不能为凹，如果遇到草图为凹的情况，需事先将草图拆分为多个凸多边形处理
+ */
+export declare class Extruder {
+    private _coord;
+    private _loop;
+    private _height;
+    private _localLine3dsCache;
+    private _line2dsCache;
+    private _planeCache;
+    boundingBox2: Box2;
+    curveBoundingBoxList: Box2[];
+    private _boundingBox3;
+    private _lpts;
+    constructor(_c: Coord3, _l: Loop, _h: number);
+    /**
+     * 草图轮廓
+     */
+    get loop(): Loop;
+    /**
+     * 平面原点
+     */
+    get origin(): Vec3;
+    /**
+     * 平面法向
+     */
+    get normal(): Vec3;
+    /**
+     * 三维包围盒
+     * */
+    get localBoundingBox3(): Box3;
+    get worldBoundingBox3(): Box3;
+    get localLine3ds(): Ln3[];
+    get line2ds(): Ln2[];
+    get height(): number;
+    get world2localMatrix(): Matrix4;
+    get local2worldMatrix(): Matrix4;
+    /**
+     * 所有角点在世界坐标系中的位置
+     */
+    get worldCorners(): types.IXYZ[];
+    /**
+     * 所有角点在草图平面局部坐标系中的位置
+     */
+    get localCorners(): types.IXYZ[];
+    /**
+     * 将mesh变换到局部坐标系下
+     */
+    transformWorldMesh2Local(mesh: types.IFlatMeshPlus): types.IFlatMeshPlus;
+    /**
+     * 将本地坐标系下的mesh转换到世界坐标系
+     */
+    transformLocalMesh2World(mesh: types.IFlatMesh): types.IFlatMesh;
+    private transformMesh;
+    getUVAt(pt: types.IXYZ): types.IXY;
+    intersect(mesh: types.IFlatMesh, localcoordinate?: boolean): boolean;
+    edgeBoxInterTri(vertex1: IVertex, vertex2: IVertex, vertex3: IVertex): boolean;
+    /**
+     * 基于pnpoly算法判断点是否在多边形内部
+     * @param point
+     * @returns
+     */
+    inside(point: types.IXY): boolean;
+    log(): void;
+}
+export declare class ExtrudeClip {
+    static tol: number;
+    /**
+     * 入口：针对硬装背景墙定制接口，所有的face以mesh数组的形式传入、传出，但是补面则以独立mesh传出
+     * @param extruder 拉伸体
+     * @param meshes 待裁切的mesh集合
+     * @param option transformed - 是否已经变换到extruder的局部坐标系之下   fill: 是否需要补面
+     * @return clipped - 和传入的mesh一一对应   patch - 补面结果
+     */
+    static multiClip(extruder: Extruder, meshes: types.IFlatMesh[], option?: {
+        transformed?: boolean;
+        fill?: boolean;
+    }): {
+        clipped: types.IFlatMesh[];
+        patch?: types.IFlatMesh;
+    };
+    /**
+     * 入口：拉伸体切割mesh。尽量在外部将mesh转换到extruder的局部坐标系内
+     * @param extruder 拉伸体
+     * @param mesh 被切割mesh对象
+     * @param keepOuter true - 保留拉伸体外部的mesh  false - 保留拉伸体内部的mesh
+     * @param option.transformed true - 输入mesh已经转换至extruder局部坐标系   false - 输入的mesh为世界坐标系，需在内部转换到局部坐标系
+     * @param option.fill true - 补面  false - 不补面
+     * @param option.faceEdgeMap  仅在multiClip时需要传入
+     */
+    static clip(extruder: Extruder, mesh: types.IFlatMeshPlus, option?: {
+        transformed?: boolean;
+        fill?: boolean;
+        faceEdgeMap?: Map<Ln2, Ln3[]>;
+    }): types.IFlatMeshPlus;
+    /**
+     * 合并相同顶点
+     * @param faceMesh 传入facemesh
+     */
+    private static mergeVertex;
+    /**
+     * 对切割后的mesh进行补面，默认mesh在extruder的局部坐标系下
+     * @param extruder 原拉伸体
+     * @param mesh 待切割mesh
+     */
+    private static fillClip;
+    /**
+     * 将三角面片插入到mesh中
+     * @param tri 三角面片
+     * @param normal 法向
+     * @param uvTransform uv矩阵
+     * @param mesh 待插入的mesh
+     */
+    private static insertTriangle;
+    /**
+     * 准备补面数据使用：计算边是否在草图轮廓线上
+     * @param line 局部坐标系下的二维草图轮廓线
+     * @param v1 起点
+     * @param v2 终点
+     * @returns
+     */
+    private static calcLinePlaneEdge;
+    /**
+     * 切割三角面片，得到切割后的三角面片
+     * 待优化：一次计算好三角面片的顶点顺序，避免矫正
+     */
+    private static cutTriangle;
+    /**
+     * 垂直于草图平面的三角面片被拉伸体切割算法
+     * @param extruder 拉伸体
+     * @param v1 顶点1
+     * @param v2 顶点2
+     * @param v3 顶点3
+     * @returns 切割后的三角面片
+     */
+    private static cutVerticalTriangle;
+    /**
+     * 三角面片被平面上的线条代表的拉伸矩形区域切割
+     * @param extruder 拉伸体
+     * @param triPlane 三角面片所在的平面
+     * @param vertex1 顶点1
+     * @param vertex2 顶点2
+     * @param vertex3 顶点3
+     * @param cutLines 代表切割区域的线条
+     * @returns 切割结果
+     */
+    private static cutVerticalTriangleByLines;
+    /**
+     * 非垂直于草图平面的三角面片被拉伸体切割算法
+     * 投影后求交，重新三角化
+     * @param extruder 拉伸体
+     * @param v1 顶点1
+     * @param v2 顶点2
+     * @param v3 顶点3
+     * @returns 切割后的三角面片
+     */
+    private static cutInVerticalTriangle;
+    /**
+     * 纠正三角面片的方向
+     * @param vtx 三角面片的顶点
+     * @param dir 基准方向: 必须是单位向量
+     */
+    private static reviseTriangleDir;
+    /**
+     * 给定三角面片的三个顶点，计算面积
+     * @param v1 顶点1
+     * @param v2 顶点2
+     * @param v3 顶点3
+     */
+    private static calcTriangleArea;
+    /**
+     * 获取归一化的三角面片法向
+     * @param v1 第一个顶点
+     * @param v2 第二个顶点
+     * @param v3 第三个顶点
+     * @returns 法向 - 单位化
+     */
+    private static getTriangleNorm;
+    /**
+     * 给定三角面片，计算其对应的uvTransform矩阵
+     * @param vertex1
+     * @param vertex2
+     * @param vertex3
+     */
+    private static getTriangleUVTransform;
+    /**
+     * 对xyz坐标进行矩阵变换
+     * @param v 坐标点
+     * @param mt 变换矩阵
+     */
+    private static transformXYZ;
+    /**
+     * 求顶点的uv值
+     * @param v 顶点坐标
+     * @param uvTransform 变换矩阵
+     */
+    private static getUVByUVTransform;
+    /**
+     * 基于分割关系求uv
+     * @param v 带求顶点
+     * @param vertex1 原三角面片顶点1
+     * @param vertex2 原三角面片顶点2
+     * @param vertex3 原三角面片顶点3
+     */
+    private static getUVBySplit;
+    private static sameDirection;
+    private static crossVector3;
+    private static V3SQLength;
+    /**
+     * 判断两个顶点是否共点
+     * @param va 顶点a
+     * @param vb 顶点b
+     * @returns
+     */
+    private static sameVertex;
+}
+export {};

package/types/algorithm/mesh/mesh_assist.d.ts

@@ -0,0 +1,21 @@
+import { Coord3 } from '../../base/coord3';
+import { types } from '../../type_define/i_types';
+export declare class MeshAssist {
+    static getEdgeNeighbourMap(faces: number[]): Map<number, number>;
+    /**
+     * 从给定的 seedF0s 出发，通过相邻关系找到满足条件的面片
+     * @param faces 原始面片集
+     * @param newFaces 找到的面片会添加到该数组中
+     * @param seedF0s 种子面片的 face 序号
+     * @param isF0Valid 面片是否满足条件
+     * @returns 返回添加的面片在原面片集中的序号（包含种子面片）
+     */
+    static pickNeighbourFaces(faces: number[], newFaces: number[], seedF0s: number[], isF0Valid: (fj0: number) => boolean, edgeNeighbourMap?: Map<number, number>): Set<number>;
+    /**
+     * 根据 nextMap，将对应点投影到 coord 平面上，生成一系列多边形
+     * @param mesh
+     * @param nextMap
+     * @param coord
+     */
+    static getLoop2ds(mesh: types.IFlatMesh, nextMap: Map<number, number[]>, coord: Coord3): types.IXY[][][];
+}

package/types/algorithm/mesh/mesh_contour.d.ts

@@ -0,0 +1,26 @@
+import { types } from '../../type_define/i_types';
+export declare class MeshContour {
+    static execute(mesh: types.IFlatMesh): types.IXY[][][];
+    static getUpMesh(mesh: types.IFlatMesh): types.IFlatMesh[];
+    /**
+     * 在投影平面上
+     * @param mesh
+     * @param point
+     */
+    static getHeightAtProjectPoint(mesh: types.IFlatMesh, point: types.IXY): number[];
+    private static _isFacesUp;
+    private static _splitMesh;
+    private static _getOuterLoop;
+    private static _idFiMap2NextMap;
+    private static _filterLoops;
+    /**
+     * 对 mesh 排序，返回的 mesh 数组按从下到上排序
+     * @param meshes
+     * @param outerLoops
+     */
+    private static _sortMeshes;
+    private static _isMeshLower;
+    private static _mixContours;
+    private static _filterSmallLoops;
+    private static _getArea;
+}