Measure.C File Reference

#include "InfoStream.h"
#include "Measure.h"
#include "Node.h"
#include "Parameters.h"
#include "Molecule.h"

Go to the source code of this file.

Functions
static int	Tcl_centerOfNumber (ClientData, Tcl_Interp *interp, int argc, const char *argv[])
static int	Tcl_centerOfMass (ClientData, Tcl_Interp *interp, int argc, const char *argv[])
static int	Tcl_radiusOfGyration (ClientData, Tcl_Interp *interp, int argc, const char *argv[])
static int	Tcl_loadCoords (ClientData, Tcl_Interp *interp, int objc, Tcl_Obj *const objv[])

Function Documentation

Tcl_centerOfMass()

static int Tcl_centerOfMass ( ClientData  , Tcl_Interp *  interp, int  argc, const char *  argv[]  )

static

Definition at line 55 of file Measure.C.

References Molecule::atommass(), Node::coords, Node::molecule, Molecule::numAtoms, Node::Object(), Vector::x, Vector::y, and Vector::z.

Referenced by Measure::createCommands().

                                                                                          {

  Node *node = Node::Object();
  Molecule *molecule = node->molecule;
  // Parameters *parameters = node->parameters;
  Vector *coordinates = node->coords;
  int numAtoms = molecule->numAtoms;

  Vector center = 0;
  BigReal totalMass = 0;
  for( int i = 0; i < numAtoms; ++i ) {
    BigReal mass = molecule->atommass(i);
    totalMass += mass;
    center += mass * coordinates[i];
  }
  center /= totalMass;

  char s[1024];
  sprintf(s,"%g %g %g", center.x, center.y, center.z);
  Tcl_SetResult(interp,s,TCL_VOLATILE);

  return TCL_OK;
}

Tcl_centerOfNumber()

static int Tcl_centerOfNumber ( ClientData  , Tcl_Interp *  interp, int  argc, const char *  argv[]  )

static

Definition at line 32 of file Measure.C.

References Node::coords, Node::molecule, Molecule::numAtoms, Node::Object(), Vector::x, Vector::y, and Vector::z.

Referenced by Measure::createCommands().

                                                                                            {

  Node *node = Node::Object();
  Molecule *molecule = node->molecule;
  // Parameters *parameters = node->parameters;
  Vector *coordinates = node->coords;
  int numAtoms = molecule->numAtoms;

  int number = 0;
  Vector center = 0;
  for( int i = 0; i < numAtoms; ++i ) {
    number += 1;
    center += coordinates[i];
  }
  center /= number;

  char s[1024];
  sprintf(s,"%g %g %g", center.x, center.y, center.z);
  Tcl_SetResult(interp,s,TCL_VOLATILE);

  return TCL_OK;
}

Tcl_loadCoords()

static int Tcl_loadCoords ( ClientData  , Tcl_Interp *  interp, int  objc, Tcl_Obj *const  objv[]  )

static

Definition at line 111 of file Measure.C.

References coords, Node::coords, Node::molecule, Molecule::numAtoms, and Node::Object().

Referenced by Measure::createCommands().

                                                                                            {
  
  if (objc < 2 || objc > 3) {
    Tcl_SetResult(interp,(char*)"loadCoords: wrong # args",TCL_VOLATILE);
    return TCL_ERROR;
  }
  Tcl_Obj * const vname = objv[1];
  Node *node = Node::Object();
  Molecule *molecule = node->molecule;
  // Parameters *parameters = node->parameters;
  const Vector *coords = node->coords;
  int numAtoms = molecule->numAtoms;
 
  if (objc == 2) {
    // get all the coordinates
    for (int i=0; i<numAtoms; i++) {
      Tcl_Obj *newlist = Tcl_NewListObj(0, NULL);
      Tcl_Obj *arrkey = Tcl_NewIntObj(i+1);
      Tcl_IncrRefCount(arrkey);
      Tcl_ListObjAppendElement(interp, newlist, 
        Tcl_NewDoubleObj(coords[i].x));
      Tcl_ListObjAppendElement(interp, newlist, 
        Tcl_NewDoubleObj(coords[i].y));
      Tcl_ListObjAppendElement(interp, newlist, 
        Tcl_NewDoubleObj(coords[i].z));
      if (!Tcl_ObjSetVar2(interp, vname, arrkey, newlist, 0))
        return TCL_ERROR;
      Tcl_DecrRefCount(arrkey);
    }
  } else {
    // third argument must be a list of indices
    int nelems;
    Tcl_Obj **elems;
    if (Tcl_ListObjGetElements(interp, objv[2], &nelems, &elems) != TCL_OK) {
      return TCL_ERROR;
    }
    for (int i=0; i<nelems; i++) {
      int ind; 
      if (Tcl_GetIntFromObj(interp, elems[i], &ind) != TCL_OK)
        return TCL_ERROR;
      Tcl_Obj *newlist = Tcl_NewListObj(0, NULL);
      Tcl_Obj *arrkey = Tcl_NewIntObj(ind);
      Tcl_IncrRefCount(arrkey);
      --ind;
      Tcl_ListObjAppendElement(interp, newlist, 
        Tcl_NewDoubleObj(coords[ind].x));
      Tcl_ListObjAppendElement(interp, newlist, 
        Tcl_NewDoubleObj(coords[ind].y));
      Tcl_ListObjAppendElement(interp, newlist, 
        Tcl_NewDoubleObj(coords[ind].z));
      if (!Tcl_ObjSetVar2(interp, vname, arrkey, newlist, 0))
        return TCL_ERROR;
      Tcl_DecrRefCount(arrkey);
    }
  }
  return TCL_OK;
}

Tcl_radiusOfGyration()

static int Tcl_radiusOfGyration ( ClientData  , Tcl_Interp *  interp, int  argc, const char *  argv[]  )

static

Definition at line 79 of file Measure.C.

References Molecule::atommass(), Node::coords, Node::molecule, Molecule::numAtoms, and Node::Object().

Referenced by Measure::createCommands().

                                                                                              {

  Node *node = Node::Object();
  Molecule *molecule = node->molecule;
  // Parameters *parameters = node->parameters;
  Vector *coordinates = node->coords;
  int numAtoms = molecule->numAtoms;

  Vector center = 0;
  BigReal totalMass = 0;
  int i;
  for( i = 0; i < numAtoms; ++i ) {
    BigReal mass = molecule->atommass(i);
    totalMass += mass;
    center += mass * coordinates[i];
  }
  center /= totalMass;

  BigReal moment = 0;
  for( i = 0; i < numAtoms; ++i ) {
    BigReal mass = molecule->atommass(i);
    moment += mass * (coordinates[i] - center).length2();
  }
  BigReal radius = sqrt(moment/totalMass);

  char s[1024];
  sprintf(s,"%g", radius);
  Tcl_SetResult(interp,s,TCL_VOLATILE);

  return TCL_OK;
}