All files / Sources/Filters/Sources/RTAnalyticSource index.js

98.18% Statements 54/55
33.33% Branches 1/3
100% Functions 5/5
98.07% Lines 51/52

Press n or j to go to the next uncovered block, b, p or k for the previous block.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144                22x   22x 21x       21x 21x                   21x             105x     21x   21x 21x 21x 21x   21x 21x 63x     21x   21x   21x 21x 21x 21x           21x 21x 1232x 1232x 1232x 1232x 211451x   211451x 211451x 211451x 40774709x 40774709x   40774709x 40774709x           40774709x         21x       21x   21x 21x     21x               1x                           22x     22x   22x   22x   22x   22x   22x 22x         1x          
import macro from 'vtk.js/Sources/macros';
import vtkImageData from 'vtk.js/Sources/Common/DataModel/ImageData';
import vtkDataArray from 'vtk.js/Sources/Common/Core/DataArray';
 
// ----------------------------------------------------------------------------
 
function vtkRTAnalyticSource(publicAPI, model) {
  // Set our className
  model.classHierarchy.push('vtkRTAnalyticSource');
 
  publicAPI.requestData = (inData, outData) => {
    Iif (model.deleted) {
      return;
    }
 
    const state = {};
    const dataset = {
      type: 'vtkImageData',
      mtime: model.mtime,
      metadata: {
        source: 'vtkRTAnalyticSource',
        state,
      },
    };
 
    // Add parameter used to create dataset as metadata.state[*]
    [
      'standardDeviation',
      'center',
      'frequency',
      'magnitude',
      'maximum',
    ].forEach((field) => {
      state[field] = [].concat(model[field]);
    });
 
    const id = vtkImageData.newInstance(dataset);
 
    id.setSpacing(1.0, 1.0, 1.0);
    id.setExtent.apply(this, model.wholeExtent);
    id.setOrigin(0.0, 0.0, 0.0);
    id.setDirection(model.dataDirection);
 
    let dims = [0, 0, 0];
    dims = dims.map(
      (_, i) => model.wholeExtent[i * 2 + 1] - model.wholeExtent[i * 2] + 1
    );
 
    const newArray = new Uint8Array(dims[0] * dims[1] * dims[2]);
    const temp2 =
      1.0 / (2.0 * model.standardDeviation * model.standardDeviation);
 
    let xval = 0;
    let yval = 0;
    let zval = 0;
    const scale = [
      1.0 / (model.wholeExtent[1] - model.wholeExtent[0]),
      1.0 / (model.wholeExtent[3] - model.wholeExtent[2]),
      1.0 / (model.wholeExtent[5] - model.wholeExtent[4]),
    ];
 
    let i = 0;
    for (let z = model.wholeExtent[4]; z <= model.wholeExtent[5]; z++) {
      zval = (model.center[2] - z) * scale[2];
      const zfactor = model.magnitude[2] * Math.cos(model.frequency[2] * zval);
      zval *= zval;
      for (let y = model.wholeExtent[2]; y <= model.wholeExtent[3]; y++) {
        yval = (model.center[1] - y) * scale[1];
        const yfactor =
          model.magnitude[1] * Math.sin(model.frequency[1] * yval);
        yval *= yval;
        for (let x = model.wholeExtent[0]; x <= model.wholeExtent[1]; x++) {
          xval = (model.center[0] - x) * scale[0];
          const sum = zval + yval + xval * xval;
          const xfactor =
            model.magnitude[0] * Math.sin(model.frequency[0] * xval);
          newArray[i] =
            model.maximum * Math.exp(-sum * temp2) +
            xfactor +
            yfactor +
            zfactor +
            model.offset;
          i++;
        }
      }
    }
 
    const da = vtkDataArray.newInstance({
      numberOfComponents: 1,
      values: newArray,
    });
    da.setName('scalars');
 
    const cpd = id.getPointData();
    cpd.setScalars(da);
 
    // Update output
    outData[0] = id;
  };
}
 
// ----------------------------------------------------------------------------
// Object factory
// ----------------------------------------------------------------------------
 
const DEFAULT_VALUES = {
  offset: 40,
  maximum: 120,
  center: [0, 0, 0],
  frequency: [60, 30, 40],
  magnitude: [10, 18, 5],
  standardDeviation: 0.5,
  wholeExtent: [-10, 10, -10, 10, -10, 10],
  dataDirection: [1, 0, 0, 0, 1, 0, 0, 0, 1],
};
 
// ----------------------------------------------------------------------------
 
export function extend(publicAPI, model, initialValues = {}) {
  Object.assign(model, DEFAULT_VALUES, initialValues);
 
  // Build VTK API
  macro.obj(publicAPI, model);
 
  macro.setGet(publicAPI, model, ['offset', 'maximum', 'standardDeviation']);
 
  macro.setGetArray(publicAPI, model, ['center', 'frequency', 'magnitude'], 3);
 
  macro.setGetArray(publicAPI, model, ['wholeExtent'], 6);
 
  macro.setGetArray(publicAPI, model, ['dataDirection'], 9);
 
  macro.algo(publicAPI, model, 0, 1);
  vtkRTAnalyticSource(publicAPI, model);
}
 
// ----------------------------------------------------------------------------
 
export const newInstance = macro.newInstance(extend, 'vtkRTAnalyticSource');
 
// ----------------------------------------------------------------------------
 
export default { newInstance, extend };