casadi 3.5.6.pre2__cp39-none-win_amd64.whl → 3.6.0__cp39-none-win_amd64.whl

casadi/include/casadi/core/runtime/casadi_jac.hpp

@@ -0,0 +1,186 @@
+//
+//    MIT No Attribution
+//
+//    Copyright (C) 2010-2023 Joel Andersson, Joris Gillis, Moritz Diehl, KU Leuven.
+//
+//    Permission is hereby granted, free of charge, to any person obtaining a copy of this
+//    software and associated documentation files (the "Software"), to deal in the Software
+//    without restriction, including without limitation the rights to use, copy, modify,
+//    merge, publish, distribute, sublicense, and/or sell copies of the Software, and to
+//    permit persons to whom the Software is furnished to do so.
+//
+//    THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,
+//    INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
+//    PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
+//    HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+//    OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
+//    SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+//
+
+
+// SYMBOL "jac_prob"
+template<typename T1>
+struct casadi_jac_prob {
+  // Number of outputs, i.e. rows of the Jacobian
+  casadi_int n_out;
+  // Number of inputs, i.e. columns of the Jacobian
+  casadi_int n_in;
+  // Number of colors
+  casadi_int n_color;
+  // Extended Jacobian sparsity
+  const casadi_int* sp_ext;
+  // Jacobian coloring
+  const casadi_int* coloring;
+  // Nominal values for inputs, if any
+  const T1* nom_in;
+  // Index mapping for outputs (i.e. Jacobian rows), if any
+  const size_t* map_out;
+  // Index mapping for inputs (i.e.  Jacobian columns), if any
+  const size_t* map_in;
+};
+// C-REPLACE "casadi_jac_prob<T1>" "struct casadi_jac_prob"
+
+// SYMBOL "jac_setup"
+template<typename T1>
+void casadi_jac_setup(casadi_jac_prob<T1>* p, const casadi_int* sp_ext,
+    const casadi_int* coloring) {
+  // Set pointers
+  p->sp_ext = sp_ext;
+  p->coloring = coloring;
+  // Dimensions are given by the sparsity patterns
+  p->n_out = sp_ext[0];
+  p->n_in = sp_ext[1];
+  p->n_color = coloring[1];
+  // The following defaults to null
+  p->nom_in = 0;
+  p->map_out = 0;
+  p->map_in = 0;
+}
+
+// SYMBOL "jac_work"
+template<typename T1>
+void casadi_jac_work(const casadi_jac_prob<T1>* p, casadi_int* sz_iw, casadi_int* sz_w) {
+  // Reset sz_w, sz_iw
+  *sz_w = *sz_iw = 0;
+  // Work vectors in data struct
+  *sz_iw += p->n_in;  // iseed
+  *sz_w += p->n_in;  // seed
+  *sz_iw += p->n_out;  // isens
+  *sz_w += p->n_out;  // sens
+  *sz_w += p->n_out;  // scal
+  *sz_iw += p->n_out;  // wrt
+  *sz_iw += p->n_out;  // nzind
+}
+
+// SYMBOL "jac_data"
+template<typename T1>
+struct casadi_jac_data {
+  // Number of seeds, sensitivities for the current color
+  casadi_int nseed, nsens;
+  // Inputs that are being seeded
+  casadi_int *iseed;
+  // Set of seeds for the seeded inputs
+  T1 *seed;
+  // Set of outputs for which sensitivities are calculated
+  casadi_int *isens;
+  // Set of values for the calculated sensitivities
+  T1 *sens;
+  // Scaling factors for calculated sensitivities
+  T1 *scal;
+  // Input corresponding to calculated sensitivities
+  casadi_int *wrt;
+  // Jacobian nonzero corresponding to calculated sensitivities
+  casadi_int *nzind;
+};
+// C-REPLACE "casadi_jac_data<T1>" "struct casadi_jac_data"
+
+// SYMBOL "jac_init"
+template<typename T1>
+void casadi_jac_init(const casadi_jac_prob<T1>* p, casadi_jac_data<T1>* d,
+    casadi_int** iw, T1** w) {
+  // Set work vectors
+  d->iseed = *iw; *iw += p->n_in;
+  d->seed = *w; *w += p->n_in;
+  d->isens = *iw; *iw += p->n_out;
+  d->sens = *w; *w += p->n_out;
+  d->scal = *w; *w += p->n_out;
+  d->wrt = *iw; *iw += p->n_out;
+  d->nzind = *iw; *iw += p->n_out;
+}
+
+// SYMBOL "jac_pre"
+template<typename T1>
+void casadi_jac_pre(const casadi_jac_prob<T1>* p, casadi_jac_data<T1>* d, casadi_int c) {
+  // Local variables
+  casadi_int i, kc, vin, vout, Jk;
+  double nom, inv_nom;
+  const casadi_int *color_colind, *color_row, *jac_colind, *jac_row;
+  // Extract sparsities
+  color_colind = p->coloring + 2;
+  color_row = color_colind + p->n_color + 1;
+  jac_colind = p->sp_ext + 2;
+  jac_row = jac_colind + p->n_in + 1;
+  // Loop over input indices for color
+  d->nseed = d->nsens = 0;
+  for (kc = color_colind[c]; kc < color_colind[c + 1]; ++kc) {
+    vin = color_row[kc];
+    // Nominal value, used as a seed for the column
+    nom = p->nom_in ? p->nom_in[vin] : 1;
+    inv_nom = 1. / nom;
+    // Collect seeds for column
+    d->seed[d->nseed] = nom;
+    d->iseed[d->nseed] = vin;
+    d->nseed++;
+    // Request corresponding outputs
+    for (Jk = jac_colind[vin]; Jk < jac_colind[vin + 1]; ++Jk) {
+      vout = jac_row[Jk];
+      d->scal[d->nsens] = inv_nom;
+      d->isens[d->nsens] = vout;
+      d->wrt[d->nsens] = vin;
+      d->nzind[d->nsens] = Jk;
+      d->nsens++;
+    }
+  }
+  // Map indices
+  if (p->map_in) {
+    for (i = 0; i < d->nseed; ++i) d->iseed[i] = p->map_in[d->iseed[i]];
+    for (i = 0; i < d->nsens; ++i) d->wrt[i] = p->map_in[d->wrt[i]];
+  }
+  if (p->map_out) {
+    for (i = 0; i < d->nsens; ++i) d->isens[i] = p->map_out[d->isens[i]];
+  }
+}
+
+// SYMBOL "jac_scale"
+template<typename T1>
+void casadi_jac_scale(const casadi_jac_prob<T1>* p, casadi_jac_data<T1>* d) {
+  // Local variables
+  casadi_int i;
+  // Scale derivatives
+  for (i = 0; i < d->nsens; ++i) d->sens[i] *= d->scal[i];
+}
+
+// SYMBOL "get_sub"
+template<typename T1>
+void casadi_get_sub(T1* sub, const casadi_int* sp_a, const T1* nz_a,
+    casadi_int rbegin, casadi_int rend, casadi_int cbegin, casadi_int cend) {
+  // Local variables
+  casadi_int nc, r, c, k;
+  const casadi_int *colind, *row;
+  // Quick return if null
+  if (sub == 0) return;
+  // Extract sparsity
+  nc = sp_a[1];
+  colind = sp_a + 2;
+  row = colind + nc + 1;
+  // Loop over columns
+  for (c = cbegin; c < cend; ++c) {
+    // Loop over nonzeros
+    for (k = colind[c]; k < colind[c + 1]; ++k) {
+      // Get row
+      r = row[k];
+      // Save, if in range
+      if (r >= rbegin && r < rend) *sub++ = nz_a[k];
+    }
+  }
+}

casadi/include/casadi/core/runtime/casadi_kkt.hpp

@@ -0,0 +1,67 @@
+//
+//    MIT No Attribution
+//
+//    Copyright (C) 2010-2023 Joel Andersson, Joris Gillis, Moritz Diehl, KU Leuven.
+//
+//    Permission is hereby granted, free of charge, to any person obtaining a copy of this
+//    software and associated documentation files (the "Software"), to deal in the Software
+//    without restriction, including without limitation the rights to use, copy, modify,
+//    merge, publish, distribute, sublicense, and/or sell copies of the Software, and to
+//    permit persons to whom the Software is furnished to do so.
+//
+//    THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,
+//    INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
+//    PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
+//    HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+//    OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
+//    SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+//
+
+
+// SYMBOL "kkt"
+template<typename T1>
+void casadi_kkt(const casadi_int* sp_kkt, T1* nz_kkt,
+    const casadi_int* sp_h, const T1* nz_h,
+    const casadi_int* sp_a, const T1* nz_a,
+    const T1* S, const T1* D, T1* w, casadi_int* iw) {
+  // Local variables
+  casadi_int i, k, j, nx, nz;
+  const casadi_int *h_colind, *h_row, *a_colind, *a_row,
+                   *kkt_colind, *kkt_row;
+  // Extract sparsities
+  a_row = (a_colind = sp_a + 2) + (nx = sp_a[1]) + 1;
+  h_row = (h_colind = sp_h + 2) + nx + 1;
+  kkt_row = (kkt_colind = sp_kkt + 2) + (nz = sp_kkt[1]) + 1;
+  // Running indices for each row of A
+  for (i = nx; i < nz; ++i) iw[i - nx] = kkt_colind[i];
+  // Reset w to zero
+  casadi_clear(w, nz);
+  // Loop over columns of [H + D_x; A]
+  for (i=0; i<nx; ++i) {
+    // Copy scaled column of H to w
+    for (k=h_colind[i]; k<h_colind[i+1]; ++k) {
+      j = h_row[k];
+      w[j] = nz_h[k] * S[i] * S[j];
+    }
+    // Copy scaled column of A to w
+    for (k=a_colind[i]; k<a_colind[i+1]; ++k) {
+      j = a_row[k] + nx;
+      w[j] = nz_a[k] * S[i] * S[j];
+    }
+    // Add D_x to diagonal
+    w[i] += D[i];
+    // Copy column to KKT
+    for (k=kkt_colind[i]; k<kkt_colind[i+1]; ++k) {
+      j = kkt_row[k];
+      nz_kkt[k] = w[j];
+      if (j >= nx) nz_kkt[iw[j - nx]++] = w[j];
+    }
+    // Zero out w
+    for (k=h_colind[i]; k<h_colind[i+1]; ++k) w[h_row[k]] = 0;
+    for (k=a_colind[i]; k<a_colind[i+1]; ++k) w[a_row[k] + nx] = 0;
+  }
+  // Copy -D_g to diagonal
+  for (i=nx; i<nz; ++i) {
+    nz_kkt[iw[i - nx]++] = -D[i];
+  }
+}

casadi/include/casadi/core/runtime/casadi_kron.hpp

@@ -1,4 +1,22 @@
-// NOLINT(legal/copyright)
+//
+//    MIT No Attribution
+//
+//    Copyright (C) 2010-2023 Joel Andersson, Joris Gillis, Moritz Diehl, KU Leuven.
+//
+//    Permission is hereby granted, free of charge, to any person obtaining a copy of this
+//    software and associated documentation files (the "Software"), to deal in the Software
+//    without restriction, including without limitation the rights to use, copy, modify,
+//    merge, publish, distribute, sublicense, and/or sell copies of the Software, and to
+//    permit persons to whom the Software is furnished to do so.
+//
+//    THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,
+//    INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
+//    PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
+//    HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+//    OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
+//    SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+//
+
 // SYMBOL "kron"
 template<typename T1>
 void casadi_kron(const T1* a, const casadi_int* sp_a, const T1* b, const casadi_int* sp_b, T1* r) {

casadi/include/casadi/core/runtime/casadi_ldl.hpp

@@ -1,4 +1,22 @@
-// NOLINT(legal/copyright)
+//
+//    MIT No Attribution
+//
+//    Copyright (C) 2010-2023 Joel Andersson, Joris Gillis, Moritz Diehl, KU Leuven.
+//
+//    Permission is hereby granted, free of charge, to any person obtaining a copy of this
+//    software and associated documentation files (the "Software"), to deal in the Software
+//    without restriction, including without limitation the rights to use, copy, modify,
+//    merge, publish, distribute, sublicense, and/or sell copies of the Software, and to
+//    permit persons to whom the Software is furnished to do so.
+//
+//    THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,
+//    INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
+//    PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
+//    HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+//    OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
+//    SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+//
+
 // SYMBOL "ldl"
 // Calculate the nonzeros of the transposed L factor (strictly lower entries only)
 // as well as D for an LDL^T factorization

casadi/include/casadi/core/runtime/casadi_logsumexp.hpp

@@ -0,0 +1,41 @@
+//
+//    MIT No Attribution
+//
+//    Copyright (C) 2010-2023 Joel Andersson, Joris Gillis, Moritz Diehl, KU Leuven.
+//
+//    Permission is hereby granted, free of charge, to any person obtaining a copy of this
+//    software and associated documentation files (the "Software"), to deal in the Software
+//    without restriction, including without limitation the rights to use, copy, modify,
+//    merge, publish, distribute, sublicense, and/or sell copies of the Software, and to
+//    permit persons to whom the Software is furnished to do so.
+//
+//    THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,
+//    INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
+//    PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
+//    HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+//    OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
+//    SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+//
+
+// SYMBOL "logsumexp"
+template<typename T1>
+T1 casadi_logsumexp(const T1* x, casadi_int n) {
+  casadi_int i, max_ind;
+  T1 max, r;
+  if (n==1) return x[0];
+  max_ind = 0;
+  max = x[0];
+  // Determine max, argmax
+  for (i=1; i<n; ++i) {
+    if (x[i]>x[0]) {
+      max = x[i];
+      max_ind = i;
+    }
+  }
+  r = 0;
+  for (i=0; i<n; ++i) {
+    if (i!=max_ind) r += exp(x[i]-max);
+  }
+  // C-REPLACE "log1p" "casadi_log1p"
+  return log1p(r)+max;
+}

casadi/include/casadi/core/runtime/casadi_low.hpp

@@ -1,4 +1,22 @@
-// NOLINT(legal/copyright)
+//
+//    MIT No Attribution
+//
+//    Copyright (C) 2010-2023 Joel Andersson, Joris Gillis, Moritz Diehl, KU Leuven.
+//
+//    Permission is hereby granted, free of charge, to any person obtaining a copy of this
+//    software and associated documentation files (the "Software"), to deal in the Software
+//    without restriction, including without limitation the rights to use, copy, modify,
+//    merge, publish, distribute, sublicense, and/or sell copies of the Software, and to
+//    permit persons to whom the Software is furnished to do so.
+//
+//    THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,
+//    INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
+//    PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
+//    HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+//    OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
+//    SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+//
+
 
 // SYMBOL "low"
 template<typename T1>

casadi/include/casadi/core/runtime/casadi_lsqr.hpp

@@ -1,7 +1,25 @@
 // C-REPLACE "fabs" "casadi_fabs"
 // C-REPLACE "sign" "casadi_sign"
 
-// NOLINT(legal/copyright)
+//
+//    MIT No Attribution
+//
+//    Copyright (C) 2010-2023 Joel Andersson, Joris Gillis, Moritz Diehl, KU Leuven.
+//
+//    Permission is hereby granted, free of charge, to any person obtaining a copy of this
+//    software and associated documentation files (the "Software"), to deal in the Software
+//    without restriction, including without limitation the rights to use, copy, modify,
+//    merge, publish, distribute, sublicense, and/or sell copies of the Software, and to
+//    permit persons to whom the Software is furnished to do so.
+//
+//    THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,
+//    INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
+//    PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
+//    HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+//    OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
+//    SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+//
+
 // SYMBOL "lsqr_sym_ortho"
 template<typename T1>
 void casadi_lsqr_sym_ortho(T1 a, T1 b, T1* cs, T1* sn, T1* rho) {
@@ -27,7 +45,25 @@ void casadi_lsqr_sym_ortho(T1 a, T1 b, T1* cs, T1* sn, T1* rho) {
     }
 }
 
-// NOLINT(legal/copyright)
+//
+//    MIT No Attribution
+//
+//    Copyright (C) 2010-2023 Joel Andersson, Joris Gillis, Moritz Diehl, KU Leuven.
+//
+//    Permission is hereby granted, free of charge, to any person obtaining a copy of this
+//    software and associated documentation files (the "Software"), to deal in the Software
+//    without restriction, including without limitation the rights to use, copy, modify,
+//    merge, publish, distribute, sublicense, and/or sell copies of the Software, and to
+//    permit persons to whom the Software is furnished to do so.
+//
+//    THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,
+//    INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
+//    PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
+//    HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+//    OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
+//    SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+//
+
 // SYMBOL "lsqr_single_solve"
 // Ref: scipy
 template<typename T1>
@@ -91,6 +127,11 @@ int casadi_lsqr_single_solve(const T1* A, T1* x, casadi_int tr, const casadi_int
     rnorm = beta;
     arnorm = alpha * beta;
 
+    if (arnorm==0.0) {
+      casadi_clear(x, m);
+      return 0;
+    }
+
     while (itn<iter_lim) {
       itn++;
       for (i=0;i<m;++i) u[i]*=-alpha;
@@ -175,7 +216,25 @@ int casadi_lsqr_single_solve(const T1* A, T1* x, casadi_int tr, const casadi_int
     return 0;
 }
 
-// NOLINT(legal/copyright)
+//
+//    MIT No Attribution
+//
+//    Copyright (C) 2010-2023 Joel Andersson, Joris Gillis, Moritz Diehl, KU Leuven.
+//
+//    Permission is hereby granted, free of charge, to any person obtaining a copy of this
+//    software and associated documentation files (the "Software"), to deal in the Software
+//    without restriction, including without limitation the rights to use, copy, modify,
+//    merge, publish, distribute, sublicense, and/or sell copies of the Software, and to
+//    permit persons to whom the Software is furnished to do so.
+//
+//    THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,
+//    INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
+//    PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
+//    HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+//    OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
+//    SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+//
+
 // SYMBOL "lsqr_solve"
 template<typename T1>
 int casadi_lsqr_solve(const T1* A, T1* x, casadi_int nrhs, casadi_int tr,

casadi/include/casadi/core/runtime/casadi_masked_norm_inf.hpp

@@ -0,0 +1,33 @@
+//
+//    MIT No Attribution
+//
+//    Copyright (C) 2010-2023 Joel Andersson, Joris Gillis, Moritz Diehl, KU Leuven.
+//
+//    Permission is hereby granted, free of charge, to any person obtaining a copy of this
+//    software and associated documentation files (the "Software"), to deal in the Software
+//    without restriction, including without limitation the rights to use, copy, modify,
+//    merge, publish, distribute, sublicense, and/or sell copies of the Software, and to
+//    permit persons to whom the Software is furnished to do so.
+//
+//    THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,
+//    INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
+//    PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
+//    HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+//    OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
+//    SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+//
+
+// SYMBOL "masked_norm_inf"
+template<typename T1>
+T1 casadi_masked_norm_inf(casadi_int n, const T1* x, const casadi_int* mask) {
+  casadi_int i;
+  T1 ret = 0;
+// C-REPLACE "fmax" "casadi_fmax"
+  for (i=0; i<n; ++i) {
+    if (mask[i]) {
+      ret = fmax(ret, fabs(x[i]));
+    }
+  }
+
+  return ret;
+}

casadi/include/casadi/core/runtime/casadi_max_viol.hpp

@@ -1,4 +1,22 @@
-// NOLINT(legal/copyright)
+//
+//    MIT No Attribution
+//
+//    Copyright (C) 2010-2023 Joel Andersson, Joris Gillis, Moritz Diehl, KU Leuven.
+//
+//    Permission is hereby granted, free of charge, to any person obtaining a copy of this
+//    software and associated documentation files (the "Software"), to deal in the Software
+//    without restriction, including without limitation the rights to use, copy, modify,
+//    merge, publish, distribute, sublicense, and/or sell copies of the Software, and to
+//    permit persons to whom the Software is furnished to do so.
+//
+//    THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,
+//    INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
+//    PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
+//    HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+//    OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
+//    SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+//
+
 // SYMBOL "max_viol"
 template<typename T1>
 T1 casadi_max_viol(casadi_int n, const T1* x, const T1* lb, const T1* ub) {
@@ -11,6 +29,7 @@ T1 casadi_max_viol(casadi_int n, const T1* x, const T1* lb, const T1* ub) {
     x_i = x ? *x++ : zero;
     lb_i = lb ? *lb++ : zero;
     ub_i = ub ? *ub++ : zero;
+    // C-REPLACE "fmax" "casadi_fmax"
     r = fmax(r, fmax(x_i-ub_i, zero));
     r = fmax(r, fmax(lb_i-x_i, zero));
   }

casadi/include/casadi/core/runtime/casadi_mmax.hpp

@@ -0,0 +1,28 @@
+//
+//    MIT No Attribution
+//
+//    Copyright (C) 2010-2023 Joel Andersson, Joris Gillis, Moritz Diehl, KU Leuven.
+//
+//    Permission is hereby granted, free of charge, to any person obtaining a copy of this
+//    software and associated documentation files (the "Software"), to deal in the Software
+//    without restriction, including without limitation the rights to use, copy, modify,
+//    merge, publish, distribute, sublicense, and/or sell copies of the Software, and to
+//    permit persons to whom the Software is furnished to do so.
+//
+//    THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,
+//    INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
+//    PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
+//    HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+//    OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
+//    SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+//
+
+
+// C-REPLACE "std::numeric_limits<double>::infinity()" "casadi_inf"
+// SYMBOL "mmax"
+template<typename T1>
+T1 casadi_mmax(const T1* x, casadi_int n, casadi_int is_dense) {
+  T1 r = is_dense ? -std::numeric_limits<double>::infinity() : 0;
+  if (!x) return r;
+  return casadi_vfmax(x, n, r);
+}

casadi/include/casadi/core/runtime/casadi_mmin.hpp

@@ -0,0 +1,29 @@
+//
+//    MIT No Attribution
+//
+//    Copyright (C) 2010-2023 Joel Andersson, Joris Gillis, Moritz Diehl, KU Leuven.
+//
+//    Permission is hereby granted, free of charge, to any person obtaining a copy of this
+//    software and associated documentation files (the "Software"), to deal in the Software
+//    without restriction, including without limitation the rights to use, copy, modify,
+//    merge, publish, distribute, sublicense, and/or sell copies of the Software, and to
+//    permit persons to whom the Software is furnished to do so.
+//
+//    THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,
+//    INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
+//    PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
+//    HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+//    OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
+//    SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+//
+
+
+// C-REPLACE "std::numeric_limits<double>::infinity()" "casadi_inf"
+// SYMBOL "mmin"
+template<typename T1>
+T1 casadi_mmin(const T1* x, casadi_int n, casadi_int is_dense) {
+  T1 r = is_dense ? std::numeric_limits<double>::infinity() : 0;
+  if (!x) return r;
+  return casadi_vfmin(x, n, r);
+}
+

casadi/include/casadi/core/runtime/casadi_mtimes.hpp

@@ -1,4 +1,22 @@
-// NOLINT(legal/copyright)
+//
+//    MIT No Attribution
+//
+//    Copyright (C) 2010-2023 Joel Andersson, Joris Gillis, Moritz Diehl, KU Leuven.
+//
+//    Permission is hereby granted, free of charge, to any person obtaining a copy of this
+//    software and associated documentation files (the "Software"), to deal in the Software
+//    without restriction, including without limitation the rights to use, copy, modify,
+//    merge, publish, distribute, sublicense, and/or sell copies of the Software, and to
+//    permit persons to whom the Software is furnished to do so.
+//
+//    THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,
+//    INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
+//    PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
+//    HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+//    OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
+//    SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+//
+
 // SYMBOL "mtimes"
 template<typename T1>
 void casadi_mtimes(const T1* x, const casadi_int* sp_x, const T1* y, const casadi_int* sp_y, T1* z, const casadi_int* sp_z, T1* w, casadi_int tr) { // NOLINT(whitespace/line_length)

casadi/include/casadi/core/runtime/casadi_mv.hpp

@@ -1,4 +1,22 @@
-// NOLINT(legal/copyright)
+//
+//    MIT No Attribution
+//
+//    Copyright (C) 2010-2023 Joel Andersson, Joris Gillis, Moritz Diehl, KU Leuven.
+//
+//    Permission is hereby granted, free of charge, to any person obtaining a copy of this
+//    software and associated documentation files (the "Software"), to deal in the Software
+//    without restriction, including without limitation the rights to use, copy, modify,
+//    merge, publish, distribute, sublicense, and/or sell copies of the Software, and to
+//    permit persons to whom the Software is furnished to do so.
+//
+//    THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,
+//    INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
+//    PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
+//    HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+//    OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
+//    SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+//
+
 // SYMBOL "mv"
 template<typename T1>
 void casadi_mv(const T1* x, const casadi_int* sp_x, const T1* y, T1* z, casadi_int tr) {

casadi/include/casadi/core/runtime/casadi_mv_dense.hpp

@@ -1,4 +1,22 @@
-// NOLINT(legal/copyright)
+//
+//    MIT No Attribution
+//
+//    Copyright (C) 2010-2023 Joel Andersson, Joris Gillis, Moritz Diehl, KU Leuven.
+//
+//    Permission is hereby granted, free of charge, to any person obtaining a copy of this
+//    software and associated documentation files (the "Software"), to deal in the Software
+//    without restriction, including without limitation the rights to use, copy, modify,
+//    merge, publish, distribute, sublicense, and/or sell copies of the Software, and to
+//    permit persons to whom the Software is furnished to do so.
+//
+//    THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,
+//    INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
+//    PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
+//    HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+//    OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
+//    SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+//
+
 // SYMBOL "mv_dense"
 template<typename T1>
 void casadi_mv_dense(const T1* x, casadi_int nrow_x, casadi_int ncol_x,