/Users/craigcornelius/Projects/SPRING Mac Release 0.2/interpol.h

Go to the documentation of this file.

// $Id: interpol.h,v 1.9 2002/01/02 01:40:46 kevin Exp $
// $Copyright: (c)2001 National Biocomputation Center, Stanford University $

#ifndef _INTERPOL_H_
#define _INTERPOL_H_

#include <stdio.h>

#include "point3d.h"
#include "reallocablearray.h"
#include "space.h"
#include "socket.h"

class point {

 public :

  Point3D p;
  Point3D f;

  point(const Point3D np, const Point3D nf) : p(np), f(nf)
           {
           }
  point(const Point3D *np, const Point3D *nf)
           {
            p = *np;
            f = *nf;
           }
 point() : p(), f()
 {
 }

 point(const point *cp)
 {
  p = cp->p;
  f = cp->f;
 }
 
 ~point()
 {
 }

 long gethash()
 {
  long res;

  res = hashdouble(p.x);
  res ^= hashdouble(p.y);
  res ^= hashdouble(p.z);
  res ^= hashdouble(f.x);
  res ^= hashdouble(f.y);
  res ^= hashdouble(f.z);
  return res;
 }
 
 void print_datas()
 {
  fprintf(stderr, "Point class.\np: %f,%f,%f\nf: %f,%f,%f\n",
                         p.x,p.y,p.z,f.x,f.y,f.z);
 }
 
 void receive_data(char *buf)
 {
  int j=0;

  receive_float_from_buf(buf + j, &(p.x));
  j+=4;
  receive_float_from_buf(buf + j, &(p.y));
  j+=4;
  receive_float_from_buf(buf + j, &(p.z));
  j+=4;
  receive_float_from_buf(buf + j, &(f.x));
  j+=4;
  receive_float_from_buf(buf + j, &(f.y));
  j+=4;
  receive_float_from_buf(buf + j, &(f.z));
  j+=4;
 }
 
 void data_to_send(char *buf)
 {
  int j=0;

  copy_float_to_buf(buf + j, p.x);
  j+=4;
  copy_float_to_buf(buf + j, p.y);
  j+=4;
  copy_float_to_buf(buf + j, p.z);
  j+=4;
  copy_float_to_buf(buf + j, f.x);
  j+=4;
  copy_float_to_buf(buf + j, f.y);
  j+=4;
  copy_float_to_buf(buf + j, f.z);
  j+=4;
 }

};

//typedef std::vector<point *> point_l;
typedef ReallocableArray<point *> point_l;

class space_i : public space {
 public:
  point_l my_points;

  space_i() { }

  virtual ~space_i() { erase_all(); }

  virtual point *isthere(Point3D *newpoint)
  {
   for(int i=0; i<my_points.NumElements(); i++) {
//    if(((*i)->p - *newpoint).Squared() < 0.00001)
        point* cur = my_points[i];
    if(cur->p == *newpoint)
      return (cur);
   }
   return 0;
  }
  
  void add_new_point(point *p)
  {
   my_points.Append(p);
  }

  virtual void print_datas()
  {
   fprintf(stderr, "Space of interpolation.\n");
   for(int i=0; i<my_points.NumElements(); i++) {
          point* cur = my_points[i];
      cur->print_datas();
        }
  }
  
  virtual void send_datas(Socket *s)
  {
   unsigned short size;
   int j;
   char *buf;
   us_t l;

   size = 24 * my_points.NumElements();
   buf = (char *) malloc(sizeof(char) * (3 + size));
   if(buf == NULL)
    return;
   l.s = htons(size);
   buf[0] = TYPE_INTERPOL;
   for(int k=0; k<2; k++)
    buf[k+1] = l.c[k];
   j=0;
   for(int i=0; i<my_points.NumElements(); i++) {
          point* cur = my_points[i];
      cur->data_to_send(buf + 3 + (24*j));
    j++;
   }
   s->Send(buf, 3 + size);
   free(buf);
  }
  
  virtual void use_new_datas(char *buf, int size)
  {
   int i,nb;
   point *p;
   
   nb = size / 24;
   
   erase_all();

   for(i=0;i<nb; i++)
   {
    p = new point();
    p->receive_data(buf);
    add_new_point(p);
    buf+=24;
   }
  }
 
  virtual void erase_all()
  {
#if 0
   for(int i=0; i<my_points.NumElements(); i++) {
          point* cur = my_points[i];
      delete cur;
   }
#endif
   my_points.DeleteAll();
  }

  long gethash()
  {
   long res = 0;

   for(int i=0; i<my_points.NumElements(); i++) {
          point* cur = my_points[i];
      res ^= (cur->gethash());
   }

   return res;
  }
  
  virtual Point3D *get_force(Point3D *pos, Point3D* res)
  {
   Point3D *tmp1,*tmp2;
   double distance, sumdistance = 0;

   if(res == NULL)
    res = new Point3D();
   
   *res = Point3D(0,0,0);
   for(int i=0; i<my_points.NumElements(); i++) {
          point* cur = my_points[i];
    tmp1 = &(cur)->p;
    distance = sqrt((*pos - *tmp1).Squared());
    if(distance == 0)
    {
     *res = *tmp1;
     return res;
    }
    tmp2 = &(cur)->f;
    *res += (1.0/distance * (*tmp2));
    sumdistance += (1.0/distance);
   }
  
   *res /= sumdistance;
   
   return res;
  }

  virtual int istouching()
  {
   return 1;
  }

};

#endif

---