Re: Uploading data
OpenLB – Open Source Lattice Boltzmann Code › Forums › on OpenLB › General Topics › Uploading data › Re: Uploading data
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?