Rewrite constraint evaluation code.

Make it much more efficient and easier to add other seed functions.

1 files changed, 11 insertions, 0 deletions

diff --git a/nlsolve.h b/nlsolve.h
index f9b69fe..426cb24 100644
--- a/nlsolve.h
+++ b/nlsolve.h
@@ -72,6 +72,17 @@ typedef struct NLSolveContext {
      */
     unsigned int (*solve_order)[2];
 
+    /**
+     * Locations of the collocation points. The equation coefficients in NLEqCallback and
+     * NLEqJacobianCallback should be evaluated at those grid positions.
+     *
+     * colloc_grid[i][j] is an array of length solve_order[i][j] and contains
+     * the collocation points for the i-th equation in the j-th direction.
+     *
+     * Set by the solver after tdi_nlsolve_context_init().
+     */
+    double *(*colloc_grid)[2];
+
 #if HAVE_OPENCL
     cl_context       ocl_ctx;
     cl_command_queue ocl_queue;