Re: Uploading data

[kolotinsky]

Dear Fabian, following your example I have written code for my case:

//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
// Initialize force

std::vector<T> spacing(3, T());
spacing[0] = 1;
spacing[1] = 1;
spacing[2] = 1;

olb::BlockVTIreader3D< T, BaseType > datareader(“forcefield.vti”,”AppendedData”);
BlockDataF3D<T, BaseType> blockfield(datareader.getBlockData());
AnalyticalF3D<T,T>* field;
field = new SpecialAnalyticalFfromBlockF3D<T,BaseType>(blockfield, datareader.getCuboid(), spacing);
sLattice.defineExternalField( superGeometry, 1,
DESCRIPTOR<T>::ExternalField::forceBeginsAt,
DESCRIPTOR<T>::ExternalField::sizeOfForce, *field );
/////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////

However, during the simulation I am receiving wrong zero or NaN values of Energy and Umax:

////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
[LatticeStatistics] step=190; t=190; uMax=0; avEnergy=nan; avRho=nan
////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////

I have verified my vti file in Paraview but everything is ok.

I also tried to write the data uploaded from my vti file using vtiWriter provided by OpenLB:

///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
olb::VTIwriter3D< T, BaseType > datawriter;
datawriter.writeData(“field”,”AppendedData”,datareader.getBlockData(),datareader.getCuboid());
//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////

Generated force.vti file contains only the following data:

/////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
<?xml version=”1.0″?>
<VTKFile type=”ImageData” version=”0.1″ byte_order=”LittleEndian”>
<ImageData WholeExtent=” 0 49 0 49 0 49 ” Origin=”0 0 0″ Spacing=”1 1 1″>
<Piece Extent=”0 49 0 49 0 49″>
<PointData>
<DataArray type=”Float32″ Name=”AppendedData” NumberOfComponents=”0″>
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////

What is wrong?