danielherrerohernando

import clr

clr.AddReference('ProtoGeometry')
from Autodesk.DesignScript.Geometry import *

clr.AddReference('RevitAPI')
from Autodesk.Revit.DB import *

clr.AddReference('RevitNodes')
import Revit

danielherrerohernando

import clr

clr.AddReference('ProtoGeometry')
from Autodesk.DesignScript.Geometry import *

clr.AddReference('RevitAPI')
from Autodesk.Revit.DB import *

clr.AddReference('RevitNodes')
import Revit

danielherrerohernando

const { writeFileSync, readFileSync } = require('fs');
const { buildSystem } = require('./systemBuilder');
const { solve } = require('./systemSolver');
const { getExpression } = require('./helpers');

const createModel = () => {

  let params = {};

  const fit = (data,degrees) => {

danielherrerohernando

const buildAugmentedMatrix = (leftMatrix, rightMatrix) => leftMatrix.map((row,i)=>row.concat(rightMatrix[i]));

const triangularize = (augmentedMatrix) => {
  const n = augmentedMatrix.length;
  for (let i=0; i<n-1; i++) {
    for (let j=i+1; j<n; j++) {
      const c = augmentedMatrix[j][i]/augmentedMatrix[i][i];
      for (let k=i+1; k<n+1; k++) {
        augmentedMatrix[j][k] = augmentedMatrix[j][k] - c*augmentedMatrix[i][k];
      }

danielherrerohernando

const buildArray = size => new Array(size).fill('');

const buildBase = degree => buildArray(degree+1).map((_,i) => x => x**i);

const buildGramMatrix = (data, base) => buildArray(base.length)
  .map((_,i) => buildArray(base.length)
  .map((_,j) => data.reduce((acc,[x]) => acc + base[i](x)*base[j](x),0)));

const buildIndependentTerm = (data, base) => buildArray(base.length)
  .map((_,i) => [].concat(data.reduce((acc,[x,y]) => acc + base[i](x)*y,0)));