src/modelidentifier.cpp

   1 /***************************************************************************
   2  *   Copyright (C) 2006 by Michael Andreen                                 *
   3  *   andreen@student.chalmers.se                                           *
   4  *                                                                         *
   5  *   This program is free software; you can redistribute it and/or modify  *
   6  *   it under the terms of the GNU General Public License as published by  *
   7  *   the Free Software Foundation; either version 2 of the License, or     *
   8  *   (at your option) any later version.                                   *
   9  *                                                                         *
  10  *   This program is distributed in the hope that it will be useful,       *
  11  *   but WITHOUT ANY WARRANTY; without even the implied warranty of        *
  12  *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
  13  *   GNU General Public License for more details.                          *
  14  *                                                                         *
  15  *   You should have received a copy of the GNU General Public License     *
  16  *   along with this program; if not, write to the                         *
  17  *   Free Software Foundation, Inc.,                                       *
  18  *   51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA          *
  19  ***************************************************************************/
  20
  21 #include "modelidentifier.h"
  22 #include "genealgorithms.h"
  23
  24 #include <doublefann.h>
  25
  26 using namespace std;
  27
  28 ModelIdentifier::ModelIdentifier(std::string ann){
  29         _ann = fann_create_from_file(ann.c_str());
  30         if(!_ann){
  31                 throw invalid_argument("Could not create network");
  32         }
  33 }
  34
  35 ModelIdentifier::~ModelIdentifier(){
  36         fann_destroy(_ann);
  37 }
  38
  39 std::map<ModelIdentifier::Model,double> ModelIdentifier::identify(const GeneOrder& go){
  40         int pos = 0;
  41         int neg = 0;
  42         for (GeneOrder::iterator g = go.begin(); g != go.end(); ++g){
  43                 if (*g >= 0)
  44                         ++pos;
  45                 else
  46                         ++neg;
  47         }
  48         map<Model,double> scores;
  49         scores[X] = 1;
  50         scores[Whirl] = -1;
  51         return scores;
  52 }