コード例 #1
0
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 ;
}