Difference between revisions of "Positioning of computational nodes"
(→Filling the domain with nodes) |
|||
| Line 5: | Line 5: | ||
The goal is to fill arbitrary domain, which we can ask if position $\b{p}$ is in the domain our outside, with nodes following the given target density function $\rho(\b{p})$. | The goal is to fill arbitrary domain, which we can ask if position $\b{p}$ is in the domain our outside, with nodes following the given target density function $\rho(\b{p})$. | ||
| − | We start with a simple algorithm based on Poisson Disk Sampling that results in a relatively tightly packed distribution of nodes. What we do is: | + | We start with a simple algorithm based on Poisson Disk Sampling (PDS) that results in a relatively tightly packed distribution of nodes. What we do is: |
* Randomly select a node within the domain. | * Randomly select a node within the domain. | ||
* Add N (depending on the dimensionality of the space) on circle/sphere around that node. | * Add N (depending on the dimensionality of the space) on circle/sphere around that node. | ||
| Line 11: | Line 11: | ||
* Each time check if any of new nodes violates minimal distance criterion, i.e. if it is positioned to close to any of the existing nodes. | * Each time check if any of new nodes violates minimal distance criterion, i.e. if it is positioned to close to any of the existing nodes. | ||
| + | Example of PDS filled domain | ||
| + | [[File:2d_poisson_disk_sampling.png|600px]] | ||
| + | |||
| + | <syntaxhighlight lang="c++" line> | ||
| + | PoissonDiskSamplingFill fill_engine; | ||
| + | // setup | ||
| + | fill_engine.iterations(1000000).optim_after(1000).use_MC(true).proximity_relax(0.8); | ||
| + | // create test domain | ||
| + | RectangleDomain<Vec2d> domain({0, 0}, {1.2, 1.5}); | ||
| + | CircleDomain<Vec2d> o1({0.4, 0.4}, 0.25); | ||
| + | domain.subtract(o1); | ||
| + | // execute fill domain with constant density | ||
| + | auto fill_density = [](Vec2d p)->double{ | ||
| + | return (std::pow(p[0]/10, 2) + 0.01); | ||
| + | }; | ||
| + | fill_engine(domain, fill_density); | ||
| + | </syntaxhighlight> | ||
| − | |||
== Relaxation of the nodal distribution == | == Relaxation of the nodal distribution == | ||
Revision as of 20:52, 15 January 2018
Although the meshless methods do not require any topological relations between nodes and even randomly distributed nodes could be used, it is well-known that using regularly distributed nodes leads to more accurate and more stable results. So despite meshless seeming robustness regarding the nodal distribution, a certain effort has to be invested into the positioning of the nodes and following discussion, to some extent, deals with this problem. In a following discussion we want to form set of algorithms that can cover arbitrary domain in n dimensions with nodes that can be further used to solve systems of PDEs. We start with algorithms for filling the domain with nodes. Second set is denoted to improving the nodal distributions and a last set of algorithms deals with the refinement of the nodal distribution. The algorithms can of course be combined.
Filling the domain with nodes
The goal is to fill arbitrary domain, which we can ask if position $\b{p}$ is in the domain our outside, with nodes following the given target density function $\rho(\b{p})$.
We start with a simple algorithm based on Poisson Disk Sampling (PDS) that results in a relatively tightly packed distribution of nodes. What we do is:
- Randomly select a node within the domain.
- Add N (depending on the dimensionality of the space) on circle/sphere around that node.
- Randomly select one of just added nodes and repeat adding.
- Each time check if any of new nodes violates minimal distance criterion, i.e. if it is positioned to close to any of the existing nodes.
Example of PDS filled domain
1 PoissonDiskSamplingFill fill_engine;
2 // setup
3 fill_engine.iterations(1000000).optim_after(1000).use_MC(true).proximity_relax(0.8);
4 // create test domain
5 RectangleDomain<Vec2d> domain({0, 0}, {1.2, 1.5});
6 CircleDomain<Vec2d> o1({0.4, 0.4}, 0.25);
7 domain.subtract(o1);
8 // execute fill domain with constant density
9 auto fill_density = [](Vec2d p)->double{
10 return (std::pow(p[0]/10, 2) + 0.01);
11 };
12 fill_engine(domain, fill_density);
