Пример #1
0
/* metropolis algorithm (with periodic checking) */
INLINE int cago_metro1(cago_t *go, real amp, real bet)
{
  int i, j;
  rv3_t xi;
  real du;
  static int cc;

  i = (int) (go->n * rnd0());
  for (j = 0; j < 3; j++) {
    xi[j] = go->x[i][j] + (rnd0() * 2.f - 1) * amp;
  }
  du = cago_depot(go, go->x, i, xi);
  if (du < 0 || rnd0() < exp(-bet * du)) {
    rv3_copy(go->x[i], xi);
    go->epot += du;
    if (++cc % 1000 == 0) { /* refresh energy */
      real ep = cago_force(go, go->x, NULL);
      if (fabs(ep - go->epot) > 0.01) {
        fprintf(stderr, "energy corruptions %g vs. %g\n", go->epot, ep);
      }
      go->epot = ep;
    }
    return 1;
  } else {
    return 0;
  }
}
Пример #2
0
/* mutate velocities by random rotation */
INLINE int md_rotv3d(rv3_t *v, int n)
{
  int i;
  real mat[3][3], v1[3], axis[3], ang;

  rv3_normalize(rv3_rnd(axis, -.5f, 1.f));
  ang = 2 * M_PI * rnd0(); /* amplitude has to be 2 * pi */
  rm3_mkrot(mat, axis, ang);
  for (i = 0; i < n; i++) {
    rm3_mulvec(v1, mat, v[i]);
    rv3_copy(v[i], v1);
  }
  return 0;
}
Пример #3
0
void randomizescrolls() {

    uchar numscrolls=sizeof(sname)/sizeof(sname[0]);
    char * tmp;
    uchar changeto=0;

    for(int i=0; i<numscrolls; i++) {
        changeto=rnd0(numscrolls-1);  // because numscrolls is 1 based.
        // move destination to tmp
        tmp=sname[changeto];
        // move current index to destination
        sname[changeto]=sname[i];
        // move tmp to current index
        sname[i]=tmp;
    }
}
Пример #4
0
objectclass *randomscroll() {

    struct sinfo {
        objectclass *(*fp)();       // Function that creates the object
        uchar prob;                 // probability of object being created.
    };

    struct sinfo sinfo[]= {
        // Put all the class create() function is here that are
        // possible candidates for this random object generator.
        // The sum of all the probabilities in this list must be 1000.
        // Think of it as XX.X% chance of being generated.
        {scroll_of_enchant_armor::create,63},
        {scroll_of_destroy_armor::create,45},
        {scroll_of_confuse_monster::create,53},
        {scroll_of_scare_monster::create,35},
        {scroll_of_remove_curse::create,65},
        {scroll_of_enchant_weapon::create,85},
        {scroll_of_create_monster::create,45},
        {scroll_of_taming::create,15},
        {scroll_of_genocide::create,15},
        {scroll_of_light::create,95},
        {scroll_of_teleportation::create,55},
        {scroll_of_gold_detection::create,33},
        {scroll_of_food_detection::create,25},
        {scroll_of_identify::create,185},
        {scroll_of_magic_mapping::create,45},
        {scroll_of_amnesia::create,35},
        {scroll_of_fire::create,48},
        {scroll_of_punishment::create,15},
        {scroll_of_charging::create,15},
        {scroll_of_blank_paper::create,28}
        // make sure the last one does not have a trailing comma
    };

    sshort prob;
    uchar count=0;

    prob=rnd0(1000);
    while((prob -= sinfo[count].prob) > 0) count++;
    return (sinfo[count].fp)();
}
Пример #5
0
int runSimulator(char *ip, int port) {

	int bear;
	int knots;
	int mps=0;
	int myTime =0;
	int lapse=defaultLapse;
	long x, y, d;
	float lat, lon;
	frm_cmd_gps_t gps;
	char buff[TRANS_MAX_BUFF_SIZE];
	size_t transLen;
	int seq =0;
	int sckt;
	struct sockaddr_in server;

	// INIT SOCKET
	printf("Incializando...\n");
	memset(&server,0,sizeof(server));

	sckt= socket(AF_INET,SOCK_DGRAM,0);

	if (sckt<0) {
			die("Error en la creacion del socket");
	}

	server.sin_family = AF_INET;

	server.sin_addr.s_addr = inet_addr(ip);
	server.sin_port = htons(port);



	y = lat2y(42.22);
	x = lon2x(-8.65);
	printf("Empezando simulacion...(lapso=%ds)\n\n",lapse);

	while (true) {

		if (myTime<=0) {
			myTime = rnd1(30);
			bear = rnd0(359);
			mps  = rnd1(MAX_MPS);
		}

		d = mps * lapse;

		// Ojo origen de rumbo PI
		y += round( d * cos(bear * DEG_RAD));
		x += round( d * sin(bear * DEG_RAD));

		lat = y2lat(y);
		lon = x2lon(x);
		knots = mps * MPS_TO_KNOTS;

//		gps.cmd = 0x11;
//		gps.len =0;
//		gps.seq_l= seq++;
//		gps.seq_s=0;
//		gps.bearing = GPS_ENCODE_BEARING(bear);
//		gps.knots = knots;
//		gps.lat_sign = lat<0?1:0;
//		gps.lat_deg = gps.lat_sign?-lat:lat;
//		gps.lat_min =htons(GPS_ENCODE_LOCMIN(lat));
//		gps.lon_sign = lon<0?1:0;
//		gps.lon_deg = gps.lon_sign?-lon:lon;
//		gps.lon_min = htons(GPS_ENCODE_LOCMIN(lon));
//
//		gps.fix = 1;
//		gps.hdop = rnd1(20);
//		gps.time.epoch = htonl(time(NULL));

		int offset=0;
		transLen = sizeof(buff);

		frame_encode_gps(seq++,bear,knots,lat,lon,1,rnd1(20),0,buff,&transLen);

		offset +=transLen;

		if (addCnxFrame) {

			transLen = sizeof(buff) - offset;
			frame_encode_cnx("1234567890",0,1,0,buff+offset,&transLen);
			offset +=transLen;

		}

		transLen = sizeof(buff);

		if ( frame_encode_transport(1,NULL,offset,buff,&transLen) ) {

//			if (frame_test_transport(buff,transLen)) {
//				printf("Trama de transporte verificada\n");
//			}

			printf("Enviando: %d lat=%f lon=%f bear=%d knots=%d time=%d\n",seq,
					lat,lon,bear,knots,ntohl(gps.time.epoch));

			sendto(sckt,buff,transLen,0,(struct sockaddr *)&server,sizeof(server));

		}

		myTime -=lapse;
		sleep(lapse);
	}

	close(sckt);
	return 0;

}
Пример #6
0
Файл: xymc.c Проект: 3ki5tj/alge
/* run md simulation */
static int run(void)
{
  int t, i, id, etot = 0, eacc = 0;
  double r[DIM], r1[DIM], ep, ep1, dep;
  double x0 = 0, y0 = 0, x1 = 0, y1 = 0, dx, dy, alf = 0;
  ad2_t *al;
  hist2_t *hs;

  al = ad2_open(-M_PI, M_PI, M_PI/9.999, -M_PI, M_PI, M_PI/9.999,
     0.001, 1.0, 0.0, 0.0);
  hs = hs2_opensqr1(-M_PI, M_PI, M_PI/9.999);

  for (i = 0; i < DIM; i++) {
    r[i] = (rnd0() * 2 - 1) * M_PI;
  }
  ep = force(r);
  id = ad2_getid(al, r[0], r[1]);

  for (t = 1; t <= nsteps; t++) {
    for (i = 0; i < DIM; i++) {
      r1[i] = WRAP( r[i] + (rnd0() - .5) * 1.0 );
    }
    ep1 = force(r1);
    dep = ep1 - ep;
    etot++;
    if (dep <= 0 || rnd0() < exp(-dep)) {
      eacc++;
      for (i = 0; i < DIM; i++) {
        r[i] = r1[i];
      }
      ep = ep1;
    }

    hs2_add1ez(hs, r[0], r[1], 0);

    /* register */
    if (t % seglen == 0) {
      static const double fcmax = 1000000000.0;

      x1 = r[0];
      y1 = r[1];
      dx = DIFF(x1, x0);
      dy = DIFF(y1, y0);
      ad2_add(al, id, dx, dy);
/*
      ad2_getcorr(al, id, &delx, &dely);
      dx += delx;
      dy += dely;
*/
      alf = 0.1;
      al->fx[id] = dblconfine(al->fx[id] + dx * alf, -fcmax, fcmax);
      al->fy[id] = dblconfine(al->fy[id] + dy * alf, -fcmax, fcmax);

      /* set the new start point */
      x0 = x1;
      y0 = y1;
      id = ad2_getid(al, x0, y0);

      if (t % nevery == 0) {
        printf("t %d, x %g, y %g\n", t, r[0], r[1]);
      }
      if (t % nreport == 0) {
        ad2_save(al, fnout);
      }
    }
  }
  printf("acc %g\n", 1.0*eacc/etot);
  ad2_close(al);
  hs2_save(hs, "xy.his", 0);
  return 0;
}
Пример #7
0
double rnd1(double s)
{
    return s*rnd0();
}
Пример #8
0
int runSimulator(char *ip, int port) {

	int bear;
	int knots;
	int mps=0;
	int myTime =0;
	int lapse=defaultLapse;
	long x, y, d;
	float lat, lon;
	frm_cmd_gps_t gps;
	char buff[TRANS_MAX_BUFF_SIZE];
	size_t transLen;
	int seq =0;
	int sckt;
	struct sockaddr_in server;
	time_t timestamp;


	PROJ *p;
	PROJ_XY xy;
	PROJ_LP lp;

	// INIT SOCKET
	printf("Incializando...\n");
	memset(&server,0,sizeof(server));

	sckt= socket(AF_INET,SOCK_DGRAM,0);

	if (sckt<0) {
			die("Error en la creacion del socket");
	}

	server.sin_family = AF_INET;

	server.sin_addr.s_addr = inet_addr(ip);
	server.sin_port = htons(port);


	//Projection initialization

	//WGS84
	p = proj_initstr("proj=merc ellps=WGS84");

	lp.phi = ORIG_LAT * D_TO_R;
	lp.lam = ORIG_LON * D_TO_R;

	xy = proj_fwd(lp,p);

	printf("Empezando simulacion...(lapso=%ds)\n\n",lapse);

	while (true) {

		if (myTime<=0) {
			myTime = rnd1(30);
			bear = rnd0(359);
			mps  = rnd1(MAX_MPS);
		}

		d = mps * lapse;

		// Ojo origen de rumbo PI
		xy.y += round( d * cos(bear * D_TO_R));
		xy.x += round( d * sin(bear * D_TO_R));

		lp = proj_inv(xy,p);

		knots = mps * MPS_TO_KNOTS;

		int offset=0;
		transLen = sizeof(buff);


		lat = lp.phi* R_TO_D;
		lon = lp.lam * R_TO_D;

		timestamp = time(NULL);

		frame_encode_gps(seq++,bear,knots,lat,lon,1,rnd1(20),timestamp,buff,&transLen);

		offset +=transLen;

		if (addCnxFrame) {

			transLen = sizeof(buff) - offset;
			frame_encode_cnx("1234567890",0,1,0,buff+offset,&transLen);
			offset +=transLen;

		}

		transLen = sizeof(buff);

		if ( frame_encode_transport(1,NULL,offset,buff,&transLen) ) {

//			if (frame_test_transport(buff,transLen)) {
//				printf("Trama de transporte verificada\n");
//			}

			printf("Enviando: %d lat=%f lon=%f bear=%d knots=%d time=%d\n",seq,
					lat,lon,bear,knots,timestamp);

			sendto(sckt,buff,transLen,0,(struct sockaddr *)&server,sizeof(server));

		}

		myTime -=lapse;
		sleep(lapse);
	}

	proj_free(p);
	close(sckt);
	return 0;

}