void Fluid2DParticlesApp::touchesMoved( TouchEvent event ) { float s = 10; const std::vector<TouchEvent::Touch>& touches = event.getTouches(); for( std::vector<TouchEvent::Touch>::const_iterator cit = touches.begin(); cit != touches.end(); ++cit ) { if( mTouchColors.find( cit->getId() ) == mTouchColors.end() ) continue; vec2 prevPos = cit->getPrevPos(); vec2 pos = cit->getPos(); float x = (pos.x/(float)getWindowWidth())*mFluid2D.resX(); float y = (pos.y/(float)getWindowHeight())*mFluid2D.resY(); vec2 dv = pos - prevPos; mFluid2D.splatVelocity( x, y, mVelScale*dv ); mFluid2D.splatRgb( x, y, mRgbScale*mTouchColors[cit->getId()] ); if( mFluid2D.isBuoyancyEnabled() ) { mFluid2D.splatDensity( x, y, mDenScale ); } for( int i = 0; i < 5; ++i ) { vec2 partPos = pos + vec2( Rand::randFloat( -s, s ), Rand::randFloat( -s, s ) ); float life = Rand::randFloat( 3.0f, 6.0f ); mParticles.append( Particle( partPos, life, mTouchColors[cit->getId()] ) ); } } }
void Fluid2DCamAppApp::mouseDrag( MouseEvent event ) { float x = (event.getX()/(float)getWindowWidth())*mFluid2D.resX(); float y = (event.getY()/(float)getWindowHeight())*mFluid2D.resY(); if( event.isLeftDown() ) { Vec2f dv = event.getPos() - mPrevPos; mFluid2D.splatVelocity( x, y, mVelScale*dv ); mFluid2D.splatDensity( x, y, mDenScale ); } mPrevPos = event.getPos(); }
void Fluid2DBasicApp::touchesMoved( TouchEvent event ) { const std::vector<TouchEvent::Touch>& touches = event.getTouches(); for( std::vector<TouchEvent::Touch>::const_iterator cit = touches.begin(); cit != touches.end(); ++cit ) { Vec2f prevPos = cit->getPrevPos(); Vec2f pos = cit->getPos(); float x = (pos.x/(float)getWindowWidth())*mFluid2D.resX(); float y = (pos.y/(float)getWindowHeight())*mFluid2D.resY(); Vec2f dv = pos - prevPos; mFluid2D.splatVelocity( x, y, mVelScale*dv ); mFluid2D.splatDensity( x, y, mDenScale ); } }
void Fluid2DTextureApp::mouseDrag( MouseEvent event ) { float x = (event.getX()/(float)getWindowWidth())*mFluid2D.resX(); float y = (event.getY()/(float)getWindowHeight())*mFluid2D.resY(); if( event.isLeftDown() ) { vec2 dv = vec2( event.getPos() ) - mPrevPos; mFluid2D.splatVelocity( x, y, mVelScale*dv ); if( mFluid2D.isBuoyancyEnabled() ) { mFluid2D.splatDensity( x, y, mDenScale ); } } mPrevPos = event.getPos(); }
void Fluid2DParticleSoupApp::mouseDrag( MouseEvent event ) { float x = (event.getX()/(float)getWindowWidth())*mFluid2D.resX(); float y = (event.getY()/(float)getWindowHeight())*mFluid2D.resY(); if( event.isLeftDown() ) { Vec2f dv = event.getPos() - mPrevPos; mFluid2D.splatVelocity( x, y, mVelScale*dv ); mFluid2D.splatRgb( x, y, mRgbScale*mColor ); if( mFluid2D.isBuoyancyEnabled() ) { mFluid2D.splatDensity( x, y, mDenScale ); } } mPrevPos = event.getPos(); }
void Fluid2DParticleSoupApp::touchesMoved( TouchEvent event ) { const std::vector<TouchEvent::Touch>& touches = event.getTouches(); for( std::vector<TouchEvent::Touch>::const_iterator cit = touches.begin(); cit != touches.end(); ++cit ) { vec2 prevPos = cit->getPrevPos(); vec2 pos = cit->getPos(); float x = (pos.x/(float)getWindowWidth())*mFluid2D.resX(); float y = (pos.y/(float)getWindowHeight())*mFluid2D.resY(); vec2 dv = pos - prevPos; mFluid2D.splatVelocity( x, y, mVelScale*dv ); mFluid2D.splatRgb( x, y, mRgbScale*mColor ); if( mFluid2D.isBuoyancyEnabled() ) { mFluid2D.splatDensity( x, y, mDenScale ); } } }
void Fluid2DParticlesApp::mouseDrag( MouseEvent event ) { float x = (event.getX()/(float)getWindowWidth())*mFluid2D.resX(); float y = (event.getY()/(float)getWindowHeight())*mFluid2D.resY(); float s = 10; if( event.isLeftDown() ) { vec2 dv = vec2( event.getPos() ) - mPrevPos; mFluid2D.splatVelocity( x, y, mVelScale*dv ); mFluid2D.splatRgb( x, y, mRgbScale*mColor ); if( mFluid2D.isBuoyancyEnabled() ) { mFluid2D.splatDensity( x, y, mDenScale ); } // for( int i = 0; i < 10; ++i ) { vec2 partPos = vec2( event.getPos() ) + vec2( Rand::randFloat( -s, s ), Rand::randFloat( -s, s ) ); float life = Rand::randFloat( 2.0f, 4.0f ); mParticles.append( Particle( partPos, life, mColor ) ); } } mPrevPos = event.getPos(); }
void Fluid2DCamAppApp::update() { if( mCapture && mCapture.checkNewFrame() ) { if( ! mTexCam ) { mTexCam = gl::Texture( mCapture.getSurface() ); } // Flip the image if( ! mFlipped ) { Surface8u srcImg = mCapture.getSurface(); mFlipped = Surface8u( srcImg.getWidth(), srcImg.getHeight(), srcImg.hasAlpha(), srcImg.getChannelOrder() ); } Surface8u srcImg = mCapture.getSurface(); mFlipped = Surface8u( srcImg.getWidth(), srcImg.getHeight(), srcImg.hasAlpha(), srcImg.getChannelOrder() ); for( int y = 0; y < mCapture.getHeight(); ++y ) { const Color8u* src = (const Color8u*)(srcImg.getData() + (y + 1)*srcImg.getRowBytes() - srcImg.getPixelInc()); Color8u* dst = (Color8u*)(mFlipped.getData() + y*mFlipped.getRowBytes()); for( int x = 0; x < mCapture.getWidth(); ++x ) { *dst = *src; ++dst; --src; } } // Create scaled image if( ! mCurScaled ) { mCurScaled = Surface8u( mFlipped.getWidth()/kFlowScale, mFlipped.getHeight()/kFlowScale, mFlipped.hasAlpha(), mFlipped.getChannelOrder() ); } ip::resize( mFlipped, &mCurScaled ); // Optical flow if( mCurScaled && mPrvScaled ) { mPrvCvData = mCurCvData; mCurCvData = cv::Mat( toOcv( Channel( mCurScaled ) ) ); if( mPrvCvData.data && mCurCvData.data ) { int pyrLvels = 3; int winSize = 3; int iters = 5; int poly_n = 7; double poly_sigma = 1.5; cv::calcOpticalFlowFarneback( mPrvCvData, mCurCvData, mFlow, 0.5, pyrLvels, 2*winSize + 1, iters, poly_n, poly_sigma, cv::OPTFLOW_FARNEBACK_GAUSSIAN ); if( mFlow.data ) { if( mFlowVectors.empty() ) { mFlowVectors.resize( mCurScaled.getWidth()*mCurScaled.getHeight() ); } //memset( &mFlowVectors[0], 0, mCurScaled.getWidth()*mCurScaled.getHeight()*sizeof( Vec2f ) ); mNumActiveFlowVectors = 0; for( int j = 0; j < mCurScaled.getHeight(); ++j ) { for( int i = 0; i < mCurScaled.getWidth(); ++i ) { const float* fptr = reinterpret_cast<float*>(mFlow.data + j*mFlow.step + i*sizeof(float)*2); // Vec2f v = Vec2f( fptr[0], fptr[1] ); if( v.lengthSquared() >= mVelThreshold ) { if( mNumActiveFlowVectors >= (int)mFlowVectors.size() ) { mFlowVectors.push_back( std::make_pair( Vec2i( i, j ), v ) ); } else { mFlowVectors[mNumActiveFlowVectors] = std::make_pair( Vec2i( i, j ), v ); } ++mNumActiveFlowVectors; } } } } } } // Update texture mTexCam.update( mFlipped ); // Save previous frame if( ! mPrvScaled ) { mPrvScaled = Surface8u( mCurScaled.getWidth(), mCurScaled.getHeight(), mCurScaled.hasAlpha(), mCurScaled.getChannelOrder() ); } memcpy( mPrvScaled.getData(), mCurScaled.getData(), mCurScaled.getHeight()*mCurScaled.getRowBytes() ); } // Update fluid float dx = (mFluid2DResX - 2)/(float)(640/kFlowScale); float dy = (mFluid2DResY - 2)/(float)(480/kFlowScale); for( int i = 0; i < mNumActiveFlowVectors; ++i ) { Vec2f P = mFlowVectors[i].first; const Vec2f& v = mFlowVectors[i].second; mFluid2D.splatDensity( P.x*dx + 1, P.y*dy + 1, mDenScale*v.lengthSquared() ); mFluid2D.splatVelocity( P.x*dx + 1, P.y*dy + 1, v*mVelScale ); } mFluid2D.step(); // Update velocity const Vec2f* srcVel0 = mFluid2D.dbgVel0().data(); const Vec2f* srcVel1 = mFluid2D.dbgVel1().data(); Colorf* dstVel0 = (Colorf*)mSurfVel0.getData(); Colorf* dstVel1 = (Colorf*)mSurfVel1.getData(); for( int j = 0; j < mFluid2DResY; ++j ) { for( int i = 0; i < mFluid2DResX; ++i ) { *dstVel0 = Colorf( srcVel0->x, srcVel0->y, 0.0f ); *dstVel1 = Colorf( srcVel1->x, srcVel1->y, 0.0f ); ++srcVel0; ++srcVel1; ++dstVel0; ++dstVel1; } } // Update Density mChanDen0 = Channel32f( mFluid2DResX, mFluid2DResY, mFluid2DResX*sizeof(float), 1, mFluid2D.dbgDen0().data() ); mChanDen1 = Channel32f( mFluid2DResX, mFluid2DResY, mFluid2DResX*sizeof(float), 1, mFluid2D.dbgDen1().data() ); mTexDen0.update( mChanDen0 ); mTexDen1.update( mChanDen1 ); // Update velocity textures mTexVel0.update( mSurfVel0 ); mTexVel1.update( mSurfVel1 ); // Update Divergence mChanDiv = Channel32f( mFluid2DResX, mFluid2DResY, mFluid2DResX*sizeof(float), 1, mFluid2D.dbgDivergence().data() ); mTexDiv.update( mChanDiv ); // Update Divergence mChanPrs = Channel32f( mFluid2DResX, mFluid2DResY, mFluid2DResX*sizeof(float), 1, mFluid2D.dbgPressure().data() ); mTexPrs.update( mChanPrs ); // Update Curl, Curl Length mChanCurl = Channel32f( mFluid2DResX, mFluid2DResY, mFluid2DResX*sizeof(float), 1, mFluid2D.dbgCurl().data() ); mTexCurl.update( mChanCurl ); mChanCurlLen = Channel32f( mFluid2DResX, mFluid2DResY, mFluid2DResX*sizeof(float), 1, mFluid2D.dbgCurlLength().data() ); mTexCurlLen.update( mChanCurlLen ); }