/* See {rdo_linear_algebra.h}. */

#define rdo_linear_algebra_C_COPYRIGHT "Copyright © 2008 Danillo Pereira and J. Stolfi, UNICAMP"
/* Last edited on 2024-12-21 11:52:04 by stolfi */ 

/* !!! Rename this module to {rdo_linalg}. !!! */

#include <stdio.h>
#include <assert.h>
#include <stdlib.h>
#include <math.h>

#include <dspmat.h>
#include <bool.h>
#include <affirm.h>

#include <rdo_geometry.h>
#include <rdo_linear_algebra.h>

#define debug_linalg (1)
  /* Define it as {1,2,...} to get increasingly detailed debugging printouts. */

void rdo_linalg_set_some_elements(double v[], int n, int sel[], double w[], int m)
  { int i, j;
    for (i = 0; i < n; i++) { v[i] = 0.0; }
    for (j = 0; j < m; j++) 
      { i = sel[j];
        if ((i >= 0) && (i < n)) { v[i] = w[j]; }
      }
  }