/* * Ray2mesh : software for geophysicists. * Compute various scores attached to the mesh cells, based on geometric   information that rays bring when the traverse the cell. * * Copyright (C) 2003, St閜hane Genaud and Marc Grunberg * * This tool is free software; you can redistribute it and/or * modify it under the terms of the GNU Library General Public * License as published by the Free Software Foundation; either * version 2 of the License, or (at your option) any later version. * * This library 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 * Library General Public License for more details. * * You should have received a copy of the GNU Library General Public * License along with this library; if not, write to the Free * Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. */#include <assert.h>#include <math.h>#include <time.h>#include <ray/raydescartes.h>#include <mesh/mesh.h>#include <mesh/layer.h>#include <mesh/cellinfo.h>#include <mesh/convert_coord.h>#include "util.h"#include "ray2mesh.h"#include "raymaths.h"		/* geo_distance */#include "score.h"#include "const_ray2mesh.h"#include "parameterization.h"/*#define DEBUG*//*#define INTERPOLATED_NODE_DEBUG */#define SHOW_CURRENT_NODE fprintf(stderr, "%s:ray2mesh():%d ray #%ld, node #%d, point %f,%f,%f\n", \				__FILE__, __LINE__, \				rayid, nbnode, x0, y0, z0); #define INTERSECT_BUG	  fprintf(stderr, "%s:ray2mesh():%d ray #%ld, intersection bug between (%f