Punishments for Analytic Potentials


For analytic potentials it may sometimes be necessary to use additional punishments to keep the potential in a predefined parameter space. If the initial contraints, given in the parameter file can not achieve this, then there is a special function which is called before every force calculation. This function can be used to check for additional contraints on analytic potentials.

The function is called apot_punish and can be found in functions.c.

/****************************************************************
 *
 * punish analytic potential for bad habits
 *
 ****************************************************************/
 
real apot_punish(real *params, real *forces)

The parameters are

real *params
&nbsp;&nbsp;&nbsp;This is the array with the new parameters, access it with params[idx[i]], where $i<\text{ndim}$.

real *forces
&nbsp;&nbsp;&nbsp;This is the global forces array, store additional punishments to forces[punish_par_p+ i], where $i<\text{ndim}$

Here is an example, where a function is punished if two parameters are too close together:

/* punish m=n for universal embedding function */
    if (strcmp(apot_table.names[i], "universal") == 0) {
      x = params[j+ 2] - params[j+ 1];
      if (fabs(x) < 1e-6) {
	forces[punish_pot_p+ i] = apot_punish_value / (x * x);
	tmpsum += apot_punish_value / (x * x);
      }
    }

The return value is used to make potfit aware of these extra punishments.