IBWriteBodySTL.c File Reference

Writes a 3D body surface to STL file. More...

#include <stdio.h>
#include <stdlib.h>
#include <basic.h>
#include <mpivars.h>
#include <immersedboundaries.h>

Go to the source code of this file.

Functions
int	IBWriteBodySTL (Body3D *body, char *filename, void *m, int rank, int *stat)

Detailed Description

Writes a 3D body surface to STL file.

Author

Debojyoti Ghosh

Definition in file IBWriteBodySTL.c.

Function Documentation

IBWriteBodySTL()

int IBWriteBodySTL	(	Body3D *	body,
		char *	filename,
		void *	m,
		int	rank,
		int *	stat
	)

Function to write a 3D surface from a STL file. See https://en.wikipedia.org/wiki/STL_(file_format) for details of the STL file format. Notes:

• This function only writes ASCII STL files.

It is assumed that all MPI ranks have the body data. The MPI rank specified as the input rank will actually write the file.

Parameters

body	3D body to write to file
filename	Filename
m	MPI object of type MPIVariables
rank	MPI rank that does the writing
stat	Status (0: success; non-0: failure)

Definition at line 23 of file IBWriteBodySTL.c.

{
  MPIVariables  *mpi = (MPIVariables*) m;

  if (mpi->rank == rank) {
    FILE *out;
    out = fopen(filename,"w");
    if (!out) {
      fprintf(stderr,"Error in IBWriteBodySTL(): Unable to open file ");
      fprintf(stderr,"%s for writing on rank %d.\n",filename,mpi->rank);
      *stat = 1;
    } else {
      fprintf(out,"solid TEST\n");
      int n;
      for (n = 0; n < body->nfacets; n++) {
        fprintf(out,"  facet normal %1.16e %1.16e %1.16e\n",
                body->surface[n].nx,body->surface[n].ny,body->surface[n].nz);
        fprintf(out,"    outer loop\n");
        fprintf(out,"      vertex %1.16e %1.16e %1.16e\n",
                body->surface[n].x1,body->surface[n].y1,body->surface[n].z1);
        fprintf(out,"      vertex %1.16e %1.16e %1.16e\n",
                body->surface[n].x2,body->surface[n].y2,body->surface[n].z2);
        fprintf(out,"      vertex %1.16e %1.16e %1.16e\n",
                body->surface[n].x3,body->surface[n].y3,body->surface[n].z3);
        fprintf(out,"    endloop\n");
        fprintf(out,"  endfacet\n");
      }
      fprintf(out,"endsolid TEST\n");
      *stat = 0;
    }
  }
  return(0);
}