admesh 0.1.0

data/ext/admesh/admesh/src/stl.h

@@ -0,0 +1,201 @@
+/*  ADMesh -- process triangulated solid meshes
+ *  Copyright (C) 1995, 1996  Anthony D. Martin <amartin@engr.csulb.edu>
+ *  Copyright (C) 2013, 2014  several contributors, see AUTHORS
+ *
+ *  This program is free software; you can redistribute it and/or modify
+ *  it under the terms of the GNU General Public License as published by
+ *  the Free Software Foundation; either version 2 of the License, or
+ *  (at your option) any later version.
+
+ *  This program is distributed in the hope that it will be useful,
+ *  but WITHOUT ANY WARRANTY; without even the implied warranty of
+ *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ *  GNU General Public License for more details.
+
+ *  You should have received a copy of the GNU General Public License along
+ *  with this program; if not, write to the Free Software Foundation, Inc.,
+ *  51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
+ *
+ *  Questions, comments, suggestions, etc to
+ *           https://github.com/admesh/admesh/issues
+ */
+
+#ifndef __admesh_stl__
+#define __admesh_stl__
+
+#include <stdio.h>
+
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+#define STL_MAX(A,B) ((A)>(B)? (A):(B))
+#define STL_MIN(A,B) ((A)<(B)? (A):(B))
+#define ABS(X)  ((X) < 0 ? -(X) : (X))
+
+#define LABEL_SIZE             80
+#define NUM_FACET_SIZE         4
+#define HEADER_SIZE            84
+#define STL_MIN_FILE_SIZE      284
+#define ASCII_LINES_PER_FACET  7
+#define SIZEOF_EDGE_SORT       24
+
+typedef struct {
+  float x;
+  float y;
+  float z;
+} stl_vertex;
+
+typedef struct {
+  float x;
+  float y;
+  float z;
+} stl_normal;
+
+typedef char stl_extra[2];
+
+typedef struct {
+  stl_normal normal;
+  stl_vertex vertex[3];
+  stl_extra  extra;
+} stl_facet;
+#define SIZEOF_STL_FACET       50
+
+typedef enum {binary, ascii, inmemory} stl_type;
+
+typedef struct {
+  stl_vertex p1;
+  stl_vertex p2;
+  int        facet_number;
+} stl_edge;
+
+typedef struct stl_hash_edge {
+  unsigned       key[6];
+  int            facet_number;
+  int            which_edge;
+  struct stl_hash_edge  *next;
+} stl_hash_edge;
+
+typedef struct {
+  int   neighbor[3];
+  char  which_vertex_not[3];
+} stl_neighbors;
+
+typedef struct {
+  int   vertex[3];
+} v_indices_struct;
+
+typedef struct {
+  char          header[81];
+  stl_type      type;
+  int           number_of_facets;
+  stl_vertex    max;
+  stl_vertex    min;
+  stl_vertex    size;
+  float         bounding_diameter;
+  float         shortest_edge;
+  float         volume;
+  unsigned      number_of_blocks;
+  int           connected_edges;
+  int           connected_facets_1_edge;
+  int           connected_facets_2_edge;
+  int           connected_facets_3_edge;
+  int           facets_w_1_bad_edge;
+  int           facets_w_2_bad_edge;
+  int           facets_w_3_bad_edge;
+  int           original_num_facets;
+  int           edges_fixed;
+  int           degenerate_facets;
+  int           facets_removed;
+  int           facets_added;
+  int           facets_reversed;
+  int           backwards_edges;
+  int           normals_fixed;
+  int           number_of_parts;
+  int           malloced;
+  int           freed;
+  int           facets_malloced;
+  int           collisions;
+  int           shared_vertices;
+  int           shared_malloced;
+} stl_stats;
+
+typedef struct {
+  FILE          *fp;
+  stl_facet     *facet_start;
+  stl_edge      *edge_start;
+  stl_hash_edge **heads;
+  stl_hash_edge *tail;
+  int           M;
+  stl_neighbors *neighbors_start;
+  v_indices_struct *v_indices;
+  stl_vertex    *v_shared;
+  stl_stats     stats;
+  char          error;
+} stl_file;
+
+
+extern void stl_open(stl_file *stl, char *file);
+extern void stl_close(stl_file *stl);
+extern void stl_stats_out(stl_file *stl, FILE *file, char *input_file);
+extern void stl_print_edges(stl_file *stl, FILE *file);
+extern void stl_print_neighbors(stl_file *stl, char *file);
+extern void stl_put_little_int(FILE *fp, int value_in);
+extern void stl_put_little_float(FILE *fp, float value_in);
+extern void stl_write_ascii(stl_file *stl, const char *file, const char *label);
+extern void stl_write_binary(stl_file *stl, const char *file, const char *label);
+extern void stl_write_binary_block(stl_file *stl, FILE *fp);
+extern void stl_check_facets_exact(stl_file *stl);
+extern void stl_check_facets_nearby(stl_file *stl, float tolerance);
+extern void stl_remove_unconnected_facets(stl_file *stl);
+extern void stl_write_vertex(stl_file *stl, int facet, int vertex);
+extern void stl_write_facet(stl_file *stl, char *label, int facet);
+extern void stl_write_edge(stl_file *stl, char *label, stl_hash_edge edge);
+extern void stl_write_neighbor(stl_file *stl, int facet);
+extern void stl_write_quad_object(stl_file *stl, char *file);
+extern void stl_verify_neighbors(stl_file *stl);
+extern void stl_fill_holes(stl_file *stl);
+extern void stl_fix_normal_directions(stl_file *stl);
+extern void stl_fix_normal_values(stl_file *stl);
+extern void stl_reverse_all_facets(stl_file *stl);
+extern void stl_translate(stl_file *stl, float x, float y, float z);
+extern void stl_translate_relative(stl_file *stl, float x, float y, float z);
+extern void stl_scale_versor(stl_file *stl, float versor[3]);
+extern void stl_scale(stl_file *stl, float factor);
+extern void stl_rotate_x(stl_file *stl, float angle);
+extern void stl_rotate_y(stl_file *stl, float angle);
+extern void stl_rotate_z(stl_file *stl, float angle);
+extern void stl_mirror_xy(stl_file *stl);
+extern void stl_mirror_yz(stl_file *stl);
+extern void stl_mirror_xz(stl_file *stl);
+extern void stl_open_merge(stl_file *stl, char *file);
+extern void stl_invalidate_shared_vertices(stl_file *stl);
+extern void stl_generate_shared_vertices(stl_file *stl);
+extern void stl_write_obj(stl_file *stl, char *file);
+extern void stl_write_off(stl_file *stl, char *file);
+extern void stl_write_dxf(stl_file *stl, char *file, char *label);
+extern void stl_write_vrml(stl_file *stl, char *file);
+extern void stl_calculate_normal(float normal[], stl_facet *facet);
+extern void stl_normalize_vector(float v[]);
+extern void stl_calculate_volume(stl_file *stl);
+
+extern void stl_repair(stl_file *stl, int fixall_flag, int exact_flag, int tolerance_flag, float tolerance, int increment_flag, float increment, int nearby_flag, int iterations, int remove_unconnected_flag, int fill_holes_flag, int normal_directions_flag, int normal_values_flag, int reverse_all_flag, int verbose_flag);
+
+extern void stl_initialize(stl_file *stl);
+extern void stl_count_facets(stl_file *stl, char *file);
+extern void stl_allocate(stl_file *stl);
+extern void stl_read(stl_file *stl, int first_facet, int first);
+extern void stl_facet_stats(stl_file *stl, stl_facet facet, int first);
+extern void stl_reallocate(stl_file *stl);
+extern void stl_add_facet(stl_file *stl, stl_facet *new_facet);
+extern void stl_get_size(stl_file *stl);
+
+extern void stl_clear_error(stl_file *stl);
+extern int stl_get_error(stl_file *stl);
+extern void stl_exit_on_error(stl_file *stl);
+
+#ifdef __cplusplus
+}
+#endif
+
+#endif

data/ext/admesh/admesh/src/stl_io.c

@@ -0,0 +1,479 @@
+/*  ADMesh -- process triangulated solid meshes
+ *  Copyright (C) 1995, 1996  Anthony D. Martin <amartin@engr.csulb.edu>
+ *  Copyright (C) 2013, 2014  several contributors, see AUTHORS
+ *
+ *  This program is free software; you can redistribute it and/or modify
+ *  it under the terms of the GNU General Public License as published by
+ *  the Free Software Foundation; either version 2 of the License, or
+ *  (at your option) any later version.
+
+ *  This program is distributed in the hope that it will be useful,
+ *  but WITHOUT ANY WARRANTY; without even the implied warranty of
+ *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ *  GNU General Public License for more details.
+
+ *  You should have received a copy of the GNU General Public License along
+ *  with this program; if not, write to the Free Software Foundation, Inc.,
+ *  51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
+ *
+ *  Questions, comments, suggestions, etc to
+ *           https://github.com/admesh/admesh/issues
+ */
+
+#include <stdlib.h>
+#include <string.h>
+#include "stl.h"
+#include "config.h"
+
+#if !defined(SEEK_SET)
+#define SEEK_SET 0
+#define SEEK_CUR 1
+#define SEEK_END 2
+#endif
+
+void
+stl_print_edges(stl_file *stl, FILE *file) {
+  int i;
+  int edges_allocated;
+
+  if (stl->error) return;
+
+  edges_allocated = stl->stats.number_of_facets * 3;
+  for(i = 0; i < edges_allocated; i++) {
+    fprintf(file, "%d, %f, %f, %f, %f, %f, %f\n",
+            stl->edge_start[i].facet_number,
+            stl->edge_start[i].p1.x, stl->edge_start[i].p1.y,
+            stl->edge_start[i].p1.z, stl->edge_start[i].p2.x,
+            stl->edge_start[i].p2.y, stl->edge_start[i].p2.z);
+  }
+}
+
+
+void
+stl_stats_out(stl_file *stl, FILE *file, char *input_file) {
+  if (stl->error) return;
+
+  /* this is here for Slic3r, without our config.h
+     it won't use this part of the code anyway */
+#ifndef VERSION
+#define VERSION "unknown"
+#endif
+  fprintf(file, "\n\
+================= Results produced by ADMesh version " VERSION " ================\n");
+  fprintf(file, "\
+Input file         : %s\n", input_file);
+  if(stl->stats.type == binary) {
+    fprintf(file, "\
+File type          : Binary STL file\n");
+  } else {
+    fprintf(file, "\
+File type          : ASCII STL file\n");
+  }
+  fprintf(file, "\
+Header             : %s\n", stl->stats.header);
+  fprintf(file, "============== Size ==============\n");
+  fprintf(file, "Min X = % f, Max X = % f\n",
+          stl->stats.min.x, stl->stats.max.x);
+  fprintf(file, "Min Y = % f, Max Y = % f\n",
+          stl->stats.min.y, stl->stats.max.y);
+  fprintf(file, "Min Z = % f, Max Z = % f\n",
+          stl->stats.min.z, stl->stats.max.z);
+
+  fprintf(file, "\
+========= Facet Status ========== Original ============ Final ====\n");
+  fprintf(file, "\
+Number of facets                 : %5d               %5d\n",
+          stl->stats.original_num_facets, stl->stats.number_of_facets);
+  fprintf(file, "\
+Facets with 1 disconnected edge  : %5d               %5d\n",
+          stl->stats.facets_w_1_bad_edge, stl->stats.connected_facets_2_edge -
+          stl->stats.connected_facets_3_edge);
+  fprintf(file, "\
+Facets with 2 disconnected edges : %5d               %5d\n",
+          stl->stats.facets_w_2_bad_edge, stl->stats.connected_facets_1_edge -
+          stl->stats.connected_facets_2_edge);
+  fprintf(file, "\
+Facets with 3 disconnected edges : %5d               %5d\n",
+          stl->stats.facets_w_3_bad_edge, stl->stats.number_of_facets -
+          stl->stats.connected_facets_1_edge);
+  fprintf(file, "\
+Total disconnected facets        : %5d               %5d\n",
+          stl->stats.facets_w_1_bad_edge + stl->stats.facets_w_2_bad_edge +
+          stl->stats.facets_w_3_bad_edge, stl->stats.number_of_facets -
+          stl->stats.connected_facets_3_edge);
+
+  fprintf(file,
+          "=== Processing Statistics ===     ===== Other Statistics =====\n");
+  fprintf(file, "\
+Number of parts       : %5d        Volume   : % f\n",
+          stl->stats.number_of_parts, stl->stats.volume);
+  fprintf(file, "\
+Degenerate facets     : %5d\n", stl->stats.degenerate_facets);
+  fprintf(file, "\
+Edges fixed           : %5d\n", stl->stats.edges_fixed);
+  fprintf(file, "\
+Facets removed        : %5d\n", stl->stats.facets_removed);
+  fprintf(file, "\
+Facets added          : %5d\n", stl->stats.facets_added);
+  fprintf(file, "\
+Facets reversed       : %5d\n", stl->stats.facets_reversed);
+  fprintf(file, "\
+Backwards edges       : %5d\n", stl->stats.backwards_edges);
+  fprintf(file, "\
+Normals fixed         : %5d\n", stl->stats.normals_fixed);
+}
+
+void
+stl_write_ascii(stl_file *stl, const char *file, const char *label) {
+  int       i;
+  FILE      *fp;
+  char      *error_msg;
+
+  if (stl->error) return;
+
+  /* Open the file */
+  fp = fopen(file, "w");
+  if(fp == NULL) {
+    error_msg = (char*)
+                malloc(81 + strlen(file)); /* Allow 80 chars+file size for message */
+    sprintf(error_msg, "stl_write_ascii: Couldn't open %s for writing",
+            file);
+    perror(error_msg);
+    free(error_msg);
+    stl->error = 1;
+    return;
+  }
+
+  fprintf(fp, "solid  %s\n", label);
+
+  for(i = 0; i < stl->stats.number_of_facets; i++) {
+    fprintf(fp, "  facet normal % .8E % .8E % .8E\n",
+            stl->facet_start[i].normal.x, stl->facet_start[i].normal.y,
+            stl->facet_start[i].normal.z);
+    fprintf(fp, "    outer loop\n");
+    fprintf(fp, "      vertex % .8E % .8E % .8E\n",
+            stl->facet_start[i].vertex[0].x, stl->facet_start[i].vertex[0].y,
+            stl->facet_start[i].vertex[0].z);
+    fprintf(fp, "      vertex % .8E % .8E % .8E\n",
+            stl->facet_start[i].vertex[1].x, stl->facet_start[i].vertex[1].y,
+            stl->facet_start[i].vertex[1].z);
+    fprintf(fp, "      vertex % .8E % .8E % .8E\n",
+            stl->facet_start[i].vertex[2].x, stl->facet_start[i].vertex[2].y,
+            stl->facet_start[i].vertex[2].z);
+    fprintf(fp, "    endloop\n");
+    fprintf(fp, "  endfacet\n");
+  }
+
+  fprintf(fp, "endsolid  %s\n", label);
+
+  fclose(fp);
+}
+
+void
+stl_print_neighbors(stl_file *stl, char *file) {
+  int i;
+  FILE *fp;
+  char *error_msg;
+
+  if (stl->error) return;
+
+  /* Open the file */
+  fp = fopen(file, "w");
+  if(fp == NULL) {
+    error_msg = (char*)
+                malloc(81 + strlen(file)); /* Allow 80 chars+file size for message */
+    sprintf(error_msg, "stl_print_neighbors: Couldn't open %s for writing",
+            file);
+    perror(error_msg);
+    free(error_msg);
+    stl->error = 1;
+    return;
+  }
+
+  for(i = 0; i < stl->stats.number_of_facets; i++) {
+    fprintf(fp, "%d, %d,%d, %d,%d, %d,%d\n",
+            i,
+            stl->neighbors_start[i].neighbor[0],
+            (int)stl->neighbors_start[i].which_vertex_not[0],
+            stl->neighbors_start[i].neighbor[1],
+            (int)stl->neighbors_start[i].which_vertex_not[1],
+            stl->neighbors_start[i].neighbor[2],
+            (int)stl->neighbors_start[i].which_vertex_not[2]);
+  }
+  fclose(fp);
+}
+
+void
+stl_put_little_int(FILE *fp, int value_in) {
+  int new_value;
+  union {
+    int  int_value;
+    char char_value[4];
+  } value;
+
+  value.int_value = value_in;
+
+  new_value  = value.char_value[0] & 0xFF;
+  new_value |= (value.char_value[1] & 0xFF) << 0x08;
+  new_value |= (value.char_value[2] & 0xFF) << 0x10;
+  new_value |= (value.char_value[3] & 0xFF) << 0x18;
+  fwrite(&new_value, sizeof(int), 1, fp);
+}
+
+void
+stl_put_little_float(FILE *fp, float value_in) {
+  int new_value;
+  union {
+    float float_value;
+    char  char_value[4];
+  } value;
+
+  value.float_value = value_in;
+
+  new_value  = value.char_value[0] & 0xFF;
+  new_value |= (value.char_value[1] & 0xFF) << 0x08;
+  new_value |= (value.char_value[2] & 0xFF) << 0x10;
+  new_value |= (value.char_value[3] & 0xFF) << 0x18;
+  fwrite(&new_value, sizeof(int), 1, fp);
+}
+
+void
+stl_write_binary_block(stl_file *stl, FILE *fp)
+{
+  int i;
+  for(i = 0; i < stl->stats.number_of_facets; i++)
+    {
+      stl_put_little_float(fp, stl->facet_start[i].normal.x);
+      stl_put_little_float(fp, stl->facet_start[i].normal.y);
+      stl_put_little_float(fp, stl->facet_start[i].normal.z);
+      stl_put_little_float(fp, stl->facet_start[i].vertex[0].x);
+      stl_put_little_float(fp, stl->facet_start[i].vertex[0].y);
+      stl_put_little_float(fp, stl->facet_start[i].vertex[0].z);
+      stl_put_little_float(fp, stl->facet_start[i].vertex[1].x);
+      stl_put_little_float(fp, stl->facet_start[i].vertex[1].y);
+      stl_put_little_float(fp, stl->facet_start[i].vertex[1].z);
+      stl_put_little_float(fp, stl->facet_start[i].vertex[2].x);
+      stl_put_little_float(fp, stl->facet_start[i].vertex[2].y);
+      stl_put_little_float(fp, stl->facet_start[i].vertex[2].z);
+      fputc(stl->facet_start[i].extra[0], fp);
+      fputc(stl->facet_start[i].extra[1], fp);
+    }
+}
+
+void
+stl_write_binary(stl_file *stl, const char *file, const char *label) {
+  FILE      *fp;
+  int       i;
+  char      *error_msg;
+
+  if (stl->error) return;
+
+  /* Open the file */
+  fp = fopen(file, "w");
+  if(fp == NULL) {
+    error_msg = (char*)
+                malloc(81 + strlen(file)); /* Allow 80 chars+file size for message */
+    sprintf(error_msg, "stl_write_binary: Couldn't open %s for writing",
+            file);
+    perror(error_msg);
+    free(error_msg);
+    stl->error = 1;
+    return;
+  }
+
+  fprintf(fp, "%s", label);
+  for(i = strlen(label); i < LABEL_SIZE; i++) putc(0, fp);
+
+  fseek(fp, LABEL_SIZE, SEEK_SET);
+
+  stl_put_little_int(fp, stl->stats.number_of_facets);
+
+  stl_write_binary_block(stl, fp);
+
+  fclose(fp);
+}
+
+void
+stl_write_vertex(stl_file *stl, int facet, int vertex) {
+  if (stl->error) return;
+  printf("  vertex %d/%d % .8E % .8E % .8E\n", vertex, facet,
+         stl->facet_start[facet].vertex[vertex].x,
+         stl->facet_start[facet].vertex[vertex].y,
+         stl->facet_start[facet].vertex[vertex].z);
+}
+
+void
+stl_write_facet(stl_file *stl, char *label, int facet) {
+  if (stl->error) return;
+  printf("facet (%d)/ %s\n", facet, label);
+  stl_write_vertex(stl, facet, 0);
+  stl_write_vertex(stl, facet, 1);
+  stl_write_vertex(stl, facet, 2);
+}
+
+void
+stl_write_edge(stl_file *stl, char *label, stl_hash_edge edge) {
+  if (stl->error) return;
+  printf("edge (%d)/(%d) %s\n", edge.facet_number, edge.which_edge, label);
+  if(edge.which_edge < 3) {
+    stl_write_vertex(stl, edge.facet_number, edge.which_edge % 3);
+    stl_write_vertex(stl, edge.facet_number, (edge.which_edge + 1) % 3);
+  } else {
+    stl_write_vertex(stl, edge.facet_number, (edge.which_edge + 1) % 3);
+    stl_write_vertex(stl, edge.facet_number, edge.which_edge % 3);
+  }
+}
+
+void
+stl_write_neighbor(stl_file *stl, int facet) {
+  if (stl->error) return;
+  printf("Neighbors %d: %d, %d, %d ;  %d, %d, %d\n", facet,
+         stl->neighbors_start[facet].neighbor[0],
+         stl->neighbors_start[facet].neighbor[1],
+         stl->neighbors_start[facet].neighbor[2],
+         stl->neighbors_start[facet].which_vertex_not[0],
+         stl->neighbors_start[facet].which_vertex_not[1],
+         stl->neighbors_start[facet].which_vertex_not[2]);
+}
+
+void
+stl_write_quad_object(stl_file *stl, char *file) {
+  FILE      *fp;
+  int       i;
+  int       j;
+  char      *error_msg;
+  stl_vertex connect_color;
+  stl_vertex uncon_1_color;
+  stl_vertex uncon_2_color;
+  stl_vertex uncon_3_color;
+  stl_vertex color;
+
+  if (stl->error) return;
+
+  /* Open the file */
+  fp = fopen(file, "w");
+  if(fp == NULL) {
+    error_msg = (char*)
+                malloc(81 + strlen(file)); /* Allow 80 chars+file size for message */
+    sprintf(error_msg, "stl_write_quad_object: Couldn't open %s for writing",
+            file);
+    perror(error_msg);
+    free(error_msg);
+    stl->error = 1;
+    return;
+  }
+
+  connect_color.x = 0.0;
+  connect_color.y = 0.0;
+  connect_color.z = 1.0;
+  uncon_1_color.x = 0.0;
+  uncon_1_color.y = 1.0;
+  uncon_1_color.z = 0.0;
+  uncon_2_color.x = 1.0;
+  uncon_2_color.y = 1.0;
+  uncon_2_color.z = 1.0;
+  uncon_3_color.x = 1.0;
+  uncon_3_color.y = 0.0;
+  uncon_3_color.z = 0.0;
+
+  fprintf(fp, "CQUAD\n");
+  for(i = 0; i < stl->stats.number_of_facets; i++) {
+    j = ((stl->neighbors_start[i].neighbor[0] == -1) +
+         (stl->neighbors_start[i].neighbor[1] == -1) +
+         (stl->neighbors_start[i].neighbor[2] == -1));
+    if(j == 0) {
+      color = connect_color;
+    } else if(j == 1) {
+      color = uncon_1_color;
+    } else if(j == 2) {
+      color = uncon_2_color;
+    } else {
+      color = uncon_3_color;
+    }
+    fprintf(fp, "%f %f %f    %1.1f %1.1f %1.1f 1\n",
+            stl->facet_start[i].vertex[0].x,
+            stl->facet_start[i].vertex[0].y,
+            stl->facet_start[i].vertex[0].z, color.x, color.y, color.z);
+    fprintf(fp, "%f %f %f    %1.1f %1.1f %1.1f 1\n",
+            stl->facet_start[i].vertex[1].x,
+            stl->facet_start[i].vertex[1].y,
+            stl->facet_start[i].vertex[1].z, color.x, color.y, color.z);
+    fprintf(fp, "%f %f %f    %1.1f %1.1f %1.1f 1\n",
+            stl->facet_start[i].vertex[2].x,
+            stl->facet_start[i].vertex[2].y,
+            stl->facet_start[i].vertex[2].z, color.x, color.y, color.z);
+    fprintf(fp, "%f %f %f    %1.1f %1.1f %1.1f 1\n",
+            stl->facet_start[i].vertex[2].x,
+            stl->facet_start[i].vertex[2].y,
+            stl->facet_start[i].vertex[2].z, color.x, color.y, color.z);
+  }
+  fclose(fp);
+}
+
+void
+stl_write_dxf(stl_file *stl, char *file, char *label) {
+  int       i;
+  FILE      *fp;
+  char      *error_msg;
+
+  if (stl->error) return;
+
+  /* Open the file */
+  fp = fopen(file, "w");
+  if(fp == NULL) {
+    error_msg = (char*)
+                malloc(81 + strlen(file)); /* Allow 80 chars+file size for message */
+    sprintf(error_msg, "stl_write_ascii: Couldn't open %s for writing",
+            file);
+    perror(error_msg);
+    free(error_msg);
+    stl->error = 1;
+    return;
+  }
+
+  fprintf(fp, "999\n%s\n", label);
+  fprintf(fp, "0\nSECTION\n2\nHEADER\n0\nENDSEC\n");
+  fprintf(fp, "0\nSECTION\n2\nTABLES\n0\nTABLE\n2\nLAYER\n70\n1\n\
+0\nLAYER\n2\n0\n70\n0\n62\n7\n6\nCONTINUOUS\n0\nENDTAB\n0\nENDSEC\n");
+  fprintf(fp, "0\nSECTION\n2\nBLOCKS\n0\nENDSEC\n");
+
+  fprintf(fp, "0\nSECTION\n2\nENTITIES\n");
+
+  for(i = 0; i < stl->stats.number_of_facets; i++) {
+    fprintf(fp, "0\n3DFACE\n8\n0\n");
+    fprintf(fp, "10\n%f\n20\n%f\n30\n%f\n",
+            stl->facet_start[i].vertex[0].x, stl->facet_start[i].vertex[0].y,
+            stl->facet_start[i].vertex[0].z);
+    fprintf(fp, "11\n%f\n21\n%f\n31\n%f\n",
+            stl->facet_start[i].vertex[1].x, stl->facet_start[i].vertex[1].y,
+            stl->facet_start[i].vertex[1].z);
+    fprintf(fp, "12\n%f\n22\n%f\n32\n%f\n",
+            stl->facet_start[i].vertex[2].x, stl->facet_start[i].vertex[2].y,
+            stl->facet_start[i].vertex[2].z);
+    fprintf(fp, "13\n%f\n23\n%f\n33\n%f\n",
+            stl->facet_start[i].vertex[2].x, stl->facet_start[i].vertex[2].y,
+            stl->facet_start[i].vertex[2].z);
+  }
+
+  fprintf(fp, "0\nENDSEC\n0\nEOF\n");
+
+  fclose(fp);
+}
+
+void
+stl_clear_error(stl_file *stl) {
+  stl->error = 0;
+}
+
+void
+stl_exit_on_error(stl_file *stl) {
+  if (!stl->error) return;
+  stl->error = 0;
+  stl_close(stl);
+  exit(1);
+}
+
+int
+stl_get_error(stl_file *stl) {
+  return stl->error;
+}