         return false;
     }
 
+#if 1
     // based on Bemap
     float Dx  = 0.0f;
     float Dy  = 0.0f;

    // ec.value   = 0;
    // ec.edge    = 0;
    ec.angle_from_bemap = 0;

    return true;
}

    }

    /* reject condition: tr(H)^2/det(H) < (r+1)^2/r */
    if( edgeval >= (edge_limit+1.0f)*(edge_limit+1.0f)/edge_limit ) {
        return false;
    }

