void VTKStencil:: apply ( FlowField & flowField, int i, int j, int k ){ if (i==1 || j==1 || k==1){return;} int flag_value=flowField.getFlags().getValue(i,j,k); FLOAT velocity[3]; FLOAT pressure; flowField.getPressureAndVelocity(pressure,velocity,i,j,k); if (flag_value & OBSTACLE_SELF){ if(output_flag == 0){ *(this->_outputFile)<<"0.0"<<std::endl; }else{ *(this->_outputFile)<<"0.0 0.0 0.0"<<std::endl; } }else{ if(output_flag == 0){ *(this->_outputFile)<< pressure <<std::endl; }else{ *(this->_outputFile)<< velocity[0] << " " << velocity[1] << " " << velocity[2] <<std::endl; } } return ; }