Refactor FEAScript and Newton-Raphson methods

- Updated the context object in FEAScriptModel to explicitly assign eikonalActivationFlag.
- Simplified the assembly function call in newtonRaphson by removing conditional logic for front propagation, ensuring consistent parameter passing.

Author: nikoscham

dist/feascript.cjs.js

@@ -90,7 +90,7 @@ export class FEAScriptModel {
       const context = {
         meshConfig: this.meshConfig,
         boundaryConditions: this.boundaryConditions,
-        eikonalActivationFlag,
+        eikonalActivationFlag: eikonalActivationFlag,
         solverMethod: this.solverMethod,
         initialSolution,
       };

dist/feascript.cjs.js.map

@@ -55,21 +55,12 @@ export function newtonRaphson(assembleMat, context, maxIterations = 100, toleran
     }
 
     // Compute Jacobian and residual matrices
-    if (assembleMat.name === "assembleFrontPropagationMat") {
-      // Pass an additional viscous parameter for front propagation
-      ({ jacobianMatrix, residualVector, nodesCoordinates } = assembleMat(
-        context.meshConfig,
-        context.boundaryConditions,
-        solutionVector, // The solution vector is required in the case of a non-linear equation
-        context.eikonalActivationFlag
-      ));
-    } else {
-      // Standard call for other assembly functions
-      ({ jacobianMatrix, residualVector, nodesCoordinates } = assembleMat(
-        context.meshConfig,
-        context.boundaryConditions
-      ));
-    }
+    ({ jacobianMatrix, residualVector, nodesCoordinates } = assembleMat(
+      context.meshConfig,
+      context.boundaryConditions,
+      solutionVector, // The solution vector is required in the case of a non-linear equation
+      context.eikonalActivationFlag
+    ));
 
     // Solve the linear system based on the specified solver method
     const linearSystemResult = solveLinearSystem(context.solverMethod, jacobianMatrix, residualVector);