brepjs-bim 0.23.0 → 0.23.1

package/dist/brepjs-bim.cjs

@@ -206,7 +206,7 @@ function fromBrepError(inner, code, message) {
 	};
 }
 //#endregion
-//#region \0@oxc-project+runtime@0.143.0/helpers/esm/usingCtx.js
+//#region \0@oxc-project+runtime@0.146.0/helpers/esm/usingCtx.js
 function _usingCtx() {
 	var r = "function" == typeof SuppressedError ? SuppressedError : function(r, e) {
 		var n = Error();
@@ -8779,7 +8779,7 @@ function writeWallGeometry(w, spec, geomSubContextId, parentPlacementId) {
 			value: "AREA"
 		},
 		ProfileName: null,
-		Position: w.ref(writeAxis2Placement2D$1(w)),
+		Position: w.ref(writeAxis2Placement2D$1(w, [thicknessM / 2, heightM / 2])),
 		XDim: w.mkType(web_ifc.IFCPOSITIVELENGTHMEASURE, thicknessM),
 		YDim: w.mkType(web_ifc.IFCPOSITIVELENGTHMEASURE, heightM)
 	});
@@ -9000,12 +9000,12 @@ function writeRoofGeometry(w, spec, solid, geomSubContextId, parentPlacementId)
 		usedFallback: false
 	};
 }
-function writeAxis2Placement2D$1(w) {
+function writeAxis2Placement2D$1(w, location = [0, 0]) {
 	const originId = w.nextId();
 	w.writeLine({
 		expressID: originId,
 		type: web_ifc.IFCCARTESIANPOINT,
-		Coordinates: [w.mkType(web_ifc.IFCLENGTHMEASURE, 0), w.mkType(web_ifc.IFCLENGTHMEASURE, 0)]
+		Coordinates: [w.mkType(web_ifc.IFCLENGTHMEASURE, location[0]), w.mkType(web_ifc.IFCLENGTHMEASURE, location[1])]
 	});
 	const id = w.nextId();
 	w.writeLine({
@@ -11510,11 +11510,10 @@ function writeOpeningGeometry(w, guid, openingSpec, wallSpec, wallPlacementId, g
 	const offsetAlongWallM = toIfcLengthM(openingSpec.offsetAlongWall);
 	const offsetFromFloorM = toIfcLengthM(openingSpec.offsetFromFloor);
 	const thicknessM = toIfcLengthM(wallSpec.thickness);
-	const wallHeightM = toIfcLengthM(wallSpec.height);
 	const placement3DId = writeAxis2Placement3D(w, [
 		offsetAlongWallM + widthM / 2,
-		thicknessM / 2,
-		-wallHeightM / 2 + offsetFromFloorM + heightM / 2
+		thicknessM,
+		offsetFromFloorM + heightM / 2
 	], [
 		0,
 		-1,
@@ -12390,6 +12389,24 @@ var SUPPORTED_SCHEMAS = [
 * MUST call {@link SpfReader.close} (it always issues `CloseModel`, preventing
 * WASM handle leaks). The handle is `using`-friendly via `Symbol.dispose`.
 */
+var IFC_Z_UP_TRANSFORMATION = [
+	1,
+	0,
+	0,
+	0,
+	0,
+	0,
+	1,
+	0,
+	0,
+	-1,
+	0,
+	0,
+	0,
+	0,
+	0,
+	1
+];
 function deleteVector(vec) {
 	vec.delete?.();
 }
@@ -12422,6 +12439,7 @@ var SpfReader = class SpfReader {
 			if (modelId >= 0) api.CloseModel(modelId);
 			return (0, brepjs.err)(importError("OPEN_MODEL_FAILED", "web-ifc OpenModel returned an invalid model id"));
 		}
+		api.SetGeometryTransformation(modelId, IFC_Z_UP_TRANSFORMATION);
 		const schemaRaw = api.GetModelSchema(modelId);
 		if (!SUPPORTED_SCHEMAS.includes(schemaRaw)) {
 			api.CloseModel(modelId);
@@ -14121,12 +14139,14 @@ function reconstructExtrusion(reader, extrusionId, scale, worldTransform, diagno
 		diagnostics.push(issue("warning", faceResult.error.code, faceResult.error.message, extrusionId));
 		return NONE;
 	}
+	const profileFrame = readProfilePosition(reader, sweptArea.value, scale);
 	const depthMm = depth * scale * 1e3;
-	const extrudeDir = readDirection(reader, ext["ExtrudedDirection"]) ?? [
+	const sweptDir = readDirection(reader, ext["ExtrudedDirection"]) ?? [
 		0,
 		0,
 		1
 	];
+	const extrudeDir = profileFrame === null ? sweptDir : rotateByTranspose(profileFrame.linear, sweptDir);
 	const extrudeVec = [
 		extrudeDir[0] * depthMm,
 		extrudeDir[1] * depthMm,
@@ -14148,7 +14168,11 @@ function reconstructExtrusion(reader, extrusionId, scale, worldTransform, diagno
 		return NONE;
 	}
 	const localFrame = readAxis2Placement3D(reader, ext["Position"], scale);
-	const placed = placeSolid(solidResult.value, localFrame, worldTransform);
+	const placed = placeSolid(solidResult.value, [
+		worldTransform,
+		localFrame,
+		profileFrame
+	]);
 	if (!placed.ok) {
 		diagnostics.push(issue("warning", placed.error.code, placed.error.message, extrusionId));
 		return NONE;
@@ -14176,16 +14200,19 @@ function reconstructRevolution(reader, revolutionId, scale, worldTransform, diag
 		return NONE;
 	}
 	const axis1 = reader.getLine(axisRef.value);
-	const center = (axis1 === null ? void 0 : readPoint(reader, axis1["Location"], scale)) ?? [
+	const sweptCenter = (axis1 === null ? void 0 : readPoint(reader, axis1["Location"], scale)) ?? [
 		0,
 		0,
 		0
 	];
-	const direction = (axis1 === null ? void 0 : readDirection(reader, axis1["Axis"])) ?? [
+	const sweptDirection = (axis1 === null ? void 0 : readDirection(reader, axis1["Axis"])) ?? [
 		0,
 		0,
 		1
 	];
+	const profileFrame = readProfilePosition(reader, sweptArea.value, scale);
+	const center = profileFrame === null ? sweptCenter : inverseRigidPoint(profileFrame, sweptCenter);
+	const direction = profileFrame === null ? sweptDirection : rotateByTranspose(profileFrame.linear, sweptDirection);
 	const revolved = (() => {
 		try {
 			var _usingCtx3 = _usingCtx();
@@ -14212,7 +14239,11 @@ function reconstructRevolution(reader, revolutionId, scale, worldTransform, diag
 		return NONE;
 	}
 	const localFrame = readAxis2Placement3D(reader, rev["Position"], scale);
-	const placed = placeSolid(revolved.value, localFrame, worldTransform);
+	const placed = placeSolid(revolved.value, [
+		worldTransform,
+		localFrame,
+		profileFrame
+	]);
 	if (!placed.ok) {
 		diagnostics.push(issue("warning", placed.error.code, placed.error.message, revolutionId));
 		return NONE;
@@ -14225,20 +14256,9 @@ function reconstructTessellated(reader, productExpressId, diagnostics) {
 		diagnostics.push(issue("warning", "GEOMETRY_RECONSTRUCTION_FAILED", "Tessellated geometry yielded no triangles", productExpressId));
 		return NONE;
 	}
-	const stl = packBinaryStl(mesh.vertices, mesh.indices, 1e3);
-	let solid = null;
-	try {
-		const cast = (0, brepjs.castShape)((0, brepjs.getKernel)().importSTL(stl.buffer));
-		if ((0, brepjs.isSolid)(cast)) {
-			const valid = (0, brepjs.validSolid)(cast);
-			if (valid.ok) solid = valid.value;
-			else cast[Symbol.dispose]();
-		} else cast[Symbol.dispose]();
-	} catch (e) {
-		diagnostics.push(issue("info", "TESSELLATION_NOT_MANIFOLD", `STL round-trip failed: ${errMsg$1(e)}`, productExpressId));
-	}
+	const solid = sewMeshToSolid(mesh, 1e3, productExpressId, diagnostics);
 	if (solid !== null) {
-		diagnostics.push(issue("info", "TESSELLATED_MANIFOLD", "Tessellated mesh recovered as a closed solid via STL round-trip", productExpressId));
+		diagnostics.push(issue("info", "TESSELLATED_MANIFOLD", "Tessellated mesh recovered as a closed solid by sewing its triangles", productExpressId));
 		return {
 			kind: "SOLID",
 			solid,
@@ -14253,6 +14273,52 @@ function reconstructTessellated(reader, productExpressId, diagnostics) {
 		diagnostic: "TESSELLATED_LOSSY"
 	};
 }
+/** Loose enough to weld float32 vertices web-ifc emits per face, in mm. */
+var MESH_SEW_TOLERANCE_MM = .001;
+/**
+* Brands a fresh kernel result and releases its pre-downcast arena slot when
+* the cast moved to a new one (occt-wasm); in-place kernels share the slot and
+* are left alone. Mirrors brepjs's internal castResultShape.
+*/
+function castFreshResult(raw) {
+	const cast = (0, brepjs.castShape)(raw);
+	const rawId = raw.id;
+	const castId = cast.wrapped.id;
+	if (!(rawId !== void 0 && castId !== void 0 ? rawId === castId : cast.wrapped === raw)) (0, brepjs.getKernel)().dispose(raw);
+	return cast;
+}
+/**
+* Sews the streamed triangles into a closed solid through the kernel's mesh
+* builder, the same path brepjs's mesh importers use. Null when the mesh is
+* open or the kernel rejects it.
+*/
+function sewMeshToSolid(mesh, scaleToMm, productExpressId, diagnostics) {
+	const kernel = (0, brepjs.getKernel)();
+	const { vertices, indices } = mesh;
+	const point = (index) => [
+		(vertices[index * 3] ?? 0) * scaleToMm,
+		(vertices[index * 3 + 1] ?? 0) * scaleToMm,
+		(vertices[index * 3 + 2] ?? 0) * scaleToMm
+	];
+	const triangles = [];
+	for (let t = 0; t + 2 < indices.length; t += 3) {
+		const face = kernel.buildTriFace(point(indices[t] ?? 0), point(indices[t + 1] ?? 0), point(indices[t + 2] ?? 0));
+		if (face !== null) triangles.push(face);
+	}
+	if (triangles.length === 0) return null;
+	try {
+		const cast = castFreshResult(kernel.sewAndSolidify(triangles, MESH_SEW_TOLERANCE_MM));
+		const valid = (0, brepjs.isSolid)(cast) ? (0, brepjs.validSolid)(cast) : null;
+		if (valid !== null && valid.ok) return valid.value;
+		cast[Symbol.dispose]();
+		return null;
+	} catch (e) {
+		diagnostics.push(issue("info", "TESSELLATION_NOT_MANIFOLD", `Mesh sewing failed: ${errMsg$1(e)}`, productExpressId));
+		return null;
+	} finally {
+		for (const face of triangles) kernel.dispose(face);
+	}
+}
 function collectMesh(reader, productExpressId) {
 	const vertices = [];
 	const indices = [];
@@ -14281,7 +14347,8 @@ function collectMesh(reader, productExpressId) {
 				}
 			}
 		} finally {
-			flatMesh.delete();
+			releaseEmbind(flatMesh.geometries);
+			releaseEmbind(flatMesh);
 		}
 	});
 	if (vertices.length === 0) return null;
@@ -14290,25 +14357,6 @@ function collectMesh(reader, productExpressId) {
 		indices: new Uint32Array(indices)
 	};
 }
-function packBinaryStl(vertices, indices, scaleToMm) {
-	const triCount = Math.floor(indices.length / 3);
-	const buffer = /* @__PURE__ */ new ArrayBuffer(84 + triCount * 50);
-	const view = new DataView(buffer);
-	view.setUint32(80, triCount, true);
-	let offset = 84;
-	for (let t = 0; t < triCount; t++) {
-		offset += 12;
-		for (let c = 0; c < 3; c++) {
-			const vi = (indices[t * 3 + c] ?? 0) * 3;
-			view.setFloat32(offset, (vertices[vi] ?? 0) * scaleToMm, true);
-			view.setFloat32(offset + 4, (vertices[vi + 1] ?? 0) * scaleToMm, true);
-			view.setFloat32(offset + 8, (vertices[vi + 2] ?? 0) * scaleToMm, true);
-			offset += 12;
-		}
-		offset += 2;
-	}
-	return new Uint8Array(buffer);
-}
 var PARAMETRIC_PROFILE_BUILDERS = /* @__PURE__ */ new Map([
 	[web_ifc.IFCRECTANGLEPROFILEDEF, (f) => ({
 		kind: "RECTANGULAR",
@@ -14452,19 +14500,99 @@ function profileToFace(profile) {
 	if (!face.ok) return (0, brepjs.err)(importError("GEOMETRY_RECONSTRUCTION_FAILED", `Profile face build failed: ${face.error.message}`));
 	return (0, brepjs.ok)(face.value);
 }
-function placeSolid(solid, localFrame, worldTransform) {
-	let current = solid;
-	for (const transform of [localFrame, worldTransform]) {
-		if (transform === null || isIdentity(transform)) continue;
-		const applied = (0, brepjs.applyMatrix)(current, transform);
-		if (!applied.ok) {
-			current[Symbol.dispose]();
-			return (0, brepjs.err)(importError("PLACEMENT_READ_FAILED", `Placement transform failed: ${applied.error.message}`));
-		}
-		current[Symbol.dispose]();
-		current = applied.value;
-	}
-	return (0, brepjs.ok)(current);
+/**
+* IfcParameterizedProfileDef.Position places the profile inside the swept
+* solid's XY plane; the writers rely on it to corner-anchor centred rectangle
+* profiles. Absent on arbitrary (polyline) profiles.
+*/
+function readProfilePosition(reader, profileExpressId, scale) {
+	const def = reader.getLine(profileExpressId);
+	const positionRef = def === null ? void 0 : asRef(def["Position"]);
+	if (positionRef === void 0) return null;
+	const position = reader.getLine(positionRef.value);
+	if (position === null) return null;
+	const [x, y] = readPoint2D(reader, position["Location"], scale) ?? [0, 0];
+	const [dx, dy] = readDirection2D(reader, position["RefDirection"]) ?? [1, 0];
+	const len = Math.hypot(dx, dy);
+	const c = len < 1e-12 ? 1 : dx / len;
+	const s = len < 1e-12 ? 0 : dy / len;
+	return {
+		linear: [
+			c,
+			-s,
+			0,
+			s,
+			c,
+			0,
+			0,
+			0,
+			1
+		],
+		translation: [
+			x,
+			y,
+			0
+		]
+	};
+}
+/**
+* Applies the placement chain (outermost first) as one composed rigid motion.
+* A second applyMatrix on an already-transformed solid can fail BRepCheck on
+* occt-wasm, so the frames are never applied one after another.
+*/
+function placeSolid(solid, frames) {
+	let composed = null;
+	for (const frame of frames) {
+		if (frame === null) continue;
+		composed = composed === null ? frame : composeRigid(composed, frame);
+	}
+	if (composed === null || isIdentity(composed)) return (0, brepjs.ok)(solid);
+	const applied = (0, brepjs.applyMatrix)(solid, composed);
+	solid[Symbol.dispose]();
+	if (!applied.ok) return (0, brepjs.err)(importError("PLACEMENT_READ_FAILED", `Placement transform failed: ${applied.error.message}`));
+	return (0, brepjs.ok)(applied.value);
+}
+/** `outer ∘ inner`: the motion that applies `inner` first, then `outer`. */
+function composeRigid(outer, inner) {
+	const [a00, a01, a02, a10, a11, a12, a20, a21, a22] = outer.linear;
+	const [b00, b01, b02, b10, b11, b12, b20, b21, b22] = inner.linear;
+	const [tx, ty, tz] = inner.translation;
+	const [ux, uy, uz] = outer.translation;
+	return {
+		linear: [
+			a00 * b00 + a01 * b10 + a02 * b20,
+			a00 * b01 + a01 * b11 + a02 * b21,
+			a00 * b02 + a01 * b12 + a02 * b22,
+			a10 * b00 + a11 * b10 + a12 * b20,
+			a10 * b01 + a11 * b11 + a12 * b21,
+			a10 * b02 + a11 * b12 + a12 * b22,
+			a20 * b00 + a21 * b10 + a22 * b20,
+			a20 * b01 + a21 * b11 + a22 * b21,
+			a20 * b02 + a21 * b12 + a22 * b22
+		],
+		translation: [
+			a00 * tx + a01 * ty + a02 * tz + ux,
+			a10 * tx + a11 * ty + a12 * tz + uy,
+			a20 * tx + a21 * ty + a22 * tz + uz
+		]
+	};
+}
+/** Rᵀ·v for a rigid frame's rotation part. */
+function rotateByTranspose(linear, v) {
+	return [
+		linear[0] * v[0] + linear[3] * v[1] + linear[6] * v[2],
+		linear[1] * v[0] + linear[4] * v[1] + linear[7] * v[2],
+		linear[2] * v[0] + linear[5] * v[1] + linear[8] * v[2]
+	];
+}
+/** Maps a point from the frame's parent space into the frame: Rᵀ·(p − t). */
+function inverseRigidPoint(frame, p) {
+	const t = frame.translation;
+	return rotateByTranspose(frame.linear, [
+		p[0] - t[0],
+		p[1] - t[1],
+		p[2] - t[2]
+	]);
 }
 function finalizeSolid(solid, entity, diagnostics) {
 	if (!(0, brepjs.isSolid)(solid)) {
@@ -14613,6 +14741,24 @@ function readPolylinePoints(reader, ref, scale) {
 	}
 	return out.length >= 3 ? out : void 0;
 }
+function readPoint2D(reader, ref, scale) {
+	const pointRef = asRef(ref);
+	if (pointRef === void 0) return void 0;
+	const pt = reader.getLine(pointRef.value);
+	if (pt === null) return void 0;
+	const coords = asMeasureArray(pt["Coordinates"]);
+	if (coords.length < 2) return void 0;
+	return [(coords[0] ?? 0) * scale * 1e3, (coords[1] ?? 0) * scale * 1e3];
+}
+function readDirection2D(reader, ref) {
+	const dirRef = asRef(ref);
+	if (dirRef === void 0) return void 0;
+	const dir = reader.getLine(dirRef.value);
+	if (dir === null) return void 0;
+	const ratios = asMeasureArray(dir["DirectionRatios"]);
+	if (ratios.length < 2) return void 0;
+	return [ratios[0] ?? 0, ratios[1] ?? 0];
+}
 function readPoint(reader, ref, scale) {
 	const pointRef = asRef(ref);
 	if (pointRef === void 0) return void 0;
@@ -14639,6 +14785,9 @@ function readDirection(reader, ref) {
 		ratios[2] ?? 0
 	];
 }
+function releaseEmbind(value) {
+	value.delete?.();
+}
 function asRef(value) {
 	if (value !== null && typeof value === "object" && "value" in value) {
 		const v = value.value;
@@ -28855,13 +29004,25 @@ function specLocalShift(el) {
 		0
 	];
 }
+/** Skip the resolver occurrence wrappers already owned by `cumulativeFrame`,
+*  then recover only the Body's outer literal Datum translations. An unexpected
+*  wrapper shape returns zero rather than risking an occurrence transform twice. */
+function bodyLocalTranslation(el, occurrenceTransformDepth) {
+	let body = el.geometry;
+	for (let i = 0; i < occurrenceTransformDepth; i++) {
+		if (body.kind !== "Translate" && body.kind !== "Rotate") return ZERO_TRANSLATION;
+		body = body.target;
+	}
+	return peelTranslates(body).total;
+}
 /** World frame of an element's own body: the walk's cumulative frame (all
 *  authored transforms, ancestors + own) composed with the body's authored axes
-*  (`axisX`/`axisZ` props) and its local origin (`origin` prop + centring
-*  shift). Under no rotation this reduces to `composedOrigin` + shift. */
-function elementBodyFrame(el, cumulativeFrame) {
+*  (`axisX`/`axisZ` props) and its local origin (`origin` prop + literal Body
+*  Datum + centring shift). Under no rotation this reduces to `composedOrigin`
+*  + shift. */
+function elementBodyFrame(el, cumulativeFrame, occurrenceTransformDepth) {
 	return frameMul(frameMul(cumulativeFrame, frameFromPlacement({
-		origin: authoredSpecOrigin(el),
+		origin: addTranslation(authoredSpecOrigin(el), bodyLocalTranslation(el, occurrenceTransformDepth)),
 		axisX: authoredAxisX(el),
 		axisZ: authoredAxisZ(el)
 	})), translationFrame(specLocalShift(el)));
@@ -28880,10 +29041,10 @@ function civilNodeFrame(el, cumulativeFrame) {
 *  `axisZ` are recomputed from the element's world frame relative to its spatial
 *  container. Stair/ramp flights, which carry their own per-flight frame, are
 *  each re-placed the same way. */
-function rotatedRoutedInput(flatInput, el, cumulativeFrame, spatialFrame) {
+function rotatedRoutedInput(flatInput, el, cumulativeFrame, occurrenceTransformDepth, spatialFrame) {
 	const toSpatial = frameInverse(spatialFrame);
 	if (Array.isArray(flatInput["flights"])) {
-		const elementFrame = elementBodyFrame(el, cumulativeFrame);
+		const elementFrame = elementBodyFrame(el, cumulativeFrame, occurrenceTransformDepth);
 		const flights = flatInput["flights"];
 		return {
 			...flatInput,
@@ -28924,7 +29085,7 @@ function rotatedRoutedInput(flatInput, el, cumulativeFrame, spatialFrame) {
 			})
 		};
 	}
-	const placed = decomposeFrame(frameMul(toSpatial, elementBodyFrame(el, cumulativeFrame)));
+	const placed = decomposeFrame(frameMul(toSpatial, elementBodyFrame(el, cumulativeFrame, occurrenceTransformDepth)));
 	return {
 		...flatInput,
 		origin: placed.origin,
@@ -29089,6 +29250,7 @@ function familiesToBim(root, options) {
 	const walk = (el, state) => {
 		const rotatedHere = state.rotated || hasRotateOp(el);
 		const cumulativeTranslationHere = addTranslation(state.cumulativeTranslation, authoredTranslation(el));
+		const occurrenceTransformDepthHere = state.occurrenceTransformDepth + el.localTransforms.length;
 		const cumulativeFrameHere = frameMul(state.cumulativeFrame, frameFromOps(el.localTransforms));
 		let proxiedHere = false;
 		let nextRotated = rotatedHere;
@@ -29152,7 +29314,7 @@ function familiesToBim(root, options) {
 			}
 			const routedInput = el.semantics?.kind === "product" ? civilProductSpecInput(el) : ("input" in route ? route.input : specInput)(el);
 			const rotatedBase = Array.isArray(routedInput["flights"]) ? specInput(el) : routedInput;
-			const placedInput = rotatedHere ? rotatedRoutedInput(rotatedBase, el, cumulativeFrameHere, state.spatialFrame) : usesAuthoredCivilHierarchy ? relativeSpecInput(routedInput, state.projectedSpatialTranslation) : routedInput;
+			const placedInput = rotatedHere ? rotatedRoutedInput(rotatedBase, el, cumulativeFrameHere, occurrenceTransformDepthHere, state.spatialFrame) : usesAuthoredCivilHierarchy ? relativeSpecInput(routedInput, state.projectedSpatialTranslation) : routedInput;
 			const parsed = route.parse(placedInput);
 			if (!parsed.ok) return parsed;
 			const added = route.add(model, parsed.value, el.keyPath);
@@ -29189,6 +29351,7 @@ function familiesToBim(root, options) {
 				rotated: nextRotated,
 				cumulativeTranslation: cumulativeTranslationHere,
 				projectedSpatialTranslation: nextProjectedSpatialTranslation,
+				occurrenceTransformDepth: occurrenceTransformDepthHere,
 				cumulativeFrame: nextCumulativeFrame,
 				spatialFrame: nextSpatialFrame
 			});
@@ -29202,6 +29365,7 @@ function familiesToBim(root, options) {
 		rotated: false,
 		cumulativeTranslation: ZERO_TRANSLATION,
 		projectedSpatialTranslation: ZERO_TRANSLATION,
+		occurrenceTransformDepth: 0,
 		cumulativeFrame: IDENTITY_FRAME,
 		spatialFrame: IDENTITY_FRAME
 	});

package/dist/brepjs-bim.js

@@ -183,7 +183,7 @@ function fromBrepError(inner, code, message) {
 	};
 }
 //#endregion
-//#region \0@oxc-project+runtime@0.143.0/helpers/esm/usingCtx.js
+//#region \0@oxc-project+runtime@0.146.0/helpers/esm/usingCtx.js
 function _usingCtx() {
 	var r = "function" == typeof SuppressedError ? SuppressedError : function(r, e) {
 		var n = Error();
@@ -8756,7 +8756,7 @@ function writeWallGeometry(w, spec, geomSubContextId, parentPlacementId) {
 			value: "AREA"
 		},
 		ProfileName: null,
-		Position: w.ref(writeAxis2Placement2D$1(w)),
+		Position: w.ref(writeAxis2Placement2D$1(w, [thicknessM / 2, heightM / 2])),
 		XDim: w.mkType(WebIFC.IFCPOSITIVELENGTHMEASURE, thicknessM),
 		YDim: w.mkType(WebIFC.IFCPOSITIVELENGTHMEASURE, heightM)
 	});
@@ -8977,12 +8977,12 @@ function writeRoofGeometry(w, spec, solid, geomSubContextId, parentPlacementId)
 		usedFallback: false
 	};
 }
-function writeAxis2Placement2D$1(w) {
+function writeAxis2Placement2D$1(w, location = [0, 0]) {
 	const originId = w.nextId();
 	w.writeLine({
 		expressID: originId,
 		type: WebIFC.IFCCARTESIANPOINT,
-		Coordinates: [w.mkType(WebIFC.IFCLENGTHMEASURE, 0), w.mkType(WebIFC.IFCLENGTHMEASURE, 0)]
+		Coordinates: [w.mkType(WebIFC.IFCLENGTHMEASURE, location[0]), w.mkType(WebIFC.IFCLENGTHMEASURE, location[1])]
 	});
 	const id = w.nextId();
 	w.writeLine({
@@ -11487,11 +11487,10 @@ function writeOpeningGeometry(w, guid, openingSpec, wallSpec, wallPlacementId, g
 	const offsetAlongWallM = toIfcLengthM(openingSpec.offsetAlongWall);
 	const offsetFromFloorM = toIfcLengthM(openingSpec.offsetFromFloor);
 	const thicknessM = toIfcLengthM(wallSpec.thickness);
-	const wallHeightM = toIfcLengthM(wallSpec.height);
 	const placement3DId = writeAxis2Placement3D(w, [
 		offsetAlongWallM + widthM / 2,
-		thicknessM / 2,
-		-wallHeightM / 2 + offsetFromFloorM + heightM / 2
+		thicknessM,
+		offsetFromFloorM + heightM / 2
 	], [
 		0,
 		-1,
@@ -12367,6 +12366,24 @@ var SUPPORTED_SCHEMAS = [
 * MUST call {@link SpfReader.close} (it always issues `CloseModel`, preventing
 * WASM handle leaks). The handle is `using`-friendly via `Symbol.dispose`.
 */
+var IFC_Z_UP_TRANSFORMATION = [
+	1,
+	0,
+	0,
+	0,
+	0,
+	0,
+	1,
+	0,
+	0,
+	-1,
+	0,
+	0,
+	0,
+	0,
+	0,
+	1
+];
 function deleteVector(vec) {
 	vec.delete?.();
 }
@@ -12399,6 +12416,7 @@ var SpfReader = class SpfReader {
 			if (modelId >= 0) api.CloseModel(modelId);
 			return err(importError("OPEN_MODEL_FAILED", "web-ifc OpenModel returned an invalid model id"));
 		}
+		api.SetGeometryTransformation(modelId, IFC_Z_UP_TRANSFORMATION);
 		const schemaRaw = api.GetModelSchema(modelId);
 		if (!SUPPORTED_SCHEMAS.includes(schemaRaw)) {
 			api.CloseModel(modelId);
@@ -14098,12 +14116,14 @@ function reconstructExtrusion(reader, extrusionId, scale, worldTransform, diagno
 		diagnostics.push(issue("warning", faceResult.error.code, faceResult.error.message, extrusionId));
 		return NONE;
 	}
+	const profileFrame = readProfilePosition(reader, sweptArea.value, scale);
 	const depthMm = depth * scale * 1e3;
-	const extrudeDir = readDirection(reader, ext["ExtrudedDirection"]) ?? [
+	const sweptDir = readDirection(reader, ext["ExtrudedDirection"]) ?? [
 		0,
 		0,
 		1
 	];
+	const extrudeDir = profileFrame === null ? sweptDir : rotateByTranspose(profileFrame.linear, sweptDir);
 	const extrudeVec = [
 		extrudeDir[0] * depthMm,
 		extrudeDir[1] * depthMm,
@@ -14125,7 +14145,11 @@ function reconstructExtrusion(reader, extrusionId, scale, worldTransform, diagno
 		return NONE;
 	}
 	const localFrame = readAxis2Placement3D(reader, ext["Position"], scale);
-	const placed = placeSolid(solidResult.value, localFrame, worldTransform);
+	const placed = placeSolid(solidResult.value, [
+		worldTransform,
+		localFrame,
+		profileFrame
+	]);
 	if (!placed.ok) {
 		diagnostics.push(issue("warning", placed.error.code, placed.error.message, extrusionId));
 		return NONE;
@@ -14153,16 +14177,19 @@ function reconstructRevolution(reader, revolutionId, scale, worldTransform, diag
 		return NONE;
 	}
 	const axis1 = reader.getLine(axisRef.value);
-	const center = (axis1 === null ? void 0 : readPoint(reader, axis1["Location"], scale)) ?? [
+	const sweptCenter = (axis1 === null ? void 0 : readPoint(reader, axis1["Location"], scale)) ?? [
 		0,
 		0,
 		0
 	];
-	const direction = (axis1 === null ? void 0 : readDirection(reader, axis1["Axis"])) ?? [
+	const sweptDirection = (axis1 === null ? void 0 : readDirection(reader, axis1["Axis"])) ?? [
 		0,
 		0,
 		1
 	];
+	const profileFrame = readProfilePosition(reader, sweptArea.value, scale);
+	const center = profileFrame === null ? sweptCenter : inverseRigidPoint(profileFrame, sweptCenter);
+	const direction = profileFrame === null ? sweptDirection : rotateByTranspose(profileFrame.linear, sweptDirection);
 	const revolved = (() => {
 		try {
 			var _usingCtx3 = _usingCtx();
@@ -14189,7 +14216,11 @@ function reconstructRevolution(reader, revolutionId, scale, worldTransform, diag
 		return NONE;
 	}
 	const localFrame = readAxis2Placement3D(reader, rev["Position"], scale);
-	const placed = placeSolid(revolved.value, localFrame, worldTransform);
+	const placed = placeSolid(revolved.value, [
+		worldTransform,
+		localFrame,
+		profileFrame
+	]);
 	if (!placed.ok) {
 		diagnostics.push(issue("warning", placed.error.code, placed.error.message, revolutionId));
 		return NONE;
@@ -14202,20 +14233,9 @@ function reconstructTessellated(reader, productExpressId, diagnostics) {
 		diagnostics.push(issue("warning", "GEOMETRY_RECONSTRUCTION_FAILED", "Tessellated geometry yielded no triangles", productExpressId));
 		return NONE;
 	}
-	const stl = packBinaryStl(mesh.vertices, mesh.indices, 1e3);
-	let solid = null;
-	try {
-		const cast = castShape(getKernel().importSTL(stl.buffer));
-		if (isSolid(cast)) {
-			const valid = validSolid(cast);
-			if (valid.ok) solid = valid.value;
-			else cast[Symbol.dispose]();
-		} else cast[Symbol.dispose]();
-	} catch (e) {
-		diagnostics.push(issue("info", "TESSELLATION_NOT_MANIFOLD", `STL round-trip failed: ${errMsg$1(e)}`, productExpressId));
-	}
+	const solid = sewMeshToSolid(mesh, 1e3, productExpressId, diagnostics);
 	if (solid !== null) {
-		diagnostics.push(issue("info", "TESSELLATED_MANIFOLD", "Tessellated mesh recovered as a closed solid via STL round-trip", productExpressId));
+		diagnostics.push(issue("info", "TESSELLATED_MANIFOLD", "Tessellated mesh recovered as a closed solid by sewing its triangles", productExpressId));
 		return {
 			kind: "SOLID",
 			solid,
@@ -14230,6 +14250,52 @@ function reconstructTessellated(reader, productExpressId, diagnostics) {
 		diagnostic: "TESSELLATED_LOSSY"
 	};
 }
+/** Loose enough to weld float32 vertices web-ifc emits per face, in mm. */
+var MESH_SEW_TOLERANCE_MM = .001;
+/**
+* Brands a fresh kernel result and releases its pre-downcast arena slot when
+* the cast moved to a new one (occt-wasm); in-place kernels share the slot and
+* are left alone. Mirrors brepjs's internal castResultShape.
+*/
+function castFreshResult(raw) {
+	const cast = castShape(raw);
+	const rawId = raw.id;
+	const castId = cast.wrapped.id;
+	if (!(rawId !== void 0 && castId !== void 0 ? rawId === castId : cast.wrapped === raw)) getKernel().dispose(raw);
+	return cast;
+}
+/**
+* Sews the streamed triangles into a closed solid through the kernel's mesh
+* builder, the same path brepjs's mesh importers use. Null when the mesh is
+* open or the kernel rejects it.
+*/
+function sewMeshToSolid(mesh, scaleToMm, productExpressId, diagnostics) {
+	const kernel = getKernel();
+	const { vertices, indices } = mesh;
+	const point = (index) => [
+		(vertices[index * 3] ?? 0) * scaleToMm,
+		(vertices[index * 3 + 1] ?? 0) * scaleToMm,
+		(vertices[index * 3 + 2] ?? 0) * scaleToMm
+	];
+	const triangles = [];
+	for (let t = 0; t + 2 < indices.length; t += 3) {
+		const face = kernel.buildTriFace(point(indices[t] ?? 0), point(indices[t + 1] ?? 0), point(indices[t + 2] ?? 0));
+		if (face !== null) triangles.push(face);
+	}
+	if (triangles.length === 0) return null;
+	try {
+		const cast = castFreshResult(kernel.sewAndSolidify(triangles, MESH_SEW_TOLERANCE_MM));
+		const valid = isSolid(cast) ? validSolid(cast) : null;
+		if (valid !== null && valid.ok) return valid.value;
+		cast[Symbol.dispose]();
+		return null;
+	} catch (e) {
+		diagnostics.push(issue("info", "TESSELLATION_NOT_MANIFOLD", `Mesh sewing failed: ${errMsg$1(e)}`, productExpressId));
+		return null;
+	} finally {
+		for (const face of triangles) kernel.dispose(face);
+	}
+}
 function collectMesh(reader, productExpressId) {
 	const vertices = [];
 	const indices = [];
@@ -14258,7 +14324,8 @@ function collectMesh(reader, productExpressId) {
 				}
 			}
 		} finally {
-			flatMesh.delete();
+			releaseEmbind(flatMesh.geometries);
+			releaseEmbind(flatMesh);
 		}
 	});
 	if (vertices.length === 0) return null;
@@ -14267,25 +14334,6 @@ function collectMesh(reader, productExpressId) {
 		indices: new Uint32Array(indices)
 	};
 }
-function packBinaryStl(vertices, indices, scaleToMm) {
-	const triCount = Math.floor(indices.length / 3);
-	const buffer = /* @__PURE__ */ new ArrayBuffer(84 + triCount * 50);
-	const view = new DataView(buffer);
-	view.setUint32(80, triCount, true);
-	let offset = 84;
-	for (let t = 0; t < triCount; t++) {
-		offset += 12;
-		for (let c = 0; c < 3; c++) {
-			const vi = (indices[t * 3 + c] ?? 0) * 3;
-			view.setFloat32(offset, (vertices[vi] ?? 0) * scaleToMm, true);
-			view.setFloat32(offset + 4, (vertices[vi + 1] ?? 0) * scaleToMm, true);
-			view.setFloat32(offset + 8, (vertices[vi + 2] ?? 0) * scaleToMm, true);
-			offset += 12;
-		}
-		offset += 2;
-	}
-	return new Uint8Array(buffer);
-}
 var PARAMETRIC_PROFILE_BUILDERS = /* @__PURE__ */ new Map([
 	[WebIFC.IFCRECTANGLEPROFILEDEF, (f) => ({
 		kind: "RECTANGULAR",
@@ -14429,19 +14477,99 @@ function profileToFace(profile) {
 	if (!face.ok) return err(importError("GEOMETRY_RECONSTRUCTION_FAILED", `Profile face build failed: ${face.error.message}`));
 	return ok(face.value);
 }
-function placeSolid(solid, localFrame, worldTransform) {
-	let current = solid;
-	for (const transform of [localFrame, worldTransform]) {
-		if (transform === null || isIdentity(transform)) continue;
-		const applied = applyMatrix(current, transform);
-		if (!applied.ok) {
-			current[Symbol.dispose]();
-			return err(importError("PLACEMENT_READ_FAILED", `Placement transform failed: ${applied.error.message}`));
-		}
-		current[Symbol.dispose]();
-		current = applied.value;
-	}
-	return ok(current);
+/**
+* IfcParameterizedProfileDef.Position places the profile inside the swept
+* solid's XY plane; the writers rely on it to corner-anchor centred rectangle
+* profiles. Absent on arbitrary (polyline) profiles.
+*/
+function readProfilePosition(reader, profileExpressId, scale) {
+	const def = reader.getLine(profileExpressId);
+	const positionRef = def === null ? void 0 : asRef(def["Position"]);
+	if (positionRef === void 0) return null;
+	const position = reader.getLine(positionRef.value);
+	if (position === null) return null;
+	const [x, y] = readPoint2D(reader, position["Location"], scale) ?? [0, 0];
+	const [dx, dy] = readDirection2D(reader, position["RefDirection"]) ?? [1, 0];
+	const len = Math.hypot(dx, dy);
+	const c = len < 1e-12 ? 1 : dx / len;
+	const s = len < 1e-12 ? 0 : dy / len;
+	return {
+		linear: [
+			c,
+			-s,
+			0,
+			s,
+			c,
+			0,
+			0,
+			0,
+			1
+		],
+		translation: [
+			x,
+			y,
+			0
+		]
+	};
+}
+/**
+* Applies the placement chain (outermost first) as one composed rigid motion.
+* A second applyMatrix on an already-transformed solid can fail BRepCheck on
+* occt-wasm, so the frames are never applied one after another.
+*/
+function placeSolid(solid, frames) {
+	let composed = null;
+	for (const frame of frames) {
+		if (frame === null) continue;
+		composed = composed === null ? frame : composeRigid(composed, frame);
+	}
+	if (composed === null || isIdentity(composed)) return ok(solid);
+	const applied = applyMatrix(solid, composed);
+	solid[Symbol.dispose]();
+	if (!applied.ok) return err(importError("PLACEMENT_READ_FAILED", `Placement transform failed: ${applied.error.message}`));
+	return ok(applied.value);
+}
+/** `outer ∘ inner`: the motion that applies `inner` first, then `outer`. */
+function composeRigid(outer, inner) {
+	const [a00, a01, a02, a10, a11, a12, a20, a21, a22] = outer.linear;
+	const [b00, b01, b02, b10, b11, b12, b20, b21, b22] = inner.linear;
+	const [tx, ty, tz] = inner.translation;
+	const [ux, uy, uz] = outer.translation;
+	return {
+		linear: [
+			a00 * b00 + a01 * b10 + a02 * b20,
+			a00 * b01 + a01 * b11 + a02 * b21,
+			a00 * b02 + a01 * b12 + a02 * b22,
+			a10 * b00 + a11 * b10 + a12 * b20,
+			a10 * b01 + a11 * b11 + a12 * b21,
+			a10 * b02 + a11 * b12 + a12 * b22,
+			a20 * b00 + a21 * b10 + a22 * b20,
+			a20 * b01 + a21 * b11 + a22 * b21,
+			a20 * b02 + a21 * b12 + a22 * b22
+		],
+		translation: [
+			a00 * tx + a01 * ty + a02 * tz + ux,
+			a10 * tx + a11 * ty + a12 * tz + uy,
+			a20 * tx + a21 * ty + a22 * tz + uz
+		]
+	};
+}
+/** Rᵀ·v for a rigid frame's rotation part. */
+function rotateByTranspose(linear, v) {
+	return [
+		linear[0] * v[0] + linear[3] * v[1] + linear[6] * v[2],
+		linear[1] * v[0] + linear[4] * v[1] + linear[7] * v[2],
+		linear[2] * v[0] + linear[5] * v[1] + linear[8] * v[2]
+	];
+}
+/** Maps a point from the frame's parent space into the frame: Rᵀ·(p − t). */
+function inverseRigidPoint(frame, p) {
+	const t = frame.translation;
+	return rotateByTranspose(frame.linear, [
+		p[0] - t[0],
+		p[1] - t[1],
+		p[2] - t[2]
+	]);
 }
 function finalizeSolid(solid, entity, diagnostics) {
 	if (!isSolid(solid)) {
@@ -14590,6 +14718,24 @@ function readPolylinePoints(reader, ref, scale) {
 	}
 	return out.length >= 3 ? out : void 0;
 }
+function readPoint2D(reader, ref, scale) {
+	const pointRef = asRef(ref);
+	if (pointRef === void 0) return void 0;
+	const pt = reader.getLine(pointRef.value);
+	if (pt === null) return void 0;
+	const coords = asMeasureArray(pt["Coordinates"]);
+	if (coords.length < 2) return void 0;
+	return [(coords[0] ?? 0) * scale * 1e3, (coords[1] ?? 0) * scale * 1e3];
+}
+function readDirection2D(reader, ref) {
+	const dirRef = asRef(ref);
+	if (dirRef === void 0) return void 0;
+	const dir = reader.getLine(dirRef.value);
+	if (dir === null) return void 0;
+	const ratios = asMeasureArray(dir["DirectionRatios"]);
+	if (ratios.length < 2) return void 0;
+	return [ratios[0] ?? 0, ratios[1] ?? 0];
+}
 function readPoint(reader, ref, scale) {
 	const pointRef = asRef(ref);
 	if (pointRef === void 0) return void 0;
@@ -14616,6 +14762,9 @@ function readDirection(reader, ref) {
 		ratios[2] ?? 0
 	];
 }
+function releaseEmbind(value) {
+	value.delete?.();
+}
 function asRef(value) {
 	if (value !== null && typeof value === "object" && "value" in value) {
 		const v = value.value;
@@ -28832,13 +28981,25 @@ function specLocalShift(el) {
 		0
 	];
 }
+/** Skip the resolver occurrence wrappers already owned by `cumulativeFrame`,
+*  then recover only the Body's outer literal Datum translations. An unexpected
+*  wrapper shape returns zero rather than risking an occurrence transform twice. */
+function bodyLocalTranslation(el, occurrenceTransformDepth) {
+	let body = el.geometry;
+	for (let i = 0; i < occurrenceTransformDepth; i++) {
+		if (body.kind !== "Translate" && body.kind !== "Rotate") return ZERO_TRANSLATION;
+		body = body.target;
+	}
+	return peelTranslates(body).total;
+}
 /** World frame of an element's own body: the walk's cumulative frame (all
 *  authored transforms, ancestors + own) composed with the body's authored axes
-*  (`axisX`/`axisZ` props) and its local origin (`origin` prop + centring
-*  shift). Under no rotation this reduces to `composedOrigin` + shift. */
-function elementBodyFrame(el, cumulativeFrame) {
+*  (`axisX`/`axisZ` props) and its local origin (`origin` prop + literal Body
+*  Datum + centring shift). Under no rotation this reduces to `composedOrigin`
+*  + shift. */
+function elementBodyFrame(el, cumulativeFrame, occurrenceTransformDepth) {
 	return frameMul(frameMul(cumulativeFrame, frameFromPlacement({
-		origin: authoredSpecOrigin(el),
+		origin: addTranslation(authoredSpecOrigin(el), bodyLocalTranslation(el, occurrenceTransformDepth)),
 		axisX: authoredAxisX(el),
 		axisZ: authoredAxisZ(el)
 	})), translationFrame(specLocalShift(el)));
@@ -28857,10 +29018,10 @@ function civilNodeFrame(el, cumulativeFrame) {
 *  `axisZ` are recomputed from the element's world frame relative to its spatial
 *  container. Stair/ramp flights, which carry their own per-flight frame, are
 *  each re-placed the same way. */
-function rotatedRoutedInput(flatInput, el, cumulativeFrame, spatialFrame) {
+function rotatedRoutedInput(flatInput, el, cumulativeFrame, occurrenceTransformDepth, spatialFrame) {
 	const toSpatial = frameInverse(spatialFrame);
 	if (Array.isArray(flatInput["flights"])) {
-		const elementFrame = elementBodyFrame(el, cumulativeFrame);
+		const elementFrame = elementBodyFrame(el, cumulativeFrame, occurrenceTransformDepth);
 		const flights = flatInput["flights"];
 		return {
 			...flatInput,
@@ -28901,7 +29062,7 @@ function rotatedRoutedInput(flatInput, el, cumulativeFrame, spatialFrame) {
 			})
 		};
 	}
-	const placed = decomposeFrame(frameMul(toSpatial, elementBodyFrame(el, cumulativeFrame)));
+	const placed = decomposeFrame(frameMul(toSpatial, elementBodyFrame(el, cumulativeFrame, occurrenceTransformDepth)));
 	return {
 		...flatInput,
 		origin: placed.origin,
@@ -29066,6 +29227,7 @@ function familiesToBim(root, options) {
 	const walk = (el, state) => {
 		const rotatedHere = state.rotated || hasRotateOp(el);
 		const cumulativeTranslationHere = addTranslation(state.cumulativeTranslation, authoredTranslation(el));
+		const occurrenceTransformDepthHere = state.occurrenceTransformDepth + el.localTransforms.length;
 		const cumulativeFrameHere = frameMul(state.cumulativeFrame, frameFromOps(el.localTransforms));
 		let proxiedHere = false;
 		let nextRotated = rotatedHere;
@@ -29129,7 +29291,7 @@ function familiesToBim(root, options) {
 			}
 			const routedInput = el.semantics?.kind === "product" ? civilProductSpecInput(el) : ("input" in route ? route.input : specInput)(el);
 			const rotatedBase = Array.isArray(routedInput["flights"]) ? specInput(el) : routedInput;
-			const placedInput = rotatedHere ? rotatedRoutedInput(rotatedBase, el, cumulativeFrameHere, state.spatialFrame) : usesAuthoredCivilHierarchy ? relativeSpecInput(routedInput, state.projectedSpatialTranslation) : routedInput;
+			const placedInput = rotatedHere ? rotatedRoutedInput(rotatedBase, el, cumulativeFrameHere, occurrenceTransformDepthHere, state.spatialFrame) : usesAuthoredCivilHierarchy ? relativeSpecInput(routedInput, state.projectedSpatialTranslation) : routedInput;
 			const parsed = route.parse(placedInput);
 			if (!parsed.ok) return parsed;
 			const added = route.add(model, parsed.value, el.keyPath);
@@ -29166,6 +29328,7 @@ function familiesToBim(root, options) {
 				rotated: nextRotated,
 				cumulativeTranslation: cumulativeTranslationHere,
 				projectedSpatialTranslation: nextProjectedSpatialTranslation,
+				occurrenceTransformDepth: occurrenceTransformDepthHere,
 				cumulativeFrame: nextCumulativeFrame,
 				spatialFrame: nextSpatialFrame
 			});
@@ -29179,6 +29342,7 @@ function familiesToBim(root, options) {
 		rotated: false,
 		cumulativeTranslation: ZERO_TRANSLATION,
 		projectedSpatialTranslation: ZERO_TRANSLATION,
+		occurrenceTransformDepth: 0,
 		cumulativeFrame: IDENTITY_FRAME,
 		spatialFrame: IDENTITY_FRAME
 	});

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "brepjs-bim",
-  "version": "0.23.0",
+  "version": "0.23.1",
   "description": "BIM layer for brepjs — IFC4-aligned parametric building elements",
   "keywords": [
     "bim",