qi-hardware/cae-tools
1
0
Fork 0
mirror of git://projects.qi-hardware.com/cae-tools.git synced 2024-12-23 03:59:00 +02:00
Code

cameo/stl.c: generate real STL (WIP)

Browse Source 
Some facets look funny in Meshlab although it shows all the normals pointing
in the right direction.
This commit is contained in:
Werner Almesberger 2013-10-14 12:32:27 -03:00
parent ebb7d028ea
commit c5f54df4a0
1 changed files with 68 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

75
cameo/stl.c
View File

@ -11,7 +11,9 @@
*/ */
#include <stdlib.h>
#include <stdio.h> #include <stdio.h>
#include <math.h>
#include "poly2d.h" #include "poly2d.h"
@ -19,20 +21,79 @@
#include "stl.h" #include "stl.h"
void stl(const char *name, const struct path *paths) static void facet(FILE *file, const struct f2d *f, int a, int b, int c,
double za, double zb, double zc, double nx, double ny, double nz)
{
double d;
d = sqrt(nx*nx+ny*ny+nz*nz);
if (!d)
d = 1;
fprintf(file, "facet normal %e %e %e\n", nx/d, ny/d, nz/d);
fprintf(file, "\touter loop\n");
fprintf(file, "\t\tvertex %e %e %e\n", f->x[a], f->y[a], za);
fprintf(file, "\t\tvertex %e %e %e\n", f->x[b], f->y[b], zb);
fprintf(file, "\t\tvertex %e %e %e\n", f->x[c], f->y[c], zc);
fprintf(file, "\tendloop\n");
fprintf(file, "endfacet\n");
}
static void surface(FILE *file, const struct f2d *f)
{
facet(file, f, 0, 1, 2, 1, 1, 1, 0, 0, 1);
facet(file, f, 2, 1, 0, 0, 0, 0, 0, 0, -1);
}
static void side(FILE *file, const struct f2d *f, int a, int b)
{
double nx, ny;
nx = f->y[b]-f->y[a];
ny = f->x[a]-f->x[b];
facet(file, f, a, b, b, 0, 0, 1, nx, ny, 0);
facet(file, f, b, a, a, 1, 1, 0, nx, ny, 0);
}
static void stl_write(FILE *file, const struct path *paths)
{ {
struct p2d *polys; struct p2d *polys;
struct f2d *faces; struct f2d *faces;
const struct f2d *f; const struct f2d *f;
int i;
fprintf(file, "solid cameo\n");
polys = paths_to_polys(paths); polys = paths_to_polys(paths);
faces = f2d_tri(polys); faces = f2d_tri(polys);
for (f = faces; f; f = f->next) for (f = faces; f; f = f->next) {
printf("%f/%f %f/%f %f/%f (%d %d %d)\n", surface(file, f);
f->x[0], f->y[0], for (i = 0; i != 3; i++)
f->x[1], f->y[1], if (f->side[i])
f->x[2], f->y[2], side(file, f, i, (i+1) % 3);
f->side[0], f->side[1], f->side[2]); }
p2d_free_all(polys); p2d_free_all(polys);
f2d_free_all(faces); f2d_free_all(faces);
fprintf(file, "endsolid cameo\n");
}
void stl(const char *name, const struct path *paths)
{
FILE *file = stdout;
if (name) {
file = fopen(name, "w");
if (!file) {
perror(name);
exit(1);
}
}
stl_write(file, paths);
if (name && fclose(file)) {
perror(name);
exit(1);
}
} }