@qe-libs/rena-wasm 0.1.11 → 0.1.12

package/README.md

@@ -1,11 +1,13 @@
 # @qe-libs/rena-wasm
 
-JavaScript/WebAssembly ENA pipeline — thin orchestration layer over
-[`@qe-libs/libqe-wasm`](../../libqe/wasm/README.md).
+JavaScript/WebAssembly ENA pipeline.
 
 Handles data parsing, unit/conversation grouping, and the full
 accumulate → normalize → center → rotate → project → node-positions pipeline.
-All math is delegated to the libqe WASM module; no C++ compilation required here.
+The math runs in WebAssembly: the ENA model code (rotations, node positions,
+window estimation) is libena, rENA's C++ layer, compiled into this package
+(`dist/libena.*`), and generic numerics and accumulation come from
+[`@qe-libs/libqe-wasm`](../../libqe/wasm/README.md).
 
 ---
 
@@ -324,9 +326,16 @@ automatically unless `factorLevels` is provided explicitly.
 
 ```bash
 npm install
+npm run build   # compiles dist/libena.* — needs emscripten (emcc) and conan
 npm test
 ```
 
+`npm run build` (`scripts/build.sh`) vendors the libena and libqe headers into
+`include/` (`../scripts/sync-headers.sh`: libena from `inst/include`, libqe from
+Conan), gets Armadillo from Conan, and builds with `emcmake`.
+`test/libena.test.js` checks that this module and libqe-wasm's own copy of the
+ENA functions give identical results.
+
 Tests cover the simple windowed pipeline (`test/ena.test.js`) and the
 context-tensor path (`test/tensor.test.js`), and check against R: infinite and
 time-based windows, flexible horizons and custom rotations

package/dist/libena.js

@@ -0,0 +1,2 @@
+async function LibENA(moduleArg={}){var Module=moduleArg;var ENVIRONMENT_IS_WEB=!!globalThis.window;var ENVIRONMENT_IS_WORKER=!!globalThis.WorkerGlobalScope;var ENVIRONMENT_IS_NODE=globalThis.process?.versions?.node&&globalThis.process?.type!="renderer";if(ENVIRONMENT_IS_NODE){const{createRequire}=await(import("node:module"));var require=createRequire(import.meta.url)}var programArgs=[];var thisProgram="./this.program";var quit_=(status,toThrow)=>{throw toThrow};var _scriptName=import.meta.url;var scriptDirectory="";function locateFile(path){if(Module["locateFile"]){return Module["locateFile"](path,scriptDirectory)}return scriptDirectory+path}var readAsync,readBinary;if(ENVIRONMENT_IS_NODE){var fs=require("node:fs");if(_scriptName.startsWith("file:")){scriptDirectory=require("node:path").dirname(require("node:url").fileURLToPath(_scriptName))+"/"}readBinary=filename=>{filename=isFileURI(filename)?new URL(filename):filename;var ret=fs.readFileSync(filename);return ret};readAsync=async(filename,binary=true)=>{filename=isFileURI(filename)?new URL(filename):filename;var ret=fs.readFileSync(filename,binary?undefined:"utf8");return ret};if(process.argv.length>1){thisProgram=process.argv[1].replace(/\\/g,"/")}programArgs=process.argv.slice(2);quit_=(status,toThrow)=>{process.exitCode=status;throw toThrow}}else if(ENVIRONMENT_IS_WEB||ENVIRONMENT_IS_WORKER){try{scriptDirectory=new URL(".",_scriptName).href}catch{}{readAsync=async url=>{var response=await fetch(url,{credentials:"same-origin"});if(response.ok){return response.arrayBuffer()}throw new Error(response.status+" : "+response.url)}}}else{}var out=console.log.bind(console);var err=console.error.bind(console);var wasmBinary;var ABORT=false;var isFileURI=filename=>filename.startsWith("file://");class EmscriptenEH{}class EmscriptenSjLj extends EmscriptenEH{}var runtimeInitialized=false;function getMemoryBuffer(){return wasmMemory.buffer}function updateMemoryViews(){if(HEAP8?.buffer?.resizable)return;var b=getMemoryBuffer();HEAP8=new Int8Array(b);HEAP16=new Int16Array(b);HEAPU8=new Uint8Array(b);HEAPU16=new Uint16Array(b);HEAP32=new Int32Array(b);HEAPU32=new Uint32Array(b);HEAPF32=new Float32Array(b);HEAPF64=new Float64Array(b);HEAP64=new BigInt64Array(b);HEAPU64=new BigUint64Array(b)}function preRun(){var preRun=Module["preRun"];if(preRun){if(typeof preRun=="function")preRun=[preRun];onPreRuns.push(...preRun)}callRuntimeCallbacks(onPreRuns)}function initRuntime(){runtimeInitialized=true;wasmExports["C"]()}function postRun(){var postRun=Module["postRun"];if(postRun){if(typeof postRun=="function")postRun=[postRun];onPostRuns.push(...postRun)}callRuntimeCallbacks(onPostRuns)}function abort(what){Module["onAbort"]?.(what);what=`Aborted(${what})`;err(what);ABORT=true;what+=". Build with -sASSERTIONS for more info.";var e=new WebAssembly.RuntimeError(what);throw e}var wasmBinaryFile;function findWasmBinary(){if(Module["locateFile"]){return locateFile("libena.wasm")}return new URL("libena.wasm",import.meta.url).href}function getBinarySync(file){if(readBinary){return readBinary(file)}throw"both async and sync fetching of the wasm failed"}async function getWasmBinary(binaryFile){if(!wasmBinary){try{var response=await readAsync(binaryFile);return new Uint8Array(response)}catch{}}return getBinarySync(binaryFile)}async function instantiateArrayBuffer(binaryFile,imports){try{var binary=await getWasmBinary(binaryFile);var instance=await WebAssembly.instantiate(binary,imports);return instance}catch(reason){err(`failed to asynchronously prepare wasm: ${reason}`);abort(reason)}}async function instantiateAsync(binary,binaryFile,imports){if(!binary&&!ENVIRONMENT_IS_NODE){try{var response=fetch(binaryFile,{credentials:"same-origin"});var instantiationResult=await WebAssembly.instantiateStreaming(response,imports);return instantiationResult}catch(reason){err(`wasm streaming compile failed: ${reason}`);err("falling back to ArrayBuffer instantiation")}}return instantiateArrayBuffer(binaryFile,imports)}function getWasmImports(){var imports={a:wasmImports};return imports}async function createWasm(){function receiveInstance(instance){wasmExports=instance.exports;assignWasmExports(wasmExports);updateMemoryViews();return wasmExports}function receiveInstantiationResult(result){return receiveInstance(result["instance"])}var info=getWasmImports();var instantiateWasm=Module["instantiateWasm"];if(instantiateWasm){return new Promise(resolve=>{instantiateWasm(info,inst=>resolve(receiveInstance(inst)))})}wasmBinaryFile??=findWasmBinary();var result=await instantiateAsync(wasmBinary,wasmBinaryFile,info);var exports=receiveInstantiationResult(result);return exports}class ExitStatus{name="ExitStatus";constructor(status){this.message=`Program terminated with exit(${status})`;this.status=status}}var HEAP8;var callRuntimeCallbacks=callbacks=>{while(callbacks.length>0){callbacks.shift()(Module)}};var onPostRuns=[];var onPreRuns=[];var noExitRuntime=true;var stackRestore=val=>__emscripten_stack_restore(val);var stackSave=()=>_emscripten_stack_get_current();var HEAPU32;class ExceptionInfo{constructor(excPtr){this.excPtr=excPtr;this.ptr=excPtr-24}set_type(type){HEAPU32[this.ptr+4>>2]=type}get_type(){return HEAPU32[this.ptr+4>>2]}set_destructor(destructor){HEAPU32[this.ptr+8>>2]=destructor}get_destructor(){return HEAPU32[this.ptr+8>>2]}set_caught(caught){caught=caught?1:0;HEAP8[this.ptr+12]=caught}get_caught(){return HEAP8[this.ptr+12]!=0}set_rethrown(rethrown){rethrown=rethrown?1:0;HEAP8[this.ptr+13]=rethrown}get_rethrown(){return HEAP8[this.ptr+13]!=0}init(type,destructor){this.set_adjusted_ptr(0);this.set_type(type);this.set_destructor(destructor)}set_adjusted_ptr(adjustedPtr){HEAPU32[this.ptr+16>>2]=adjustedPtr}get_adjusted_ptr(){return HEAPU32[this.ptr+16>>2]}}var uncaughtExceptionCount=0;var __Unwind_RaiseException=ex=>{abort()};var ___cxa_throw=(ptr,type,destructor)=>{var info=new ExceptionInfo(ptr);info.init(type,destructor);uncaughtExceptionCount++;__Unwind_RaiseException(ptr)};var __abort_js=()=>abort("");var HEAPU8;var AsciiToString=ptr=>{var str="";while(1){var ch=HEAPU8[ptr++];if(!ch)return str;str+=String.fromCharCode(ch)}};var awaitingDependencies={};var registeredTypes={};var typeDependencies={};class BindingError extends Error{constructor(message){super(message);this.name="BindingError"}}var throwBindingError=message=>{throw new BindingError(message)};function sharedRegisterType(rawType,registeredInstance,options={}){var name=registeredInstance.name;if(!rawType){throwBindingError(`type "${name}" must have a positive integer typeid pointer`)}if(registeredTypes.hasOwnProperty(rawType)){if(options.ignoreDuplicateRegistrations){return}else{throwBindingError(`Cannot register type '${name}' twice`)}}registeredTypes[rawType]=registeredInstance;delete typeDependencies[rawType];if(awaitingDependencies.hasOwnProperty(rawType)){var callbacks=awaitingDependencies[rawType];delete awaitingDependencies[rawType];callbacks.forEach(cb=>cb())}}function registerType(rawType,registeredInstance,options={}){return sharedRegisterType(rawType,registeredInstance,options)}var HEAP16;var HEAPU16;var HEAP32;var HEAP64;var HEAPU64;var integerReadValueFromPointer=(name,width,signed)=>{switch(width){case 1:return signed?pointer=>HEAP8[pointer]:pointer=>HEAPU8[pointer];case 2:return signed?pointer=>HEAP16[pointer>>1]:pointer=>HEAPU16[pointer>>1];case 4:return signed?pointer=>HEAP32[pointer>>2]:pointer=>HEAPU32[pointer>>2];case 8:return signed?pointer=>HEAP64[pointer>>3]:pointer=>HEAPU64[pointer>>3];default:throw new TypeError(`invalid integer width (${width}): ${name}`)}};var __embind_register_bigint=(primitiveType,name,size,minRange,maxRange)=>{name=AsciiToString(name);const isUnsignedType=minRange===0n;let fromWireType=value=>value;if(isUnsignedType){const bitSize=size*8;fromWireType=value=>BigInt.asUintN(bitSize,value);maxRange=fromWireType(maxRange)}registerType(primitiveType,{name,fromWireType,toWireType:(destructors,value)=>{if(typeof value=="number"){value=BigInt(value)}return value},readValueFromPointer:integerReadValueFromPointer(name,size,!isUnsignedType),destructorFunction:null})};var __embind_register_bool=(rawType,name,trueValue,falseValue)=>{name=AsciiToString(name);registerType(rawType,{name,fromWireType:function(wt){return!!wt},toWireType:function(destructors,o){return o?trueValue:falseValue},readValueFromPointer:function(pointer){return this.fromWireType(HEAPU8[pointer])},destructorFunction:null})};var emval_freelist=[];var emval_handles=[0,1,,1,null,1,true,1,false,1];var __emval_decref=handle=>{if(handle>9&&0===--emval_handles[handle+1]){var value=emval_handles[handle];emval_handles[handle]=undefined;emval_freelist.push(handle)}};var Emval={toValue:handle=>{if(!handle){throwBindingError(`Cannot use deleted val. handle = ${handle}`)}return emval_handles[handle]},toHandle:value=>{switch(value){case undefined:return 2;case null:return 4;case true:return 6;case false:return 8;default:{const handle=emval_freelist.pop()||emval_handles.length;emval_handles[handle]=value;emval_handles[handle+1]=1;return handle}}}};function readPointer(pointer){return this.fromWireType(HEAPU32[pointer>>2])}var EmValType={name:"emscripten::val",fromWireType:handle=>{var rv=Emval.toValue(handle);__emval_decref(handle);return rv},toWireType:(destructors,value)=>Emval.toHandle(value),readValueFromPointer:readPointer,destructorFunction:null};var __embind_register_emval=rawType=>registerType(rawType,EmValType);var HEAPF32;var HEAPF64;var floatReadValueFromPointer=(name,width)=>{switch(width){case 4:return function(pointer){return this.fromWireType(HEAPF32[pointer>>2])};case 8:return function(pointer){return this.fromWireType(HEAPF64[pointer>>3])};default:throw new TypeError(`invalid float width (${width}): ${name}`)}};var __embind_register_float=(rawType,name,size)=>{name=AsciiToString(name);registerType(rawType,{name,fromWireType:value=>value,toWireType:(destructors,value)=>value,readValueFromPointer:floatReadValueFromPointer(name,size),destructorFunction:null})};var createNamedFunction=(name,func)=>Object.defineProperty(func,"name",{value:name});var runDestructors=destructors=>{while(destructors.length){var ptr=destructors.pop();var del=destructors.pop();del(ptr)}};function usesDestructorStack(argTypes){for(var i=1;i<argTypes.length;++i){if(argTypes[i]!==null&&argTypes[i].destructorFunction===undefined){return true}}return false}function argsUseStackAlloc(argTypes){for(var i=1;i<argTypes.length;++i){if(argTypes[i]!==null&&argTypes[i].argStackAlloc){return true}}return false}function createJsInvoker(argTypes,isClassMethodFunc,returns,isAsync){var needsDestructorStack=usesDestructorStack(argTypes);var argsNeedStack=argsUseStackAlloc(argTypes);var useStackFrame=argsNeedStack&&!isAsync&&!needsDestructorStack;if(argsNeedStack&&!useStackFrame){needsDestructorStack=true}var argCount=argTypes.length-2;var argsList=[];var argsListWired=["fn"];if(isClassMethodFunc){argsListWired.push("thisWired")}for(var i=0;i<argCount;++i){argsList.push(`arg${i}`);argsListWired.push(`arg${i}Wired`)}argsList=argsList.join();argsListWired=argsListWired.join();var invokerFnBody=`return function (${argsList}) {\n`;if(needsDestructorStack){invokerFnBody+="var destructors = [];\n"}if(useStackFrame){invokerFnBody+="var sp = stackSave();\ntry {\n"}var dtorStack=needsDestructorStack?"destructors":"null";var args1=["humanName","throwBindingError","invoker","fn","runDestructors","fromRetWire","toClassParamWire"];if(useStackFrame){args1.push("stackSave","stackRestore")}if(isClassMethodFunc){invokerFnBody+=`var thisWired = toClassParamWire(${dtorStack}, this);\n`}for(var i=0;i<argCount;++i){var argName=`toArg${i}Wire`;invokerFnBody+=`var arg${i}Wired = ${argName}(${dtorStack}, arg${i});\n`;args1.push(argName)}invokerFnBody+=(returns||isAsync?"var rv = ":"")+`invoker(${argsListWired});\n`;if(useStackFrame){invokerFnBody+="} finally {\nstackRestore(sp);\n}\n"}if(needsDestructorStack){invokerFnBody+="runDestructors(destructors);\n"}else{for(var i=isClassMethodFunc?1:2;i<argTypes.length;++i){var paramName=i===1?"thisWired":`arg${i-2}Wired`;if(argTypes[i].destructorFunction!==null){invokerFnBody+=`${paramName}_dtor(${paramName});\n`;args1.push(`${paramName}_dtor`)}}}if(returns){invokerFnBody+="var ret = fromRetWire(rv);\n"+"return ret;\n"}else{}invokerFnBody+="}\n";return new Function(args1,invokerFnBody)}function craftInvokerFunction(humanName,argTypes,classType,cppInvokerFunc,cppTargetFunc,isAsync){var argCount=argTypes.length;if(argCount<2){throwBindingError("argTypes array size mismatch! Must at least get return value and receiver (this) types!")}var isClassMethodFunc=argTypes[1]!==null&&classType!==null;var needsDestructorStack=usesDestructorStack(argTypes);var argsNeedStack=argsUseStackAlloc(argTypes);var useStackFrame=argsNeedStack&&!isAsync&&!needsDestructorStack;if(argsNeedStack&&!useStackFrame){needsDestructorStack=true}var returns=!argTypes[0].isVoid;var retType=argTypes[0];var instType=argTypes[1];var closureArgs=[humanName,throwBindingError,cppInvokerFunc,cppTargetFunc,runDestructors,retType.fromWireType.bind(retType),instType?.toWireType.bind(instType)];if(useStackFrame){closureArgs.push(stackSave,stackRestore)}for(var i=2;i<argCount;++i){var argType=argTypes[i];closureArgs.push(argType.toWireType.bind(argType))}if(!needsDestructorStack){for(var i=isClassMethodFunc?1:2;i<argTypes.length;++i){if(argTypes[i].destructorFunction!==null){closureArgs.push(argTypes[i].destructorFunction)}}}let invokerFactory=createJsInvoker(argTypes,isClassMethodFunc,returns,isAsync);var invokerFn=invokerFactory(...closureArgs);return createNamedFunction(humanName,invokerFn)}var ensureOverloadTable=(proto,methodName,humanName)=>{if(undefined===proto[methodName].overloadTable){var prevFunc=proto[methodName];proto[methodName]=function(...args){if(!proto[methodName].overloadTable.hasOwnProperty(args.length)){throwBindingError(`Function '${humanName}' called with an invalid number of arguments (${args.length}) - expects one of (${proto[methodName].overloadTable})!`)}return proto[methodName].overloadTable[args.length].apply(this,args)};proto[methodName].overloadTable=[];proto[methodName].overloadTable[prevFunc.argCount]=prevFunc}};var exposePublicSymbol=(name,value,numArguments)=>{if(Module.hasOwnProperty(name)){if(undefined===numArguments||undefined!==Module[name].overloadTable&&undefined!==Module[name].overloadTable[numArguments]){throwBindingError(`Cannot register public name '${name}' twice`)}ensureOverloadTable(Module,name,name);if(Module[name].overloadTable.hasOwnProperty(numArguments)){throwBindingError(`Cannot register multiple overloads of a function with the same number of arguments (${numArguments})!`)}Module[name].overloadTable[numArguments]=value}else{Module[name]=value;Module[name].argCount=numArguments}};var heap32VectorToArray=(count,firstElement)=>{var array=[];for(var i=0;i<count;i++){array.push(HEAPU32[firstElement+i*4>>2])}return array};class InternalError extends Error{constructor(message){super(message);this.name="InternalError"}}var throwInternalError=message=>{throw new InternalError(message)};var replacePublicSymbol=(name,value,numArguments)=>{if(!Module.hasOwnProperty(name)){throwInternalError("Replacing nonexistent public symbol")}if(undefined!==Module[name].overloadTable&&undefined!==numArguments){Module[name].overloadTable[numArguments]=value}else{Module[name]=value;Module[name].argCount=numArguments}};var wasmTableMirror=[];var getWasmTableEntry=funcPtr=>{var func=wasmTableMirror[funcPtr];if(!func){wasmTableMirror[funcPtr]=func=wasmTable.get(funcPtr)}return func};var embind__requireFunction=(signature,rawFunction,isAsync=false)=>{signature=AsciiToString(signature);function makeDynCaller(){var rtn=getWasmTableEntry(rawFunction);return rtn}var fp=makeDynCaller();if(typeof fp!="function"){throwBindingError(`unknown function pointer with signature ${signature}: ${rawFunction}`)}return fp};class UnboundTypeError extends Error{}var getTypeName=type=>{var ptr=___getTypeName(type);var rv=AsciiToString(ptr);_free(ptr);return rv};var throwUnboundTypeError=(message,types)=>{var unboundTypes=[];var seen={};function visit(type){if(seen[type]){return}if(registeredTypes[type]){return}if(typeDependencies[type]){typeDependencies[type].forEach(visit);return}unboundTypes.push(type);seen[type]=true}types.forEach(visit);throw new UnboundTypeError(`${message}: `+unboundTypes.map(getTypeName).join([", "]))};var whenDependentTypesAreResolved=(myTypes,dependentTypes,getTypeConverters)=>{myTypes.forEach(type=>typeDependencies[type]=dependentTypes);function onComplete(typeConverters){var myTypeConverters=getTypeConverters(typeConverters);if(myTypeConverters.length!==myTypes.length){throwInternalError("Mismatched type converter count")}for(var i=0;i<myTypes.length;++i){registerType(myTypes[i],myTypeConverters[i])}}var typeConverters=new Array(dependentTypes.length);var unregisteredTypes=[];var registered=0;for(let[i,dt]of dependentTypes.entries()){if(registeredTypes.hasOwnProperty(dt)){typeConverters[i]=registeredTypes[dt]}else{unregisteredTypes.push(dt);if(!awaitingDependencies.hasOwnProperty(dt)){awaitingDependencies[dt]=[]}awaitingDependencies[dt].push(()=>{typeConverters[i]=registeredTypes[dt];++registered;if(registered===unregisteredTypes.length){onComplete(typeConverters)}})}}if(0===unregisteredTypes.length){onComplete(typeConverters)}};var getFunctionName=signature=>{signature=signature.trim();const argsIndex=signature.indexOf("(");if(argsIndex===-1)return signature;return signature.slice(0,argsIndex)};var __embind_register_function=(name,argCount,rawArgTypesAddr,signature,rawInvoker,fn,isAsync,isNonnullReturn)=>{var argTypes=heap32VectorToArray(argCount,rawArgTypesAddr);name=AsciiToString(name);name=getFunctionName(name);rawInvoker=embind__requireFunction(signature,rawInvoker,isAsync);exposePublicSymbol(name,function(){throwUnboundTypeError(`Cannot call ${name} due to unbound types`,argTypes)},argCount-1);whenDependentTypesAreResolved([],argTypes,argTypes=>{var invokerArgsArray=[argTypes[0],null].concat(argTypes.slice(1));replacePublicSymbol(name,craftInvokerFunction(name,invokerArgsArray,null,rawInvoker,fn,isAsync),argCount-1);return[]})};var __embind_register_integer=(primitiveType,name,size,minRange,maxRange)=>{name=AsciiToString(name);const isUnsignedType=minRange===0;let fromWireType=value=>value;if(isUnsignedType){var bitshift=32-8*size;fromWireType=value=>value<<bitshift>>>bitshift;maxRange=fromWireType(maxRange)}registerType(primitiveType,{name,fromWireType,toWireType:(destructors,value)=>value,readValueFromPointer:integerReadValueFromPointer(name,size,minRange!==0),destructorFunction:null})};var __embind_register_memory_view=(rawType,dataTypeIndex,name)=>{var typeMapping=[Int8Array,Uint8Array,Int16Array,Uint16Array,Int32Array,Uint32Array,Float32Array,Float64Array,BigInt64Array,BigUint64Array];var TA=typeMapping[dataTypeIndex];function decodeMemoryView(handle){var size=HEAPU32[handle>>2];var data=HEAPU32[handle+4>>2];return new TA(HEAP8.buffer,data,size)}name=AsciiToString(name);registerType(rawType,{name,fromWireType:decodeMemoryView,readValueFromPointer:decodeMemoryView},{ignoreDuplicateRegistrations:true})};var stringToUTF8Array=(str,heap,outIdx,maxBytesToWrite)=>{if(!(maxBytesToWrite>0))return 0;var startIdx=outIdx;var endIdx=outIdx+maxBytesToWrite-1;for(var i=0;i<str.length;++i){var u=str.codePointAt(i);if(u<=127){if(outIdx>=endIdx)break;heap[outIdx++]=u}else if(u<=2047){if(outIdx+1>=endIdx)break;heap[outIdx++]=192|u>>6;heap[outIdx++]=128|u&63}else if(u<=65535){if(outIdx+2>=endIdx)break;heap[outIdx++]=224|u>>12;heap[outIdx++]=128|u>>6&63;heap[outIdx++]=128|u&63}else{if(outIdx+3>=endIdx)break;heap[outIdx++]=240|u>>18;heap[outIdx++]=128|u>>12&63;heap[outIdx++]=128|u>>6&63;heap[outIdx++]=128|u&63;i++}}heap[outIdx]=0;return outIdx-startIdx};var stringToUTF8=(str,outPtr,maxBytesToWrite)=>stringToUTF8Array(str,HEAPU8,outPtr,maxBytesToWrite);var lengthBytesUTF8=str=>{var len=0;for(var i=0;i<str.length;++i){var c=str.charCodeAt(i);if(c<=127){len++}else if(c<=2047){len+=2}else if(c>=55296&&c<=57343){len+=4;++i}else{len+=3}}return len};var UTF8Decoder=globalThis.TextDecoder&&new TextDecoder;var findStringEnd=(heapOrArray,idx,maxBytesToRead,ignoreNul)=>{var maxIdx=idx+maxBytesToRead;if(ignoreNul)return maxIdx;while(heapOrArray[idx]&&!(idx>=maxIdx))++idx;return idx};var UTF8ArrayToString=(heapOrArray,idx=0,maxBytesToRead,ignoreNul)=>{var endPtr=findStringEnd(heapOrArray,idx,maxBytesToRead,ignoreNul);if(endPtr-idx>16&&heapOrArray.buffer&&UTF8Decoder){return UTF8Decoder.decode(heapOrArray.subarray(idx,endPtr))}var str="";while(idx<endPtr){var u0=heapOrArray[idx++];if(!(u0&128)){str+=String.fromCharCode(u0);continue}var u1=heapOrArray[idx++]&63;if((u0&224)==192){str+=String.fromCharCode((u0&31)<<6|u1);continue}var u2=heapOrArray[idx++]&63;if((u0&240)==224){u0=(u0&15)<<12|u1<<6|u2}else{u0=(u0&7)<<18|u1<<12|u2<<6|heapOrArray[idx++]&63}if(u0<65536){str+=String.fromCharCode(u0)}else{var ch=u0-65536;str+=String.fromCharCode(55296|ch>>10,56320|ch&1023)}}return str};var UTF8ToString=(ptr,maxBytesToRead,ignoreNul)=>ptr?UTF8ArrayToString(HEAPU8,ptr,maxBytesToRead,ignoreNul):"";var __embind_register_std_string=(rawType,name)=>{name=AsciiToString(name);var stdStringIsUTF8=true;registerType(rawType,{name,fromWireType(value){var length=HEAPU32[value>>2];var payload=value+4;var str;if(stdStringIsUTF8){str=UTF8ToString(payload,length,true)}else{str="";for(var i=0;i<length;++i){str+=String.fromCharCode(HEAPU8[payload+i])}}_free(value);return str},toWireType(destructors,value){if(value instanceof ArrayBuffer){value=new Uint8Array(value)}var length;var valueIsOfTypeString=typeof value=="string";if(!(valueIsOfTypeString||ArrayBuffer.isView(value)&&value.BYTES_PER_ELEMENT==1)){throwBindingError("Cannot pass non-string to std::string")}if(stdStringIsUTF8&&valueIsOfTypeString){length=lengthBytesUTF8(value)}else{length=value.length}var base=_malloc(4+length+1);var ptr=base+4;HEAPU32[base>>2]=length;if(valueIsOfTypeString){if(stdStringIsUTF8){stringToUTF8(value,ptr,length+1)}else{for(var i=0;i<length;++i){var charCode=value.charCodeAt(i);if(charCode>255){_free(base);throwBindingError("String has UTF-16 code units that do not fit in 8 bits")}HEAPU8[ptr+i]=charCode}}}else{HEAPU8.set(value,ptr)}if(destructors!==null){destructors.push(_free,base)}return base},readValueFromPointer:readPointer,destructorFunction(ptr){_free(ptr)}})};var UTF16Decoder=globalThis.TextDecoder?new TextDecoder("utf-16le"):undefined;var UTF16ToString=(ptr,maxBytesToRead,ignoreNul)=>{var idx=ptr>>1;var endIdx=findStringEnd(HEAPU16,idx,maxBytesToRead/2,ignoreNul);if(endIdx-idx>16&&UTF16Decoder)return UTF16Decoder.decode(HEAPU16.subarray(idx,endIdx));var str="";for(var i=idx;i<endIdx;++i){var codeUnit=HEAPU16[i];str+=String.fromCharCode(codeUnit)}return str};var stringToUTF16=(str,outPtr,maxBytesToWrite=2147483647)=>{if(maxBytesToWrite<2)return 0;maxBytesToWrite-=2;var startPtr=outPtr;var numCharsToWrite=maxBytesToWrite<str.length*2?maxBytesToWrite/2:str.length;for(var i=0;i<numCharsToWrite;++i){var codeUnit=str.charCodeAt(i);HEAP16[outPtr>>1]=codeUnit;outPtr+=2}HEAP16[outPtr>>1]=0;return outPtr-startPtr};var lengthBytesUTF16=str=>str.length*2;var UTF32ToString=(ptr,maxBytesToRead,ignoreNul)=>{var str="";var startIdx=ptr>>2;for(var i=0;!(i>=maxBytesToRead/4);i++){var utf32=HEAPU32[startIdx+i];if(!utf32&&!ignoreNul)break;str+=String.fromCodePoint(utf32)}return str};var stringToUTF32=(str,outPtr,maxBytesToWrite=2147483647)=>{if(maxBytesToWrite<4)return 0;var startPtr=outPtr;var endPtr=startPtr+maxBytesToWrite-4;for(var i=0;i<str.length;++i){var codePoint=str.codePointAt(i);if(codePoint>65535){i++}HEAP32[outPtr>>2]=codePoint;outPtr+=4;if(outPtr+4>endPtr)break}HEAP32[outPtr>>2]=0;return outPtr-startPtr};var lengthBytesUTF32=str=>{var len=0;for(var i=0;i<str.length;++i){var codePoint=str.codePointAt(i);if(codePoint>65535){i++}len+=4}return len};var __embind_register_std_wstring=(rawType,charSize,name)=>{name=AsciiToString(name);var decodeString,encodeString,lengthBytesUTF;if(charSize===2){decodeString=UTF16ToString;encodeString=stringToUTF16;lengthBytesUTF=lengthBytesUTF16}else{decodeString=UTF32ToString;encodeString=stringToUTF32;lengthBytesUTF=lengthBytesUTF32}registerType(rawType,{name,fromWireType:value=>{var length=HEAPU32[value>>2];var str=decodeString(value+4,length*charSize,true);_free(value);return str},toWireType:(destructors,value)=>{if(!(typeof value=="string")){throwBindingError(`Cannot pass non-string to C++ string type ${name}`)}var length=lengthBytesUTF(value);var ptr=_malloc(4+length+charSize);HEAPU32[ptr>>2]=length/charSize;encodeString(value,ptr+4,length+charSize);if(destructors!==null){destructors.push(_free,ptr)}return ptr},readValueFromPointer:readPointer,destructorFunction(ptr){_free(ptr)}})};var __embind_register_void=(rawType,name)=>{name=AsciiToString(name);registerType(rawType,{isVoid:true,name,fromWireType:()=>undefined,toWireType:(destructors,o)=>undefined})};var emval_methodCallers=[];var emval_addMethodCaller=caller=>{var id=emval_methodCallers.length;emval_methodCallers.push(caller);return id};var requireRegisteredType=(rawType,humanName)=>{var impl=registeredTypes[rawType];if(undefined===impl){throwBindingError(`${humanName} has unknown type ${getTypeName(rawType)}`)}return impl};var emval_lookupTypes=(argCount,argTypes)=>{var a=new Array(argCount);for(var i=0;i<argCount;++i){a[i]=requireRegisteredType(HEAPU32[argTypes+i*4>>2],`parameter ${i}`)}return a};var emval_returnValue=(toReturnWire,destructorsRef,handle)=>{var destructors=[];var result=toReturnWire(destructors,handle);if(destructors.length){HEAPU32[destructorsRef>>2]=Emval.toHandle(destructors)}return result};var emval_symbols={};var getStringOrSymbol=address=>{var symbol=emval_symbols[address];if(symbol===undefined){return AsciiToString(address)}return symbol};var __emval_create_invoker=(argCount,argTypesPtr,kind)=>{var GenericWireTypeSize=8;var[retType,...argTypes]=emval_lookupTypes(argCount,argTypesPtr);var toReturnWire=retType.toWireType.bind(retType);var argFromPtr=argTypes.map(type=>type.readValueFromPointer.bind(type));argCount--;var captures={toValue:Emval.toValue};var args=argFromPtr.map((argFromPtr,i)=>{var captureName=`argFromPtr${i}`;captures[captureName]=argFromPtr;return`${captureName}(args${i?"+"+i*GenericWireTypeSize:""})`});var functionBody;switch(kind){case 0:functionBody="toValue(handle)";break;case 2:functionBody="new (toValue(handle))";break;case 3:functionBody="";break;case 1:captures["getStringOrSymbol"]=getStringOrSymbol;functionBody="toValue(handle)[getStringOrSymbol(methodName)]";break}functionBody+=`(${args})`;if(!retType.isVoid){captures["toReturnWire"]=toReturnWire;captures["emval_returnValue"]=emval_returnValue;functionBody=`return emval_returnValue(toReturnWire, destructorsRef, ${functionBody})`}functionBody=`return function (handle, methodName, destructorsRef, args) {\n${functionBody}\n}`;var invokerFunction=new Function(Object.keys(captures),functionBody)(...Object.values(captures));var functionName=`methodCaller<(${argTypes.map(t=>t.name)}) => ${retType.name}>`;return emval_addMethodCaller(createNamedFunction(functionName,invokerFunction))};var __emval_get_property=(handle,key)=>{handle=Emval.toValue(handle);key=Emval.toValue(key);return Emval.toHandle(handle[key])};var __emval_incref=handle=>{if(handle>9){emval_handles[handle+1]+=1}};var __emval_invoke=(caller,handle,methodName,destructorsRef,args)=>emval_methodCallers[caller](handle,methodName,destructorsRef,args);var __emval_is_string=handle=>{handle=Emval.toValue(handle);return typeof handle=="string"};var __emval_new_array=()=>Emval.toHandle([]);var __emval_new_cstring=v=>Emval.toHandle(getStringOrSymbol(v));var __emval_new_object=()=>Emval.toHandle({});var __emval_run_destructors=handle=>{var destructors=Emval.toValue(handle);runDestructors(destructors);__emval_decref(handle)};var __emval_set_property=(handle,key,value)=>{handle=Emval.toValue(handle);key=Emval.toValue(key);value=Emval.toValue(value);handle[key]=value};var __tzset_js=(timezone,daylight,std_name,dst_name)=>{var currentYear=(new Date).getFullYear();var winter=new Date(currentYear,0,1);var summer=new Date(currentYear,6,1);var winterOffset=winter.getTimezoneOffset();var summerOffset=summer.getTimezoneOffset();var stdTimezoneOffset=Math.max(winterOffset,summerOffset);HEAPU32[timezone>>2]=stdTimezoneOffset*60;HEAP32[daylight>>2]=Number(winterOffset!=summerOffset);var extractZone=timezoneOffset=>{var sign=timezoneOffset>=0?"-":"+";var absOffset=Math.abs(timezoneOffset);var hours=String(Math.floor(absOffset/60)).padStart(2,"0");var minutes=String(absOffset%60).padStart(2,"0");return`UTC${sign}${hours}${minutes}`};var winterName=extractZone(winterOffset);var summerName=extractZone(summerOffset);if(summerOffset<winterOffset){stringToUTF8(winterName,std_name,17);stringToUTF8(summerName,dst_name,17)}else{stringToUTF8(winterName,dst_name,17);stringToUTF8(summerName,std_name,17)}};var getHeapMax=()=>2147483648;var alignMemory=(size,alignment)=>Math.ceil(size/alignment)*alignment;var growMemory=size=>{var oldHeapSize=wasmMemory.buffer.byteLength;var pages=(size-oldHeapSize+65535)/65536|0;try{wasmMemory.grow(pages);updateMemoryViews();return 1}catch(e){}};var _emscripten_resize_heap=requestedSize=>{var oldSize=HEAPU8.length;requestedSize>>>=0;var maxHeapSize=getHeapMax();if(requestedSize>maxHeapSize){return false}for(var cutDown=1;cutDown<=4;cutDown*=2){var overGrownHeapSize=oldSize*(1+.2/cutDown);overGrownHeapSize=Math.min(overGrownHeapSize,requestedSize+100663296);var newSize=Math.min(maxHeapSize,alignMemory(Math.max(requestedSize,overGrownHeapSize),65536));var replacement=growMemory(newSize);if(replacement){return true}}return false};var ENV={};var getExecutableName=()=>thisProgram;var getEnvStrings=()=>{if(!getEnvStrings.strings){var lang=(globalThis.navigator?.language??"C").replace("-","_")+".UTF-8";var env={USER:"web_user",LOGNAME:"web_user",PATH:"/",PWD:"/",HOME:"/home/web_user",LANG:lang,_:getExecutableName()};for(var x in ENV){if(ENV[x]===undefined)delete env[x];else env[x]=ENV[x]}var strings=[];for(var x in env){strings.push(`${x}=${env[x]}`)}getEnvStrings.strings=strings}return getEnvStrings.strings};var _environ_get=(__environ,environ_buf)=>{var bufSize=0;var envp=0;for(var string of getEnvStrings()){var ptr=environ_buf+bufSize;HEAPU32[__environ+envp>>2]=ptr;bufSize+=stringToUTF8(string,ptr,Infinity)+1;envp+=4}return 0};var _environ_sizes_get=(penviron_count,penviron_buf_size)=>{var strings=getEnvStrings();HEAPU32[penviron_count>>2]=strings.length;var bufSize=0;for(var string of strings){bufSize+=lengthBytesUTF8(string)+1}HEAPU32[penviron_buf_size>>2]=bufSize;return 0};{if(Module["noExitRuntime"])noExitRuntime=Module["noExitRuntime"];if(Module["print"])out=Module["print"];if(Module["printErr"])err=Module["printErr"];if(Module["arguments"])programArgs=Module["arguments"];if(Module["thisProgram"])thisProgram=Module["thisProgram"];var preInit=Module["preInit"];if(preInit){if(typeof preInit=="function")Module["preInit"]=preInit=[preInit];while(preInit.length>0){preInit.shift()()}}}var ___getTypeName,_free,_malloc,__emscripten_stack_restore,_emscripten_stack_get_current,memory,__indirect_function_table,wasmMemory,wasmTable;function assignWasmExports(wasmExports){___getTypeName=wasmExports["D"];_free=wasmExports["E"];_malloc=wasmExports["G"];__emscripten_stack_restore=wasmExports["H"];_emscripten_stack_get_current=wasmExports["I"];memory=wasmMemory=wasmExports["B"];__indirect_function_table=wasmTable=wasmExports["F"]}var wasmImports={b:___cxa_throw,u:__abort_js,q:__embind_register_bigint,y:__embind_register_bool,w:__embind_register_emval,p:__embind_register_float,g:__embind_register_function,l:__embind_register_integer,h:__embind_register_memory_view,x:__embind_register_std_string,o:__embind_register_std_wstring,z:__embind_register_void,f:__emval_create_invoker,a:__emval_decref,m:__emval_get_property,k:__emval_incref,e:__emval_invoke,A:__emval_is_string,j:__emval_new_array,c:__emval_new_cstring,n:__emval_new_object,d:__emval_run_destructors,i:__emval_set_property,t:__tzset_js,v:_emscripten_resize_heap,r:_environ_get,s:_environ_sizes_get};async function run(){preRun();var setStatus=Module["setStatus"];if(setStatus){setStatus("Running...");await new Promise(resolve=>setTimeout(resolve,1));setTimeout(setStatus,1,"")}if(ABORT)return;initRuntime();Module["onRuntimeInitialized"]?.();postRun()}var wasmExports;wasmExports=await createWasm();await run();
+;return Module}export default LibENA;

package/package.json

@@ -1,13 +1,16 @@
 {
   "name": "@qe-libs/rena-wasm",
-  "version": "0.1.11",
-  "description": "JavaScript/WebAssembly ENA pipeline — thin orchestration layer over @qe-libs/libqe-wasm",
+  "version": "0.1.12",
+  "description": "JavaScript/WebAssembly ENA pipeline — libena (ENA rotations, node positions) compiled to WASM, over @qe-libs/libqe-wasm",
   "main": "src/index.js",
   "type": "module",
   "files": [
+    "dist/libena.js",
+    "dist/libena.wasm",
     "src/"
   ],
   "scripts": {
+    "build": "sh scripts/build.sh",
     "test": "node --experimental-vm-modules node_modules/.bin/jest"
   },
   "publishConfig": {

package/src/data.js

@@ -17,6 +17,32 @@ export function rowKey(row, cols) {
     return cols.map(c => row[c]).join('__');
 }
 
+// rowKey() is the label callers see (and pass back as unitsUsed), but joining
+// with '__' is ambiguous: {a: 'x__y', b: 'z'} and {a: 'x', b: 'y__z'} both give
+// 'x__y__z'. Track the exact values behind each label and refuse to merge
+// different ones under the same label.
+function checkedKey(seen, row, cols, what) {
+    const key   = rowKey(row, cols);
+    const exact = JSON.stringify(cols.map(c => row[c] ?? null));
+    const prev  = seen.get(key);
+    if (prev === undefined) seen.set(key, exact);
+    else if (prev !== exact) {
+        throw new Error(`Different ${what} values both map to "${key}" because a value ` +
+                        `contains "__"; rename the values or use different ${what} columns.`);
+    }
+    return key;
+}
+
+function requireColumns(kind, cols, available) {
+    if (!Array.isArray(cols) || cols.length === 0 || !cols.every(c => typeof c === 'string')) {
+        throw new Error(`${kind} must be a non-empty array of column names`);
+    }
+    const missing = cols.filter(c => !available.has(c));
+    if (missing.length) {
+        throw new Error(`${kind}: column(s) not found in the data: ${missing.join(', ')}`);
+    }
+}
+
 /**
  * Parse tabular data into ENA-ready structures.
  *
@@ -55,9 +81,19 @@ export function parseData(rows, codes, unitCols, convoCols, unitsUsed = null) {
     const convoGroups = new Map();   // convoIdx → [rowIdx, ...]
 
     // Metadata columns: every column that is NOT a code, unit, or convo column.
-    const excludedCols = new Set([...codes, ...unitCols, ...convoCols]);
     const allCols      = nRows > 0 ? Object.keys(rows[0]) : [];
+    // A misspelled column used to read as undefined everywhere: all-zero codes,
+    // or every row in one unit called "undefined".
+    if (nRows > 0) {
+        const available = new Set(allCols);
+        requireColumns('codes', codes, available);
+        requireColumns('units', unitCols, available);
+        requireColumns('conversations', convoCols, available);
+    }
+    const excludedCols = new Set([...codes, ...unitCols, ...convoCols]);
     const metaCols     = allCols.filter(c => !excludedCols.has(c));
+    const unitSeen     = new Map();   // label → exact values (see checkedKey)
+    const convoSeen    = new Map();
 
     // Per-unit metadata rows — populated on first occurrence of each unit.
     const metaData = [];
@@ -67,13 +103,20 @@ export function parseData(rows, codes, unitCols, convoCols, unitsUsed = null) {
     for (let i = 0; i < nRows; i++) {
         const row = rows[i];
 
-        // Code values
+        // Code values: numbers, numeric strings or booleans. A missing or
+        // non-numeric value is an error rather than a silent 0.
         for (let c = 0; c < nCodes; c++) {
-            codeMatrix[i * nCodes + c] = Number(row[codes[c]]) || 0;
+            const raw = row[codes[c]];
+            const v   = typeof raw === 'string' && raw.trim() === '' ? NaN : Number(raw);
+            if (!Number.isFinite(v)) {
+                throw new Error(`Code "${codes[c]}" has a missing or non-numeric value ` +
+                                `(${JSON.stringify(raw ?? null)}) in row ${i + 1}`);
+            }
+            codeMatrix[i * nCodes + c] = v;
         }
 
         // Unit index
-        const uKey = rowKey(row, unitCols);
+        const uKey = checkedKey(unitSeen, row, unitCols, 'unit');
         if (used && !used.has(uKey)) {
             unitOf[i] = -1;
         } else if (!unitIndex.has(uKey)) {
@@ -87,7 +130,7 @@ export function parseData(rows, codes, unitCols, convoCols, unitsUsed = null) {
         if (unitOf[i] !== -1) unitOf[i] = unitIndex.get(uKey);
 
         // Conversation index
-        const cKey = rowKey(row, convoCols);
+        const cKey = checkedKey(convoSeen, row, convoCols, 'conversation');
         if (!convoIndex.has(cKey)) convoIndex.set(cKey, convoIndex.size);
         const cIdx = convoIndex.get(cKey);
         convoOf[i] = cIdx;

package/src/index.js

@@ -2,9 +2,11 @@
  * @qe-libs/rena-wasm
  *
  * JavaScript/WebAssembly ENA pipeline.
- * Thin orchestration layer over @qe-libs/libqe-wasm — handles data parsing,
- * unit/conversation grouping, and the full accumulate→normalize→center→
- * rotate→project→node-positions pipeline.
+ * Handles data parsing, unit/conversation grouping, and the full
+ * accumulate→normalize→center→rotate→project→node-positions pipeline.  The
+ * ENA math (rotations, node positions, CCD) is libena, compiled into this
+ * package (dist/libena.*) with tma's accumulation (libtma); generic numerics
+ * come from @qe-libs/libqe-wasm.
  *
  * Output structure mirrors R's ena.set object (without R-specific S3 class
  * attributes and without metadata columns prepended to every matrix).
@@ -41,6 +43,7 @@
  */
 
 import loadLibQE from '@qe-libs/libqe-wasm';
+import createLibENA from '../dist/libena.js';
 import { parseData } from './data.js';
 import {
     sphereNorm, scaleNetworks, center,
@@ -294,18 +297,42 @@ function applyCodeMask(networks, codeMask, nCodes, nUnits, nConnections, ordered
 
 // Rows of a matrix given as an array of row arrays, or as flat row-major data
 // with `cols` columns.
-function matrixRows(m, cols) {
+function matrixRows(m, cols, name = 'rotation matrix') {
     if (m == null) return null;
+    let rows;
     if (Array.isArray(m) && (m.length === 0 || Array.isArray(m[0]) || ArrayBuffer.isView(m[0]))) {
-        return m.map(r => Array.from(r, Number));
+        rows = m.map(r => Array.from(r, Number));
+    } else {
+        const flat = Array.from(m, Number);
+        if (!cols) throw new Error(`a flat ${name} needs its column count`);
+        if (flat.length % cols !== 0) {
+            throw new Error(`${name} has ${flat.length} values, not a multiple of ${cols} columns`);
+        }
+        rows = [];
+        for (let i = 0; i < flat.length; i += cols) rows.push(flat.slice(i, i + cols));
+    }
+    // Ragged rows used to be laid out at a fixed stride: a long row overwrote
+    // the start of the next one, a short one left zeros.
+    const width = cols || (rows[0] ? rows[0].length : 0);
+    for (const r of rows) {
+        if (r.length !== width) throw new Error(`${name} rows must all have ${width} values`);
+        if (!r.every(Number.isFinite)) throw new Error(`${name} contains a non-numeric value`);
     }
-    const flat = Array.from(m, Number);
-    if (!cols) throw new Error('a flat rotation matrix needs its column count');
-    const rows = [];
-    for (let i = 0; i < flat.length; i += cols) rows.push(flat.slice(i, i + cols));
     return rows;
 }
 
+// The libqe statistics read nRows × dims values from `points`; a mismatch
+// used to read past the array (leftover WASM heap memory) — libqe now throws,
+// this gives the reason.
+function checkPoints(name, points, nRows, dims) {
+    for (const [what, n] of [['row count', nRows], ['dims', dims]]) {
+        if (!Number.isInteger(n) || n < 0) throw new Error(`${name}: ${what} must be a non-negative integer`);
+    }
+    if (points == null || points.length !== nRows * dims) {
+        throw new Error(`${name}: expected ${nRows} × ${dims} = ${nRows * dims} values, got ${points?.length ?? 0}`);
+    }
+}
+
 /**
  * Validate a supplied rotation set (= R's ena.rotation.set, e.g. another
  * model's rotation) against this model's codes, and reorder it to this
@@ -330,7 +357,7 @@ function resolveRotationSet(rs, codes, nConnections, ordered) {
         );
     }
     const rotCols = rot[0].length;
-    const nodes   = matrixRows(rs.nodes, rotCols);
+    const nodes   = matrixRows(rs.nodes, rotCols, 'rotationSet.nodes');
     const center  = rs.centerVec ? Array.from(rs.centerVec, Number) : null;
     if (!nodes || nodes.length !== codes.length) {
         throw new Error('rotationSet.nodes must have one row per code');
@@ -585,6 +612,53 @@ function runCustomRotation(qe, rawNetworks, lineWeights, ownCentered, ownCenterV
     });
 }
 
+// ── WASM modules ─────────────────────────────────────────────────────────────
+
+// Functions compiled into dist/libena.*: libena's (inst/include/libena; moved
+// out of libqe in libqe's phase 4a split) and libtma's accumulation (tma's
+// inst/include/libtma; phase 5).  libqe-wasm still carries its own copies
+// until libqe 0.2.0, so these names must route to dist/libena.
+const LIBENA_FUNCTIONS = new Set([
+    'ena_correlation',
+    'node_positions',
+    'directed_node_positions',
+    'directed_node_positions_combine_pairs',
+    'ena_svd',
+    'deflate',
+    'orthogonal_svd',
+    'complete_rotation',
+    'means_rotation',
+    'generalized_means_rotation',
+    'ccd_window',
+    // libtma
+    'connection_matrix',
+    'accumulate_stanza',
+    'row_connections',
+    'rolling_window_sum',
+    'flat_index',
+    'accumulate_unit',
+    'accumulate_unit_with_rows',
+    'accumulate_tensor_unit',
+    'aggregate_row_connections',
+    'finalize_row_connections',
+]);
+
+/**
+ * Load libqe-wasm and libena and present them as one module object: the
+ * pipeline (pipeline.js, tensor.js, ...) takes a single `qe`, so libena's
+ * functions are routed to libena and everything else to libqe-wasm.
+ */
+async function loadModules() {
+    const [libqe, libena] = await Promise.all([loadLibQE(), createLibENA()]);
+    return new Proxy(libqe, {
+        get(target, prop, receiver) {
+            return LIBENA_FUNCTIONS.has(prop)
+                ? libena[prop]
+                : Reflect.get(target, prop, receiver);
+        },
+    });
+}
+
 // ── main factory ─────────────────────────────────────────────────────────────
 
 /**
@@ -593,7 +667,7 @@ function runCustomRotation(qe, rawNetworks, lineWeights, ownCentered, ownCenterV
  * @returns {Promise<{fit: Function, accumulate: Function}>}
  */
 export default async function loadENA() {
-    const qe = await loadLibQE();
+    const qe = await loadModules();
 
     const api = {
         /**
@@ -902,6 +976,7 @@ export default async function loadENA() {
          *   rows=dims, cols=3 — columns: mean, lower CI, upper CI
          */
         confInts(points, nUnits, dims, confLevel = 0.95) {
+            checkPoints('confInts', points, nUnits, dims);
             return qe.mean_ci(points, nUnits, dims, confLevel);
         },
 
@@ -917,6 +992,7 @@ export default async function loadENA() {
          *   rows=dims, cols=2 — columns: lower fence, upper fence
          */
         outlierInts(points, nUnits, dims, iqrFactor = 1.5) {
+            checkPoints('outlierInts', points, nUnits, dims);
             return qe.outlier_ci(points, nUnits, dims, iqrFactor);
         },
 
@@ -933,6 +1009,8 @@ export default async function loadENA() {
          *             U, pvalue_u, effect_r, medians }}
          */
         compareGroups(g1Points, nG1, g2Points, nG2, dims) {
+            checkPoints('compareGroups (group 1)', g1Points, nG1, dims);
+            checkPoints('compareGroups (group 2)', g2Points, nG2, dims);
             return qe.group_stats(g1Points, nG1, dims, g2Points, nG2, dims);
         },
 
@@ -1124,7 +1202,11 @@ export default async function loadENA() {
          */
         defaultTensor,
 
-        /** The underlying libqe WASM module, for custom low-level pipelines. */
+        /**
+         * The underlying WASM functions, for custom low-level pipelines:
+         * libqe-wasm's, with libena's (rotations, node positions, CCD) from
+         * this package.
+         */
         qe,
     };
 

package/src/tensor.js

@@ -82,8 +82,16 @@ export function buildContextLookup(rows, rowIndices, factors, factorLevels) {
         const row = rows[rowIndices[i]];
         times[i] = i;  // default: row index as time
         for (let f = 0; f < nFactor; f++) {
-            const col = factors[f];
-            contextLookup[i * nFactor + f] = factorLevels[col][row[col]] ?? 0;
+            const col    = factors[f];
+            const levels = factorLevels[col];
+            // Own properties only: a value like "constructor" must not resolve
+            // to an inherited Object.prototype member. A value with no level
+            // used to fall back to level 0 silently.
+            if (!levels || !Object.prototype.hasOwnProperty.call(levels, row[col])) {
+                throw new Error(`tensor factor "${col}" has no level for value ` +
+                                `${JSON.stringify(row[col] ?? null)}`);
+            }
+            contextLookup[i * nFactor + f] = levels[row[col]];
         }
     }
 
@@ -221,7 +229,9 @@ export function parseTimes(rows, col, endCol = null) {
 export function inferFactorLevels(rows, factors) {
     const levels = {};
     for (const col of factors) {
-        levels[col] = {};
+        // No prototype: values such as "constructor" or "__proto__" are
+        // ordinary keys here.
+        levels[col] = Object.create(null);
         let idx = 0;
         for (const row of rows) {
             const v = row[col];
@@ -233,6 +243,39 @@ export function inferFactorLevels(rows, factors) {
     return levels;
 }
 
+/**
+ * Check a tensor definition before it reaches libqe, which indexes the
+ * tensor with these values (an index out of range used to read past it).
+ */
+export function validateTensor(dims, tensorData, factors, factorLevels) {
+    if (!Array.isArray(dims) && !ArrayBuffer.isView(dims)) {
+        throw new Error('tensor dims must be an array');
+    }
+    const d = Array.from(dims);
+    if (d.length === 0 || !d.every(x => Number.isInteger(x) && x > 0)) {
+        throw new Error('tensor dims must be positive integers');
+    }
+    const size = d.reduce((a, b) => a * b, 1);
+    if (size !== tensorData.length) {
+        throw new Error(`tensor data has ${tensorData.length} values but dims ` +
+                        `[${d.join(', ')}] need ${size}`);
+    }
+    if (factors.length > 0 && d.length !== factors.length + 1) {
+        throw new Error(`tensor dims need one axis per factor (${factors.length}) ` +
+                        `plus the weight/window axis; got ${d.length}`);
+    }
+    factors.forEach((col, f) => {
+        const levels = factorLevels[col];
+        if (!levels) throw new Error(`no factor levels for tensor factor "${col}"`);
+        for (const [value, level] of Object.entries(levels)) {
+            if (!Number.isInteger(level) || level < 0 || level >= d[f]) {
+                throw new Error(`level ${level} of "${col}" = ${JSON.stringify(value)} ` +
+                                `is outside its axis (size ${d[f]})`);
+            }
+        }
+    });
+}
+
 // libqe's apply_tensor_unit reads the IS_DEFAULT window as an int, so an
 // infinite window (or one too large for int32) is clamped to INT32_MAX — still
 // longer than any conversation, in lines or in seconds.
@@ -354,6 +397,7 @@ export function accumulateTensor(qe, rows, codeMatrix, nRows, nCodes, nUnits,
 
     // Resolve factor levels (infer if not provided)
     const factorLevels = tensorDef.factorLevels ?? inferFactorLevels(rows, factors);
+    validateTensor(dims, tensorData, factors, factorLevels);
     const allTimes = timesCol ? parseTimes(rows, timesCol, timesEndCol || null) : null;
 
     const nConnections = ordered