void LBNSCommunicator::setVelocityComponent(
int i,
int j,
int k,
int stencilIndex,
int component,
const int *offset,
const int *flip,
double velocity){
int index=index2array(i,j,k,component,stencilIndex);
//open();
// _logComm<<"i_j_k_index:"<<i
// <<" "<<j<<" "<<k
// <<" stencil:"<<stencilIndex
// <<" component:"<<component
// <<" index:"<<index
// <<" value:"<<velocity
// <<" mid:"<<_mid<<std::endl;
LBNSData data;
data.value=velocity;
data.offset[0]=offset[0];
data.offset[1]=offset[1];
data.offset[2]=offset[2];
data.flip[0]=flip[0];
data.flip[1]=flip[1];
data.flip[2]=flip[2];
_velocities[index].push_back(data);
}
FLOAT& ScalarField::getScalar ( int i, int j, int k ) {
return _data [ index2array ( i, j, k ) ];
}
FLOAT* VectorField::getVector ( int i, int j, int k ) {
return &_data[index2array( i, j, k )];
}
int & IntScalarField::getValue ( int i, int j, int k ) {
return _data[ index2array (i, j, k) ];
}
