Method InterpolateTuple

Testing/GenericBridge/vtkBridgeCell.cxx:463–498

------------------------------------------------------------------------------ Description: Interpolate the attribute `a' at local position `pcoords' of the cell into `val'. \pre a_exists: a!=0 \pre a_is_point_centered: a->GetCentering()==vtkPointCentered \pre clamped_point: pcoords[0]>=0 && pcoords[0]<=1 && pcoords[1]>=0 && pcoords[1]<=1 && pcoords[2]>=0 && pcoords[2]<=1 \pre val_exists: val!=0 \

// \pre val_exists: val!=0
// \pre valid_size: sizeof(val)==a->GetNumberOfComponents()
void vtkBridgeCell::InterpolateTuple(vtkGenericAttribute* a, double pcoords[3], double* val)
{
  assert("pre: a_exists" && a != nullptr);
  assert("pre: a_is_point_centered" && a->GetCentering() == vtkPointCentered);
  assert("pre: clamped_point" &&
    (pcoords[0] >= 0 && pcoords[0] <= 1 && pcoords[1] >= 0 && pcoords[1] <= 1 && pcoords[2] >= 0 &&
      pcoords[2] <= 1));
  assert("pre: val_exists" && val != nullptr);

  vtkBridgeAttribute* ba = static_cast<vtkBridgeAttribute*>(a);

  int componentCount = a->GetNumberOfComponents();
  int ptCount = this->GetNumberOfPoints();

  if (a->GetCentering() == vtkPointCentered)
  {
    this->AllocateWeights();
    this->InterpolationFunctions(pcoords, this->Weights);

    memset(val, 0, sizeof(double) * componentCount);
    for (int pt = 0; pt < ptCount; ++pt)
    {
      ba->Data->GetArray(ba->AttributeNumber)
        ->GetTuple(this->Cell->GetPointId(pt), ba->InternalTuple);
      for (int component = 0; component < componentCount; ++component)
      {
        val[component] += ba->InternalTuple[component] * this->Weights[pt];
      }
    }
  }
  else // cell centered
  {
    // not need to interpolate
    ba->Data->GetArray(ba->AttributeNumber)->GetTuple(this->GetId(), val);
  }
}

//------------------------------------------------------------------------------
// Description: