コード例 #1
0
void DiffuseParticleSimulation::
        addDiffuseParticles(FragmentedVector<DiffuseParticle> &particles) {
    _diffuseParticles.reserve(_diffuseParticles.size() + particles.size());
    for (unsigned int i = 0; i < particles.size(); i++) {
        _diffuseParticles.push_back(particles[i]);
    }
}
コード例 #2
0
void IsotropicParticleMesher::_addPointsToScalarFieldAccelerator(FragmentedVector<MarkerParticle> &points,
                                                                 ImplicitSurfaceScalarField &field) {
	bool isThresholdSet = _scalarFieldAccelerator->isMaxScalarFieldValueThresholdSet();
	bool origThreshold = _scalarFieldAccelerator->getMaxScalarFieldValueThreshold();
	_scalarFieldAccelerator->setMaxScalarFieldValueThreshold(_maxScalarFieldValueThreshold);

	int n = _maxParticlesPerScalarFieldAddition;
	std::vector<vmath::vec3> positions;
    positions.reserve(fmin(n, points.size()));

    for (int startidx = 0; startidx < (int)points.size(); startidx += n) {
        int endidx = startidx + n - 1;
        if (endidx >= (int)points.size()) {
            endidx = points.size() - 1;
        }

        positions.clear();
        for (int i = startidx; i <= endidx; i++) {
            positions.push_back(points[i].position);
        }

        _scalarFieldAccelerator->addPoints(positions, field);
    }

    if (!isThresholdSet) {
    	_scalarFieldAccelerator->setMaxScalarFieldValueThreshold();
    } else {
    	_scalarFieldAccelerator->setMaxScalarFieldValueThreshold(origThreshold);
    }
}
コード例 #3
0
void DiffuseParticleSimulation::
        setDiffuseParticles(FragmentedVector<DiffuseParticle> &particles) {
    _diffuseParticles.clear();
    _diffuseParticles.shrink_to_fit();
    _diffuseParticles.reserve(particles.size());
    for (unsigned int i = 0; i < particles.size(); i++) {
        _diffuseParticles.push_back(particles[i]);
    }
}
コード例 #4
0
TriangleMesh IsotropicParticleMesher::_polygonizeAll(FragmentedVector<MarkerParticle> &particles, 
	                                                 FluidMaterialGrid &materialGrid) {
	int subd = _subdivisionLevel;
	int width = _isize*subd;
	int height = _jsize*subd;
	int depth = _ksize*subd;
	double dx = _dx / (double)subd;

	ImplicitSurfaceScalarField field(width + 1, height + 1, depth + 1, dx);

	int origsubd = materialGrid.getSubdivisionLevel();
	materialGrid.setSubdivisionLevel(subd);
	field.setMaterialGrid(materialGrid);
	materialGrid.setSubdivisionLevel(origsubd);

	field.setPointRadius(_particleRadius);

    for (unsigned int i = 0; i < particles.size(); i++) {
        field.addPoint(particles[i].position);
    }

    Polygonizer3d polygonizer(&field);
    polygonizer.polygonizeSurface();

    return polygonizer.getTriangleMesh();
}
コード例 #5
0
void IsotropicParticleMesher::_computeSliceScalarField(int startidx, int endidx, 
	                                                   FragmentedVector<MarkerParticle> &markerParticles,
	                                                   FluidMaterialGrid &materialGrid,
	                                                   ImplicitSurfaceScalarField &field) {
	
	FragmentedVector<vmath::vec3> sliceParticles;
	_getSliceParticles(startidx, endidx, markerParticles, sliceParticles);

	int width, height, depth;
	field.getGridDimensions(&width, &height, &depth);

	FluidMaterialGrid sliceMaterialGrid(width - 1, height - 1, depth - 1);
	_getSliceMaterialGrid(startidx, endidx, materialGrid, sliceMaterialGrid);

	field.setMaterialGrid(sliceMaterialGrid);

	vmath::vec3 fieldOffset = _getSliceGridPositionOffset(startidx, endidx);
	field.setOffset(fieldOffset);
	field.setPointRadius(_particleRadius);

	for (unsigned int i = 0; i < sliceParticles.size(); i++) {
        field.addPoint(sliceParticles[i]);
    }

    _updateScalarFieldSeam(startidx, endidx, field);
}
コード例 #6
0
void IsotropicParticleMesher::_addPointsToScalarField(FragmentedVector<MarkerParticle> &points,
	                                                  ImplicitSurfaceScalarField &field) {
	if (_isScalarFieldAcceleratorSet) {
		_addPointsToScalarFieldAccelerator(points, field);
	} else {
		for (unsigned int i = 0; i < points.size(); i++) {
	        field.addPoint(points[i].position);
	    }
	}
}
コード例 #7
0
void IsotropicParticleMesher::_getSliceParticles(int startidx, int endidx, 
	                                             FragmentedVector<MarkerParticle> &markerParticles,
		                                         FragmentedVector<vmath::vec3> &sliceParticles) {
	AABB bbox = _getSliceAABB(startidx, endidx);
	for (unsigned int i = 0; i < markerParticles.size(); i++) {
		if (bbox.isPointInside(markerParticles[i].position)) {
			sliceParticles.push_back(markerParticles[i].position);
		}
	}
}