scax-engine 0.1.0

package/dist/types/surfaces/spherical-image.d.ts

@@ -0,0 +1,42 @@
+import { Vector3 } from "three";
+import Ray from "../ray/ray";
+import Surface from "./surface";
+export type SphericalImageSurfaceProps = {
+    type: "spherical-image";
+    name: string;
+    r: number;
+    position: {
+        x: number;
+        y: number;
+        z: number;
+    };
+    tilt: {
+        x: number;
+        y: number;
+    };
+    retina_extra_after: boolean;
+};
+/**
+ * 곡면 표현만 하고 굴절 효과는 없습니다.
+ * 망막 위치를 표현하는 데 사용됩니다.
+ */
+export default class SphericalImageSurface extends Surface {
+    private r;
+    private retina_extra_after;
+    private hitPoints;
+    constructor(props: SphericalImageSurfaceProps);
+    /**
+     * 반경 값이 비정상이면 평면(z = position.z)으로 처리합니다.
+     */
+    private isPlanar;
+    /**
+     * 구면 중심: 꼭지점(position)에서 반경만큼 z축 이동한 점
+     */
+    private sphereCenter;
+    getHitPoints(): Vector3[];
+    clearHitPoints(): void;
+    clearTraceHistory(): void;
+    incident(ray: Ray): Vector3 | null;
+    refract(ray: Ray): Ray | null;
+}
+//# sourceMappingURL=spherical-image.d.ts.map

package/dist/types/surfaces/spherical-surface.d.ts

@@ -0,0 +1,58 @@
+import { Vector3 } from "three";
+import { RefractiveIndexSpec } from "../optics/refractive-index";
+import Ray from "../ray/ray";
+import Surface from "./surface";
+export type SphericalSurfaceProps = {
+    type: "spherical";
+    name: string;
+    r: number;
+    position: {
+        x: number;
+        y: number;
+        z: number;
+    };
+    tilt: {
+        x: number;
+        y: number;
+    };
+    n_before: RefractiveIndexSpec;
+    n_after: RefractiveIndexSpec;
+};
+export default class SphericalSurface extends Surface {
+    private r;
+    private n_before;
+    private n_after;
+    constructor(props: SphericalSurfaceProps);
+    private refractiveIndicesForRay;
+    /**
+     * 반경이 너무 크거나 비정상 값이면 평면으로 간주합니다.
+     * (legacy 코드의 planar fallback 동작을 그대로 반영)
+     */
+    private isPlanar;
+    private worldQuaternion;
+    private worldPointToLocal;
+    private localPointToWorld;
+    private worldDirToLocal;
+    private localDirToWorld;
+    /**
+     * 구면의 중심점입니다.
+     * 이 프로젝트의 구면은 +Z 축을 기준으로 배치되며,
+     * 중심은 꼭지점(position)에서 반경만큼 Z 방향으로 이동한 위치입니다.
+     */
+    private sphereCenterLocal;
+    /**
+     * 굴절 계산에 사용할 "2번째 매질 방향" 법선을 계산합니다.
+     * - 평면: +Z 법선 사용
+     * - 구면: 반경 부호에 따라 중심-입사점 벡터 방향을 맞춰 사용
+     */
+    private normalIntoSecondMediumLocal;
+    /**
+     * 주어진 XY에서 구면의 z 위치를 계산합니다.
+     * - 반경이 매우 큰 경우(평면)에는 평면 z를 반환합니다.
+     * - 구면 정의역 밖이면 null을 반환합니다.
+     */
+    private surfaceZAtXY;
+    incident(ray: Ray): Vector3 | null;
+    refract(ray: Ray): Ray | null;
+}
+//# sourceMappingURL=spherical-surface.d.ts.map

package/dist/types/surfaces/st-surface.d.ts

@@ -0,0 +1,91 @@
+import { Vector3 } from "three";
+import { RefractiveIndexSpec } from "../optics/refractive-index";
+import Ray from "../ray/ray";
+import Surface from "./surface";
+export type STSurfaceProps = {
+    type: "compound";
+    name: string;
+    position: {
+        x: number;
+        y: number;
+        z: number;
+    };
+    tilt: {
+        x: number;
+        y: number;
+    };
+    r?: number;
+    s: number;
+    c: number;
+    ax: number;
+    n_before: RefractiveIndexSpec;
+    n: RefractiveIndexSpec;
+    n_after: RefractiveIndexSpec;
+    referencePoint?: {
+        x: number;
+        y: number;
+        z: number;
+    };
+    thickness?: number;
+};
+export default class STSurface extends Surface {
+    private static readonly POWER_EPS_D;
+    private s;
+    private c;
+    private ax;
+    private n_before;
+    private n;
+    private n_after;
+    private thickness;
+    /** 구면. vertex z = position.z + thickness (+z 쪽, 광이 나가는 쪽). */
+    private sphericalSurface;
+    /** 토릭(원통). vertex z = position.z (−z 쪽, +z 진행 시 광이 먼저 맞는 쪽). 없으면 null. */
+    private toricSurface;
+    private sphericalRadiusMm;
+    private toricRadiusPerpMm;
+    constructor(props: STSurfaceProps);
+    /**
+     * 디옵터(D)로부터 반경(mm)을 계산합니다.
+     *
+     * power(D) = (n2 - n1) / R(m)
+     *   -> R(mm) = 1000 * (n2 - n1) / power(D)
+     *
+     * 굴절력이 사실상 0이면 평면으로 간주하기 위해 +Infinity를 반환합니다.
+     */
+    private radiusFromPower;
+    private refractiveIndexAtD;
+    /**
+     * ST 구면 측: z가 더 큰 꼭지( +z 진행 시 출사 쪽 ).
+     */
+    private buildSphericalSurface;
+    /**
+     * ST 토릭 측: cylinder가 있을 때만. z가 더 작은 꼭지( +z 진행 시 입사 쪽 ).
+     */
+    private buildToricSurface;
+    private applyChromaticIndicesToSubSurfaces;
+    /**
+     * 구면/토릭 곡면의 z 교차(토릭이 구면을 관통) 방지를 위해
+     * 샘플링 영역에서 필요한 최소 중심두께를 계산합니다.
+     */
+    private optimizeThickness;
+    /**
+     * 기준점(referencePoint.z)으로부터 토릭 꼭지(z)까지의 최소 간격을 확보합니다.
+     * - 기준점과 반대 방향으로 현재 토릭 꼭지가 놓인 쪽(sign)을 유지합니다.
+     * - 최소 간격은 "토릭–구면 거리(thickness) + 안전여유"입니다.
+     */
+    private optimizeToricVertexZFromReference;
+    private samplingRadiusMm;
+    /**
+     * 구면 측 꼭지점 기준 sag(mm)
+     */
+    private sphericalSagAtXY;
+    /**
+     * 토릭 측 꼭지점 기준 sag(mm)
+     */
+    private toricSagAtXY;
+    private isOpticallyNeutral;
+    refract(ray: Ray): Ray | null;
+    incident(ray: Ray): Vector3 | null;
+    getOptimizedThicknessMm(): number;
+}
+//# sourceMappingURL=st-surface.d.ts.map

package/dist/types/surfaces/surface.d.ts

@@ -0,0 +1,32 @@
+import { Vector2, Vector3 } from "three";
+import Ray from "../ray/ray";
+export type SurfaceProps = {
+    type: string;
+    name: string;
+    position: {
+        x: number;
+        y: number;
+        z: number;
+    };
+    tilt: {
+        x: number;
+        y: number;
+    };
+};
+export default abstract class Surface {
+    protected type: string;
+    protected name: string;
+    protected incidentRays: Ray[];
+    protected refractedRays: Ray[];
+    protected position: Vector3;
+    protected tilt: Vector2;
+    protected meridians: {
+        angle: number;
+        d: number;
+    }[];
+    constructor(props: SurfaceProps);
+    abstract incident(ray: Ray): Vector3 | null;
+    abstract refract(ray: Ray): Ray | null;
+    clearTraceHistory(): void;
+}
+//# sourceMappingURL=surface.d.ts.map

package/dist/types/surfaces/toric-surface.d.ts

@@ -0,0 +1,68 @@
+import { Vector3 } from "three";
+import { RefractiveIndexSpec } from "../optics/refractive-index";
+import Ray from "../ray/ray";
+import Surface from "./surface";
+export type ToricSurfaceProps = {
+    type: "toric";
+    name: string;
+    position: {
+        x: number;
+        y: number;
+        z: number;
+    };
+    tilt: {
+        x: number;
+        y: number;
+    };
+    r_axis: number;
+    r_perp: number;
+    axis_deg?: number;
+    n_before: RefractiveIndexSpec;
+    n_after: RefractiveIndexSpec;
+};
+export default class ToricSurface extends Surface {
+    private r_axis;
+    private r_perp;
+    private n_before;
+    private n_after;
+    private axisDeg;
+    constructor(props: ToricSurfaceProps);
+    private refractiveIndicesForRay;
+    /**
+     * Toric 면의 축(meridian) 회전을 위해 사용하는 삼각함수 값입니다.
+     * 축(axis) 회전은 axisDeg(도 단위)를 사용합니다.
+     */
+    private axisTrig;
+    private worldQuaternion;
+    private worldPointToLocal;
+    private localPointToWorld;
+    private worldDirToLocal;
+    private localDirToWorld;
+    /**
+     * 월드 좌표계 (x, y)를 토릭 로컬 좌표계 (u, v)로 변환합니다.
+     * - u: 축 방향 meridian
+     * - v: 축에 수직인 meridian
+     */
+    private toLocalUV;
+    /**
+     * 로컬 좌표계에서 계산한 sag 미분(dz/du, dz/dv)을
+     * 월드 좌표계의 기울기(dz/dx, dz/dy)로 다시 매핑합니다.
+     */
+    private localDerivativesToWorld;
+    /**
+     * 반경으로부터 곡률(1/R)을 계산합니다.
+     * - 반경이 너무 크거나 무한대면 평면으로 간주하여 0 반환
+     * - 반경이 0에 너무 가까우면 비정상 값으로 NaN 반환
+     */
+    private curvature;
+    /**
+     * 주어진 XY에서 토릭면의 기하정보를 계산합니다.
+     * - sag: 꼭지점 대비 z 높이
+     * - dzdx/dzdy: 면 기울기
+     * - normal: 2번째 매질 방향을 만들 때 사용할 기본 법선
+     */
+    private geometryAtXY;
+    incident(ray: Ray): Vector3 | null;
+    refract(ray: Ray): Ray | null;
+}
+//# sourceMappingURL=toric-surface.d.ts.map

package/dist/types/utils/deg-to-tabo.d.ts

@@ -0,0 +1,3 @@
+declare const TABOToDeg: (TABOAngle: number) => number;
+export default TABOToDeg;
+//# sourceMappingURL=deg-to-tabo.d.ts.map

package/dist/types/utils/tabo-to-deg.d.ts

@@ -0,0 +1,3 @@
+declare const DegToTABO: (degree: number) => number;
+export default DegToTABO;
+//# sourceMappingURL=tabo-to-deg.d.ts.map

package/package.json

@@ -0,0 +1,64 @@
+{
+  "name": "scax-engine",
+  "version": "0.1.0",
+  "description": "Optical calculation engine for an eye model, with ESM, CJS, and UMD builds.",
+  "main": "./dist/scax-engine.cjs",
+  "module": "./dist/scax-engine.js",
+  "types": "./dist/types/index.d.ts",
+  "files": [
+    "dist",
+    "!dist/**/*.map"
+  ],
+  "scripts": {
+    "clean": "rimraf dist",
+    "build": "npm run clean && tsc -p tsconfig.types.json && rollup -c rollup.config.mjs",
+    "test": "vitest run",
+    "test:watch": "vitest",
+    "prepublishOnly": "npm run test && npm run build"
+  },
+  "repository": {
+    "type": "git",
+    "url": "git+https://github.com/vizyman/scax-engine.git"
+  },
+  "keywords": [
+    "math",
+    "calculator",
+    "typescript",
+    "library"
+  ],
+  "author": "vizyman",
+  "license": "MIT",
+  "type": "module",
+  "exports": {
+    ".": {
+      "types": "./dist/types/index.d.ts",
+      "import": "./dist/scax-engine.js",
+      "require": "./dist/scax-engine.cjs",
+      "default": "./dist/scax-engine.js"
+    }
+  },
+  "unpkg": "./dist/scax-engine.umd.js",
+  "jsdelivr": "./dist/scax-engine.umd.js",
+  "bugs": {
+    "url": "https://github.com/vizyman/scax-engine/issues"
+  },
+  "homepage": "https://github.com/vizyman/scax-engine#readme",
+  "engines": {
+    "node": ">=20"
+  },
+  "devDependencies": {
+    "@rollup/plugin-commonjs": "^29.0.2",
+    "@rollup/plugin-node-resolve": "^16.0.3",
+    "@rollup/plugin-typescript": "^12.3.0",
+    "@types/node": "^25.6.0",
+    "@types/three": "^0.184.0",
+    "rimraf": "^6.1.3",
+    "rollup": "^4.60.2",
+    "tslib": "^2.8.1",
+    "typescript": "^6.0.3",
+    "vitest": "^4.1.5"
+  },
+  "dependencies": {
+    "three": "^0.184.0"
+  }
+}