QPMS
/
qpms
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Malá illustrace effektivity řazení r-ek 

Former-commit-id: 12578f504607542b500ce4287cb6341bdb119056
This commit is contained in:
Marek Nečada 2018-08-26 20:28:46 +00:00
parent e84483b2e8
commit 4d039b6031
1 changed files with 27 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

27
tests/lattice/2d_3lat_rcount.c Normal file
View File

@ -0,0 +1,27 @@
#include <qpms/lattices.h>
#include <stdio.h>
#include <float.h>
void dump_points2d_rordered(const points2d_rordered_t *ps, char *filename) {
FILE *f = fopen(filename, "w");
for (size_t i = 0; i < ps->nrs; ++i) {
fprintf(f, "# r = %.16g\n", ps->rs[i]);
for (ptrdiff_t j = ps->r_offsets[i]; j < ps->r_offsets[i+1]; ++j)
fprintf(f, "%.16g %.16g\n", ps->base[j].x, ps->base[j].y);
}
fclose(f);
}
int main() {
triangular_lattice_gen_t *g = triangular_lattice_gen_init(1, TRIANGULAR_HORIZONTAL, false,0);
triangular_lattice_gen_extend_to_steps(g, 1000);
for(size_t i = 0; i < g->ps.nrs; ++i) {
printf("%zd %.16g %td\n", i, g->ps.rs[i], g->ps.r_offsets[i+1]);
}
triangular_lattice_gen_free(g);
return 0;
}