---

Atom.h

Go to the documentation of this file.

/***************************************************************************
 *cr                                                                       
 *cr            (C) Copyright 1995-2008 The Board of Trustees of the           
 *cr                        University of Illinois                       
 *cr                         All Rights Reserved                        
 *cr                                                                   
 ***************************************************************************/

/***************************************************************************
 * RCS INFORMATION:
 *
 *        $RCSfile: Atom.h,v $
 *        $Author: johns $        $Locker:  $                $State: Exp $
 *        $Revision: 1.58 $        $Date: 2008/03/27 19:36:34 $
 *
 ***************************************************************************
 * DESCRIPTION:
 *
 * Structure which holds data for one atom.
 *
 ***************************************************************************/
#ifndef ATOM_H
#define ATOM_H

#include <string.h>
#include <stdlib.h>
#include "utilities.h"

// maximum number of bonds allowed to other atoms
#define MAXATOMBONDS      12

// Atom type flags
#define ATOMNORMAL      0
#define ATOMPROTEINBACK 1
#define ATOMNUCLEICBACK 2
#define ATOMHYDROGEN    3

// Residue type flags
#define RESNOTHING      0
#define RESPROTEIN      1
#define RESNUCLEIC      2
#define RESWATERS       3

struct MolAtom {
  // XXX contents of the Atom structure are ordered specifically so 
  // that the compiler will pack it efficiently.
 
  // items that make this particular atom unique and are absolutely 
  // needed to link it up to the rest of the structure, or are speed-critical
  short nameindex;              
  short typeindex;              
  int uniq_resid;               
  int bondTo[MAXATOMBONDS];     
  signed char bonds;            
  signed char atomicnumber;     
  signed char altlocindex;      
  char insertionstr[2];         

  // items which could potentially be moved into other data structures 
  // to save memory, but are presently kept here for extra simplicity or speed
  signed char chainindex;       
  short segnameindex;           
  int resid;                    
  short resnameindex;           

  // ATOMNORMAL, ATOMPROTEINBACK, ATOMNUCLEICBACK, ATOMHYDROGEN
  signed char atomType;         

  // RESNOTHING, RESPROTEIN, RESNUCLEIC, RESWATERS
  signed char residueType;      



  void init(int n, int theresid, char *insertion) { 
    int i;

    uniq_resid = 0; // don't know yet, found in BaseMolecule
    bonds = 0;
    resid = theresid;
    strncpy(insertionstr, insertion, 2);        insertionstr[1] = '\0';
    nameindex = typeindex = resnameindex = segnameindex = altlocindex = (-1);
    atomicnumber = (-1);

    for (i=0; i<MAXATOMBONDS; i++) {
      bondTo[i] = -1;
    }
    atomType = ATOMNORMAL;
    residueType = RESNOTHING;
  }
  
  int add_bond(int a, int type) {
    if(bonds >= MAXATOMBONDS) // fail
      return -1;

    bondTo[(int) bonds] = a;
    if (type == ATOMPROTEINBACK || type == ATOMNUCLEICBACK)
      atomType = type;
    bonds++;
    return 0;
  }

  int bonded(int a) {
    int i;
    for (i=0; i < bonds; i++) {
      if (bondTo[i] == a) {
        return TRUE;
      }
    }

    return FALSE;
  }

};

#endif

---