Commit 8b908e44 authored by Uwe Schulzweida's avatar Uwe Schulzweida
Browse files

clipping update.

parent e0c91feb
......@@ -41,7 +41,7 @@
#include "ensure_array_size.h"
#include "utils.h"
//#define VERBOSE
//#define YAC_VERBOSE_CLIPPING
static double const tol = 1.0e-12;
......@@ -126,7 +126,7 @@ void yac_compute_overlap_areas (unsigned N,
//partial_areas[n] = pole_area (overlap_buffer[n]);
}
#ifdef VERBOSE
#ifdef YAC_VERBOSE_CLIPPING
for (unsigned n = 0; n < N; n++)
printf("overlap area : %lf\n", partial_areas[n]);
#endif
......@@ -244,7 +244,7 @@ void yac_compute_concave_overlap_areas (unsigned N,
}
}
#ifdef VERBOSE
#ifdef YAC_VERBOSE_CLIPPING
for (unsigned n = 0; n < N; n++)
printf("overlap area %i: %lf \n", n, partial_areas[n]);
#endif
......@@ -1071,7 +1071,7 @@ void yac_correct_weights ( unsigned nSourceCells, double * weight ) {
unsigned iter;
double weight_diff;
#ifdef VERBOSE
#ifdef YAC_VERBOSE_CLIPPING
double weight_sum;
#endif
......@@ -1082,7 +1082,7 @@ void yac_correct_weights ( unsigned nSourceCells, double * weight ) {
for ( n = 0; n < nSourceCells; n++ )
weight_diff -= weight[n];
#ifdef VERBOSE
#ifdef YAC_VERBOSE_CLIPPING
for ( n = 0; n < nSourceCells; n++ )
weight_sum += weight[n];
......@@ -1098,7 +1098,7 @@ void yac_correct_weights ( unsigned nSourceCells, double * weight ) {
for ( n = 0; n < nSourceCells; n++ )
weight[n] += weight[n] * weight_diff;
}
#ifdef VERBOSE
#ifdef YAC_VERBOSE_CLIPPING
if ( fabs(weight_diff) > tol )
printf ("weight sum is %.15f \n", weight_sum);
#endif
......
......@@ -92,7 +92,7 @@ void yac_free_grid_cell(struct grid_cell * cell) {
yac_init_grid_cell(cell);
}
#ifdef DEBUG
#ifdef YAC_DEBUG_GRID_CELL
void print_grid_cell(FILE * stream, struct grid_cell cell, char * name) {
char * out = NULL;
......
......@@ -76,7 +76,7 @@ void yac_copy_grid_cell(struct grid_cell in_cell, struct grid_cell * out_cell);
*/
void yac_free_grid_cell(struct grid_cell * cell);
#ifdef DEBUG
#ifdef YAC_DEBUG_GRIC_CELL
/**
* prints out info about a grid_cell object and reinitialised, used for debugging
* interpolation_method_conserv and interpolation_method_patch
......
......@@ -1481,6 +1481,8 @@ int yac_gcxlatc(struct edge edge_a, struct edge edge_b,
* p and q will be identically, but only the p bits will be set
**/
#pragma intel optimization_level 0
int yac_gcxlatc_vec(double a[3], double b[3], double c[3], double d[3],
double p[3], double q[3]) {
......
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