rndpoint.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"
/*****************************************************************************/
 
/*****************************************************************************/
int nloptRandomPoint(
  double *p,
  double *low,
  double *up,
  unsigned int n,
  MERTWI *mt
) {
  if(p==NULL || low==NULL || up==NULL) return(1);
  if(n==0) return(0);
 
  if(mt==NULL) {
    for(unsigned int i=0; i<n; i++) {
      double dif=up[i]-low[i];
      if(dif<=0.0) p[i]=low[i]; else p[i]=low[i]+drand()*dif;
    }
  } else {
    for(unsigned int i=0; i<n; i++) {
      double dif=up[i]-low[i];
      if(dif<=0.0) p[i]=low[i]; else p[i]=low[i]+mertwiRandomDouble1(mt)*dif;
    }
  }
 
  return(0);
}
/*****************************************************************************/
 
/*****************************************************************************/
int nloptGaussianPoint(
  double *p,
  double *mean,
  double *sd,
  double *low,
  double *up,
  unsigned int n,
  MERTWI *mt
) {
  if(p==NULL || mean==NULL || sd==NULL) return(1);
  if(n==0) return(0);
 
  for(unsigned int i=0; i<n; i++) {
    if(mt==NULL) p[i]=mean[i]+sd[i]*drandGaussian();
    else p[i]=mean[i]+sd[i]*mertwiRandomGaussian(mt);
    if(low!=NULL && p[i]<low[i]) p[i]=low[i];
    if(up!=NULL  && p[i]>up[i])  p[i]=up[i];
  }
 
  return(0);
}
/*****************************************************************************/
 
/*****************************************************************************/