sensitivity_filter.c

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#ifdef _OPENMP
#include <omp.h>
#endif
#include "../headers/function_header.h"

extern struct _PROBLEM *problem;
extern struct _FLAGS *flags;

void density_filter(FEM *actfem)
{

  int i, j, ID;
  double H, r, dis;

  PROBLEM *prob;
  ELEMENT *actele1, *actele2;
  NODE *actnode;

  prob = problem;
  /*--------------------------------------------------------- cal df / ds */

  for (i = 0; i < actfem->nelx; i++)
  {
    actele1 = &(actfem->element[i]);
    r = actele1->filter->r;
    /*-----------initialization------------------*/
      actele1->filobj = 0.0;
    /*summation of density filter*/
    for (j = 0; j < actele1->filter->numele; j++)
    {
      ID = actele1->filter->ID[j];
      dis = actele1->filter->dis[j];
      actele2 = &(actfem->element[ID]);
      
        H = r - dis;
        actele1->filobj += actele2->obj * H * actele2->filsensi;
    }
  }
}

obj...raw sensitivity of objective (or constraint) function
filobj...filtered sensitivity