constraints.c

Go to the documentation of this file.

/*****************************************************************************/
#include "tpcclibConfig.h"
/*****************************************************************************/
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <string.h>
/*****************************************************************************/
#include "tpcextensions.h"
#include "tpcrand.h"
/*****************************************************************************/
#include "tpcnlopt.h"
/*****************************************************************************/
 
/*****************************************************************************/
unsigned int nloptCheckParameters(
  unsigned int n,
  double *pLower,
  double *pUpper,
  double *p,
  double *pAccept,
  double *penalty
) {
  unsigned int pi, nAccept=0;
  double range;
 
  if(penalty!=NULL) *penalty=1.0;
  if(n<1 || p==NULL || pLower==NULL || pUpper==NULL)
    return(nAccept);
  for(pi=0; pi<n; pi++) {
    range=pUpper[pi]-pLower[pi]; if(range<1.0E-30) range=1.0;
    if(p[pi]<pLower[pi]) {
      if(pAccept!=NULL) pAccept[pi]=pLower[pi];
      if(penalty!=NULL) *penalty += (pLower[pi]-p[pi])/range;
    } else if(p[pi]>pUpper[pi]) {
      if(pAccept!=NULL) pAccept[pi]=pUpper[pi];
      if(penalty!=NULL) *penalty += (p[pi]-pUpper[pi])/range;
    } else {
      if(pAccept!=NULL) pAccept[pi]=p[pi];
      nAccept++;
    }
  }
  return nAccept;
}
/*****************************************************************************/
 
/*****************************************************************************/
unsigned int nloptForceLimits(
  unsigned int n,
  double *pLower,
  double *pUpper,
  double *p
) {
  unsigned int pi, nAccept=0;
 
  if(n<1 || p==NULL || pLower==NULL || pUpper==NULL) return(nAccept);
  for(pi=0; pi<n; pi++) {
    if(!(p[pi]>=pLower[pi])) {
      p[pi]=pLower[pi];
    } else if(!(p[pi]<=pUpper[pi])) {
      p[pi]=pUpper[pi];
    } else {
      nAccept++;
    }
  }
  return(nAccept);
}
/*****************************************************************************/
 
/*****************************************************************************/