/* See {rdo_radiosity.h} */
#define rdo_radiosity_C_COPYRIGHT "Copyright © 2008 Danillo Pereira and J. Stolfi, UNICAMP"
/* Last edited on 2024-12-21 11:51:56 by stolfi */ 

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

#include <frgb.h>
#include <argparser.h>
#include <affirm.h>
#include <dspmat.h>

#include <rdo_geometry.h>
#include <rdo_solid.h>
#include <rdo_camera.h>
#include <rdo_radiosity.h>
#include <rdo_scene.h>

void rdo_radiosity_initialize_basis_radiance(double E[], double L[], int n)
  { int i;
    for(i = 0; i < n; i++) { L[i] = E[i]; }
  } 

void rdo_radiosity_basis_transfer_iteration(double E[], dspmat_t *R, double L[], int n)
  { assert(R->rows == n);
    assert(R->cols == n);
    double *c = NOTNULL(malloc(n * sizeof(double)));
    dspmat_map_col(R, L, n, c, n);     
    int i;
    for(i = 0; i < n; i++) { L[i] = E[i] + c[i]; }
    free(c);
  }