GitLab

This issue was created automatically from an original Mantis Issue. Further discussion may take place here.

---

Dear,

There seems to be a bug in vtkDataSetAttributes when using SetActiveAttribute (int index, int attributeType).

The following works fine (using a name): for(i = 0; i < vec_DataArray.size(); ++i) { sprintf(name, "PCTEST%i", i);

pPolyData->GetPointData()->SetActiveAttribute(name, vtkDataSetAttributes::VECTORS);

double* d = pPolyData->GetPointData()->GetVectors()->GetTuple3(1);

}

when using a number instead of a name:

for(i = 0; i < vec_DataArray.size(); ++i) { sprintf(name, "PCTEST%i", i);

pPolyData->GetPointData()->SetActiveAttribute(i, vtkDataSetAttributes::VECTORS);

double* d = pPolyData->GetPointData()->GetVectors()->GetTuple3(1);

}

it doesn’t change the active vectorfield.

Best wishes, Per

PS: Complete c++ test code:

// testAddArray.cpp : Defines the entry point for the console application. //

#include <iostream> #include <tchar.h>

#include <vtkPolyData.h> #include <vtkDataSetAttributes.h> #include <vtkAssignAttribute.h> #include <vtkPolyDataReader.h> #include <vtkPolyDataWriter.h> #include <vtkDoubleArray.h> #include <vtkPointData.h>

#pragma comment(lib, "vtkHybrid.lib") #pragma comment(lib, "vtkIO.lib") #pragma comment(lib, "vtkGraphics.lib") #pragma comment(lib, "vtkFiltering.lib") #pragma comment(lib, "vtkCommon.lib")

using namespace std;

#include <vector> void main(int argc, _TCHAR* argv[]) {

vtkPolyDataReader *pReader = vtkPolyDataReader::New(); vtkPolyData *pPolyData = vtkPolyData::New(); vtkPolyDataWriter pWriter = vtkPolyDataWriter::New(); vector<vtkDataArray > vec_DataArray; int i; char name[255];

pReader->SetFileName("fran_cut.vtk"); // 'fran_cut.vtk' from the vtk test set pReader->Update(); pPolyData->DeepCopy(pReader->GetOutput());

cerr << "pPolyData->GetNumberOfPoints() = " << pPolyData->GetNumberOfPoints() << "\n";

// Make vector field vec_DataArray.resize(3); for(i = 0; i < vec_DataArray.size(); ++i) { vec_DataArray[i] = vtkDoubleArray::New();

vtkDataArray* v = vec_DataArray[i];

sprintf(name, "PCTEST%i", i);
v->SetName(name);
v->SetNumberOfComponents(3);
v->SetNumberOfTuples(pPolyData->GetNumberOfPoints());

for(int j = 0; j < v->GetNumberOfTuples(); ++j)
{
  const double d = i + j / 1.e4;
  v->SetTuple3(j, d, d, d);
}

}

// Add vectors to dataset // // 11111111111111111 // cerr << "..........1st test...........\n";

cerr << " Shows that many field sets may be added to a dataset.\n";

pPolyData->DeepCopy(pReader->GetOutput());

for(i = 0; i < vec_DataArray.size(); ++i) { pPolyData->GetPointData()->AddArray(vec_DataArray[i]); }

pWriter->SetInput(pPolyData); pWriter->SetFileName("fran_cut_with_test_vectors_1.vtk"); pWriter->Write(); // Write data set with 'FieldData'

// // 22222222222222222 // cerr << "..........2nd test...........\n";

cerr << " Shows that it's not possible to add more the one vector-attribute to a dataset.\n";

pPolyData->DeepCopy(pReader->GetOutput());

for(i = 0; i < vec_DataArray.size(); ++i) { pPolyData->GetPointData()->SetVectors(vec_DataArray[i]); }

// check that all vectors are in dataset! for(i = 0; i < vec_DataArray.size(); ++i) { sprintf(name, "PCTEST%i", i);

pPolyData->GetPointData()->SetActiveVectors(name);

if (pPolyData->GetPointData()->GetVectors()==NULL)
{
  cerr << "no vector with name '" << name << "'....continue\n";
  continue;
}

double* d = pPolyData->GetPointData()->GetVectors()->GetTuple3(1);
cerr << name << " " 
  << d[0] << " "
  << d[1] << " "
  << d[2] << " "
  << "\n";

}

pWriter->SetInput(pPolyData); pWriter->SetFileName("fran_cut_with_test_vectors_2.vtk"); pWriter->Write(); // Write data set with 'FieldData'

// // 33333333333333333 // cerr << "..........3rd test...........\n";

cerr << " Shows how easy field data may be converted/assigned as vector-attribute to the dataset.\n";

pReader->SetFileName("fran_cut_with_test_vectors_1.vtk"); pReader->Update();

pPolyData->DeepCopy(pReader->GetOutput());

cerr << (pPolyData->GetPointData()->GetVectors() != NULL) << " - should be zero as not vector field is assigned/present\n";

for(i = 0; i < vec_DataArray.size(); ++i) { sprintf(name, "PCTEST%i", i);

pPolyData->GetPointData()->SetActiveAttribute(name, vtkDataSetAttributes::VECTORS);

double* d = pPolyData->GetPointData()->GetVectors()->GetTuple3(1);
cerr << name << " " 
  << d[0] << " "
  << d[1] << " "
  << d[2] << " "
  << "\n";

cerr << "just to show that the number also can be used. BUT IT DOESN'T WORK\n";
pPolyData->GetPointData()->SetActiveAttribute((i+1)%vec_DataArray.size(), vtkDataSetAttributes::VECTORS); // just to show that the number also can be used!

d = pPolyData->GetPointData()->GetVectors()->GetTuple3(1);
cerr << (i+1)%vec_DataArray.size() << " " 
  << d[0] << " "
  << d[1] << " "
  << d[2] << " "
  << "\n\n";

}

// // 44444444444444444 // cerr << "..........4th test...........\n";

cerr << " Shows the same as 3rd test but uses vtkAssignAttribute.\n";

vtkAssignAttribute* aa = vtkAssignAttribute::New(); aa->SetInput(pReader->GetOutput());

for(i = 0; i < vec_DataArray.size(); ++i) { sprintf(name, "PCTEST%i", i);

aa->Assign(name, vtkDataSetAttributes::VECTORS,
  vtkAssignAttribute::POINT_DATA); // POINT_DATA: Attribute found in POINT_DATA not cell data
aa->Update();

try
{
  double* d = aa->GetOutput()->GetPointData()->GetVectors()->GetTuple3(1);
  cerr << name << " " 
    << d[0] << " "
    << d[1] << " "
    << d[2] << " "
    << "\n";
}
catch (...) {cerr << "exception!!!\n";}

}

#define VTKRELEASE(x) x->Delete()

VTKRELEASE(pReader); VTKRELEASE(pPolyData); VTKRELEASE(pWriter);

for(i = 0; i < vec_DataArray.size(); ++i) VTKRELEASE(vec_DataArray[i]);

return;

}