Commit e8fc10c4 authored by Jonathan Lambrechts's avatar Jonathan Lambrechts
Browse files

pure c

parent b2350837
Pipeline #9852 failed with stages
in 2 minutes and 43 seconds
......@@ -169,7 +169,7 @@ def _mesh_bgrid(domain: _geometry.Domain, mesh_size: _geometry.MeshSizeCallback,
data = np.c_[xtri.swapaxes(1,2).reshape(-1,9),v.flatten()[tri]]
view = gmsh.view.add("tri")
gmsh.view.add_list_data(view, "ST", tri.shape[0], data.flatten())
gmsh.fltk.run()
#gmsh.fltk.run()
field = gmsh.model.mesh.field.add("PostView")
gmsh.model.mesh.field.setNumber(field, "ViewTag",view)
gmsh.model.mesh.field.setAsBackgroundMesh(field)
......
......@@ -19,10 +19,10 @@
# see <http://www.gnu.org/licenses/>.
cmake_minimum_required(VERSION 3.9)
project(seamsh C CXX)
project(seamsh C)
add_library(seamsh SHARED
seamsh.c
polymesh.cc
polymesh.c
sort.c
robustPredicates.cpp
robustPredicates.c
)
......@@ -123,8 +123,6 @@
#include <fpu_control.h>
#endif /* LINUX */
namespace robustPredicates
{
/* On some machines, the exact arithmetic routines might be defeated by the */
/* use of internal extended precision floating-point registers. Sometimes */
......@@ -4845,5 +4843,3 @@ REAL orient4d(REAL* pa, REAL* pb, REAL* pc, REAL* pd, REAL* pe,
return orient4dadapt(pa, pb, pc, pd, pe,
aheight, bheight, cheight, dheight, eheight, permanent);
}
} // end namespace
......@@ -7,12 +7,10 @@
#define ROBUST_PREDICATES_H
// namespace necessary to avoid conflicts with predicates used by Tetgen
namespace robustPredicates {
double exactinit(int filter, double maxx, double maxy, double maxz);
double incircle(double *pa, double *pb, double *pc, double *pd);
double insphere(double *pa, double *pb, double *pc, double *pd, double *pe);
double orient2d(double *pa, double *pb, double *pc);
double orient3d(double *pa, double *pb, double *pc, double *pd);
} // namespace robustPredicates
double exactinit(int filter, double maxx, double maxy, double maxz);
double incircle(double *pa, double *pb, double *pc, double *pd);
double insphere(double *pa, double *pb, double *pc, double *pd, double *pe);
double orient2d(double *pa, double *pb, double *pc);
double orient3d(double *pa, double *pb, double *pc, double *pd);
#endif
......@@ -93,10 +93,8 @@ static void _vectorRemoveFlag(void *m, const int *flag, int size) {
}
_vectorPop(m,size*r);
}
extern "C" {
void quicksort (void *const pbase, size_t total_elems, size_t size,
int (*cmp)(const void*, const void*, void*), void *arg);
}
#define vector_remove_flag(v,f,repeat) _vectorRemoveFlag((void*)v,f,repeat*sizeof(*v))
#define vector_size(v) (_vectorSize((void*)v)/sizeof(*v))
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment