---

gnodederived.h

Go to the documentation of this file.

//****************************************************
//  April, 1993, University of Illinois
// Copyright (C) 1993, 1994 Tianlin Wang
/* Copyright (C) 1994-2003 Matvec Development Team. 

  This program is free software; you can redistribute it and/or
  modify it under the terms of the GNU Library General Public
  License as published by the Free Software Foundation; either
  version 2 of the License, or (at your option) any later version.
  
  This program is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  Library General Public License for more details.
    
  You should have received a copy of the GNU Library General Public
  License along with this library; if not, write to the Free
  Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
  MA 02111-1307, USA 
*/

#ifndef gnodederived_h
#define gnodederived_h
#include <iostream>
#include <ext/hash_map>
#include <map>
#include <set>
#include "exception.h"
#include "gnodestuff.h"
#include "safe_vectors.h"

namespace matvec {
  class AlleleStateNode:public GNode{
  public:
          unsigned alleleState,acceptedAlleleState,candidateAlleleState;
          AlleleStateNode(void){
                  alleleStateVector.name = "allelStateVector";
                  acceptedAlleleStateVector.name = "acceptedAlleleStateVector";
          }
          
          SafeSTLVector<unsigned> alleleStateVector;
          SafeSTLVector<unsigned> acceptedAlleleStateVector;  
          inline unsigned getState(void) {
                  return alleleState;
          };
          inline unsigned getOldState(void) {
                  return acceptedAlleleState;
          };
          inline unsigned getAcceptedAlleleState(void){
                  return alleleStateVector[acceptedAlleleState]-1;
          };
          inline unsigned getMyAlleleState(void){
                  return alleleStateVector[alleleState]-1;
          };
          inline unsigned getWeight(void) {
                  return alleleStateVector.size();
          };
          void reset(int i){
                  alleleState = i;
          };
          void resetAndSwitch(int i){
                  acceptedAlleleState = alleleState; 
                  alleleState = i;
          };
          bool incr();
          ~AlleleStateNode(void) {;}
  };
        /*! \class AlleleStateNode gnodederived.h inc/gnodederived.h
   * \brief This is the allele state node class.
   * 
   * AlleleStateNode class has methods relevant to an allele state node
   */ 
  class AlleleOriginNode:public GNode{
  public:
    unsigned alleleOrigin,oldAlleleOrigin;
    SafeSTLVector<unsigned> alleleOriginVector;
    inline unsigned getState(void) {
      return alleleOrigin;
    };
    inline unsigned getOldState(void) {
      return oldAlleleOrigin;
    };
    inline unsigned getMyAlleleOrigin(void){
      return alleleOriginVector[alleleOrigin];
    };
    inline unsigned getWeight(void) {
      return alleleOriginVector.size();
    };
    void reset(int i){
      alleleOrigin = i;
    };
    void resetAndSwitch(int i){
      oldAlleleOrigin = alleleOrigin; 
      alleleOrigin = i;
    };
    bool incr();
    ~AlleleOriginNode(void) {;}
  };
  /*! \class AlleleOriginNode gnodederived.h inc/gnodederived.h
   * \brief This is the allele origin node class.
   * 
   * AlleleOriginNode class has methods relevant to an allele origin node
   */ 
  class MaternalPaternalAlleles {
  public:
    unsigned maternal;
    unsigned paternal;
        bool operator< (const MaternalPaternalAlleles y) const;
  };
  /*! \class MaternalPaternalAlleles gnodederived.h inc/gnodederived.h
   * \brief This is the building block class for the GenotypeNode class.
   *
   * MaternalPaternalAlleles class defines the maternal and paternal 
   * alleles of a genotype.
   */ 
  class GenotypeNode:public GNode{
public:
          unsigned genotypeState,candidateGenotypeState,acceptedGenotypeState;
          GenotypeNode(void){
                  genotypeVector.name = "genotypeVector";
                  acceptedGenotypeVector.name = "acceptedGenotypeVector";
          }
          SafeSTLVector<MaternalPaternalAlleles> genotypeVector;
          SafeSTLVector<unsigned> genotypeCount;
          SafeSTLVector<MaternalPaternalAlleles> acceptedGenotypeVector;
          
          inline unsigned getState(void) {
                  return genotypeState;
          };
          inline unsigned getOldState(void) {
                  return acceptedGenotypeState;
          };
          inline unsigned getAcceptedMatState(void){
                  return acceptedGenotypeVector[acceptedGenotypeState].maternal;
          };
          inline unsigned getAcceptedPatState(void){
                  return acceptedGenotypeVector[acceptedGenotypeState].paternal;
          };
          inline unsigned getmState(void){
                  return genotypeVector[genotypeState].maternal;
          };
          inline unsigned getpState(void){
                  return genotypeVector[genotypeState].paternal;
          };
          inline unsigned getWeight(void) {
                  return genotypeVector.size();
          };
          void reset(int i){
                  genotypeState = i;
          };
          void resetAndSwitch(int i){
                  acceptedGenotypeState = genotypeState; 
                  acceptedGenotypeVector = genotypeVector;
                  genotypeState = i;
          };
          bool incr();
          ~GenotypeNode(void) {;}
  };
  /*! \class GenotypeNode gnodederived.h inc/gnodederived.h
   * \brief This is the genotype node class.
   *
   * GenotypeNode class has methods relevant to a genotype node
   */ 
}////// end of namespace matvec 
#endif

---