MAP

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%   This program is free software; you can redistribute it and/or modify

%   it under the terms of the GNU 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 General Public License for more details.

%

%   You can obtain a copy of the GNU General Public License from

%   http://www.gnu.org/copyleft/gpl.html or by writing to

%   Free Software Foundation, Inc.,675 Mass Ave, Cambridge, MA 02139, USA.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

classdef cHMM

    %HMMCLASS Handles all of the HTK related gubbins

    %   Please see the documentation located in a separate file for further

    %   information.

    %%  *********************************************************

    %  properties                      _   _

    %                                 | | (_)

    %  _ __  _ __ ___  _ __   ___ _ __| |_ _  ___  ___

    % | '_ \| '__/ _ \| '_ \ / _ \ '__| __| |/ _ \/ __|

    % | |_) | | | (_) | |_) |  __/ |  | |_| |  __/\__ \

    % | .__/|_|  \___/| .__/ \___|_|   \__|_|\___||___/

    % | |             | |

    % |_|             |_|

    %************************************************************

    %% **********************************************************

    % Public properties - can be set by user

    %************************************************************

    properties(Access = public)

        hmmFolder

        paramType               = 'USER_D_A'; %DELTAS and ACCELERATIONS

        numCoeff                = 27; %9*3 THIS IS FOR PROBABILITY MODEL (not high spont+low spont which would be 18*3=54)

        HERestDataPath      = fullfile(pwd, 'def', 'HERest_digit');

        binPath             = fullfile(pwd, 'def', 'bin');

        configFile          = fullfile(pwd, 'def', 'config_STANDARD');

        trainWordListFile   = fullfile(pwd, 'def', 'Grammar_digit', 'words3');

        testWordListFile    = fullfile(pwd, 'def', 'Grammar_digit', 'wordsNoSevenZero');

        wordNetFile         = fullfile(pwd, 'def', 'Grammar_digit', 'wdnetNoSP.slf');

        dictFile            = fullfile(pwd, 'def', 'Grammar_digit', 'noSevenZeroDict');

    end

    %% **********************************************************

    % Dependent - never set by user. Only calculated when needed

    %************************************************************

    properties(Dependent = true)

        protoFile              % = fullfile(pwd, 'def', 'proto_RobANonly_9'); %probability only

    end

    %%  *********************************************************

    % methods        _   _               _

    %               | | | |             | |

    % _ __ ___   ___| |_| |__   ___   __| |___

    %| '_ ` _ \ / _ \ __| '_ \ / _ \ / _` / __|

    %| | | | | |  __/ |_| | | | (_) | (_| \__ \

    %|_| |_| |_|\___|\__|_| |_|\___/ \__,_|___/

    %************************************************************

    methods

        %% **********************************************************

        % Constructor

        %************************************************************

        function obj = cHMM(hmmFolder)

            if nargin > 0

                obj.hmmFolder = hmmFolder;

            end

        end % ------ OF CONSTRUCTOR

        %% **********************************************************

        % genProto - generate task specific prototype

        %************************************************************

        function genProto(obj)

            % models_1mixsil.exe - takes input (hmmdef) and copies it making a

            % basis of one, two three etc. etc.

            obj.models_1mixsilMat(fullfile(obj.hmmFolder,'hmm0','hmmdef'), fullfile(obj.hmmFolder,'hmm0','models'));

        end % ------ OF GENPROTO

        %% **********************************************************

        % istrained

        %************************************************************

        function boolans = istrained(obj)

            boolans = numel(dir(fullfile(obj.hmmFolder,'hmm36','models')));

        end

        %% **********************************************************

        % Train

        %************************************************************

        function train(obj, trainFeatureFolder)

            % Most of the following code block can be replaced by doing a

            % find and replace across the  code below. I didn't

            % want to mess with the working code too much, so I just

            % copied the object properties needed into the variable names

            % expected.

            ED_CMDFILE1 = fullfile(obj.HERestDataPath, 'sil1.hed');

            ED_CMDFILE2 = fullfile(obj.HERestDataPath, 'mix2_16.hed');

            ED_CMDFILE3 = fullfile(obj.HERestDataPath, 'mix3_16.hed');

            ED_CMDFILE4 = fullfile(obj.HERestDataPath, 'mix5_16.hed');

            ED_CMDFILE5 = fullfile(obj.HERestDataPath, 'mix7_16.hed');

            NUM_COEF        = obj.numCoeff;

            PAR_TYPE        = obj.paramType;

            LIST_FILE       = fullfile(obj.hmmFolder, 'tmp.list');

            word_list       = obj.trainWordListFile;%fullfile(obj.grammarPath, 'words');

            word_listSP     = word_list; % for use in hmm4 onwards - UGLY HACK NOW SP ABANDONED

            proto           = obj.protoFile;%Does not exist on disk just yet probably - see a few lines down

            config          = obj.configFile;

            train_list      = fullfile(trainFeatureFolder, 'list.scp');

            labels          = fullfile(trainFeatureFolder, 'labels.mlf');

            labelssp        = labels; % for use in hmm4 onwards - UGLY HACK NOW SP ABANDONED

            hmm_dir         = obj.hmmFolder;

            FEAT_ROOT       = trainFeatureFolder;

            % Now for the actual HMM training code

            mkdir(hmm_dir)

            for I = 0:36

                h = fullfile(hmm_dir,['hmm' num2str(I)]);

                mkdir(h);

            end

            obj.makeProtoHmm(proto, obj.paramType, obj.numCoeff, 18);

            fid = fopen(train_list,'r');

            disp(train_list)

            S = textscan(fid,'%s','Delimiter','\n');

            fclose(fid);

            fid = fopen(LIST_FILE,'w');

            for I = 1:size(S{1},1)

                str  = fullfile(FEAT_ROOT,S{1}{I});

                fprintf(fid,'%s\n',str);

            end;

            fclose(fid);

            % HCompV just gets the vfloor stuff out so we can begin approximating

            cmd = ['"HCompV" -T 2 -D -C "' config '" -o hmmdef -f 0.01 -m -S "' LIST_FILE '" -M "' hmm_dir filesep 'hmm0" "' proto '"'];

            system(cmd);

            %cmd = [BINDIR filesep 'macro' binExt ' ' num2str(NUM_COEF) ' ' PAR_TYPE ' ' '"' hmm_dir filesep 'hmm0' filesep 'vFloors' '" "' hmm_dir filesep 'hmm0' filesep 'macros' '"'];

            %system(cmd);

            obj.macroMat(NUM_COEF,PAR_TYPE, fullfile(hmm_dir, 'hmm0', 'vFloors'), fullfile(hmm_dir, 'hmm0', 'macros'));

            %MAKE THE INITIAL MODEL PROTOTYPE

            genProto(obj);

            disp('Seed HMM successfully Produced.....');

            %Training

            for I = 1:3

                disp(I)

                j = I-1;

                cmd = ['HERest -D -C ' config ' -I ' labels ' -t 250.0 150.0 1000.0 -S ' LIST_FILE  ' -H ' hmm_dir filesep 'hmm' num2str(j)  filesep 'macros -H ' hmm_dir filesep 'hmm' num2str(j) filesep 'models -M ' hmm_dir filesep 'hmm' num2str(I) ' ' word_list];

                %                 disp(cmd)

                system(cmd);

            end

            disp('3 iterations complete');

            rmdir ([hmm_dir filesep 'hmm4'],'s')

            copyfile ([hmm_dir filesep 'hmm3'], [hmm_dir filesep 'hmm4'])

            % The following command takes state 3 from the silence model

            % and appends it to the end of the model as state 2 of the

            % short pause model.

            % Original:

            % cmd = [BINDIR filesep 'spmodel_gen' binExt ' ' hmm_dir filesep 'hmm3' filesep 'models ' hmm_dir filesep 'hmm4' filesep 'models'];

            % system(cmd);

            % New:

            obj.spmodel_genMat(fullfile(hmm_dir,'hmm3','models'), fullfile(hmm_dir,'hmm4','models'));

            cmd = ['HHEd  -T 2 -H ' hmm_dir filesep 'hmm4' filesep 'macros -H ' hmm_dir filesep 'hmm4' filesep 'models -M ' hmm_dir filesep 'hmm5 ' ED_CMDFILE1 ' ' word_listSP ];

            system(cmd);

            disp ('SP model fixed')

            % after the spmodel_gen command - the word_list is changed to

            % word_listSP. The sp model is just ignored currently

            for I = 6:8

                disp(I)

                j = I-1;

                cmd = ['HERest -C ' config ' -I ' labelssp ' -S ' LIST_FILE  ' -H ' hmm_dir filesep 'hmm' num2str(j)  filesep 'macros -H ' hmm_dir filesep 'hmm' num2str(j) filesep 'models -M ' hmm_dir filesep 'hmm' num2str(I) ' ' word_listSP];

                system(cmd);

            end

            disp('6 iterations complete');

            cmd = ['HHEd  -T 2 -H ' hmm_dir filesep 'hmm8' filesep 'macros -H ' hmm_dir filesep 'hmm8' filesep 'models -M ' hmm_dir filesep 'hmm9 ' ED_CMDFILE2 ' ' word_listSP ];

            system(cmd);

            disp ('2 gaussians per mixture')

            for I = 10:12

                disp(I)

                j = I-1;

                cmd = ['HERest -C ' config ' -I ' labelssp ' -S ' LIST_FILE  ' -H ' hmm_dir filesep 'hmm' num2str(j)  filesep 'macros -H ' hmm_dir filesep 'hmm' num2str(j) filesep 'models -M ' hmm_dir filesep 'hmm' num2str(I) ' ' word_listSP];

                system(cmd);

            end

            disp ('9 iterations completed')

            cmd = ['HHEd  -T 2 -H ' hmm_dir filesep 'hmm12'  filesep 'macros -H ' hmm_dir filesep 'hmm12' filesep 'models -M ' hmm_dir filesep 'hmm13 ' ED_CMDFILE3 ' ' word_listSP ];

            system(cmd);

            disp ('3 gaussians per mixture')

            for I = 14:20

                disp(I)

                j = I-1;

                cmd = ['HERest -C ' config ' -I ' labelssp ' -S ' LIST_FILE  ' -H ' hmm_dir filesep 'hmm' num2str(j)  filesep 'macros -H ' hmm_dir filesep 'hmm' num2str(j) filesep 'models -M ' hmm_dir filesep 'hmm' num2str(I) ' ' word_listSP];

                system(cmd);

            end

            disp ('16 iterations completed')

            cmd = ['HHEd  -T 2 -H ' hmm_dir filesep 'hmm20'  filesep 'macros -H ' hmm_dir filesep 'hmm20' filesep 'models -M ' hmm_dir filesep 'hmm21 ' ED_CMDFILE4 ' ' word_listSP ];

            system(cmd);

            disp ('5 gaussians per mixture')

            for I = 22:28

                disp(I)

                j = I-1;

                cmd = ['HERest -C ' config ' -I ' labelssp ' -S ' LIST_FILE  ' -H ' hmm_dir filesep 'hmm' num2str(j)  filesep 'macros -H ' hmm_dir filesep 'hmm' num2str(j) filesep 'models -M ' hmm_dir filesep 'hmm' num2str(I) ' ' word_listSP];

                system(cmd);

            end

            disp ('23 iterations completed')

            cmd = ['HHEd  -T 2 -H ' hmm_dir filesep 'hmm28'  filesep 'macros -H ' hmm_dir filesep 'hmm28' filesep 'models -M ' hmm_dir filesep 'hmm29 ' ED_CMDFILE5 ' ' word_listSP ];

            system(cmd);

            disp ('7 gaussians per mixture')

            for I = 30:36

                disp(I)

                j = I-1;

                cmd = ['HERest -C ' config ' -I ' labelssp ' -S ' LIST_FILE  ' -H ' hmm_dir filesep 'hmm' num2str(j)  filesep 'macros -H ' hmm_dir filesep 'hmm' num2str(j) filesep 'models -M ' hmm_dir filesep 'hmm' num2str(I) ' ' word_listSP];

                system(cmd);

            end

            disp ('30 iterations completed')

        end % ------ OF TRAIN

        %% **********************************************************

        % Test

        %************************************************************

        function test(obj, testFeatureFolder)

            flags = '-p 0.0 -s 0.0';

            test_word_list = obj.testWordListFile;

            net = obj.wordNetFile;

            dict = obj.dictFile;

            LIST_FILE       = fullfile(testFeatureFolder, 'tmp.list');

            config          = obj.configFile;

            disp ('Now testing with 7 mixture HMMs')

            mod_file = fullfile(obj.hmmFolder, 'hmm36', 'models');

            mac_file = fullfile(obj.hmmFolder, 'hmm36', 'macros');

            RESULTS_DIR = testFeatureFolder;

            TEST_FEAT_ROOT  = testFeatureFolder;

            N1list =  fullfile(testFeatureFolder, 'list.scp');

            ftest = fopen(N1list,'r');

            S = textscan(ftest,'%s','Delimiter','\n');

            fclose(ftest);

            flist = fopen(LIST_FILE,'w');

            for I = 1:size(S{1},1)

                str  = fullfile(TEST_FEAT_ROOT,S{1}{I});

                fprintf(flist,'%s\n',str);

            end;

            fclose(flist);

            cmd = ['HVite -D -H ' mac_file ' -H ' mod_file ' -S ' LIST_FILE ' -C ' config ' -w ' net ' -l ''*'' -i ' RESULTS_DIR filesep 'result.mlf ' flags ' ' dict ' ' test_word_list];

            system(cmd);

        end % ------ OF TEST

        %% **********************************************************

        % Get methods determining feature vector related gubbins

        %************************************************************

        function value = get.protoFile(obj)

            value = fullfile(obj.hmmFolder, 'proto_AutoGen');

        end

    end % ------ OF METHODS

    %%  *********************************************************

    %      _        _   _                       _   _               _

    %     | |      | | (_)                     | | | |             | |

    %  ___| |_ __ _| |_ _  ___   _ __ ___   ___| |_| |__   ___   __| |___

    % / __| __/ _` | __| |/ __| | '_ ` _ \ / _ \ __| '_ \ / _ \ / _` / __|

    % \__ \ || (_| | |_| | (__  | | | | | |  __/ |_| | | | (_) | (_| \__ \

    % |___/\__\__,_|\__|_|\___| |_| |_| |_|\___|\__|_| |_|\___/ \__,_|___/

    %************************************************************

    methods(Static)

        %% **********************************************************

        % macroMat (matlab port of macro C code)

        % ported by NC - Nov 2011

        %************************************************************

        function macroMat(VECSIZE, PARAMETER_TYPE, infile, outfile)

          % This function takes the vFloors file that is created after

          % invoking HCompV and generates a Macro file required for further

          % HMM training

            ofp = fopen(outfile,'w');

            fprintf(ofp, '~o\n<STREAMINFO> 1 %d\n', VECSIZE);

            fprintf(ofp, '<VECSIZE> %d\n', VECSIZE);

            fprintf(ofp, '<NULLD>\n<%s>\n', PARAMETER_TYPE);

            fp = fopen(infile);

            tline = fgets(fp);

            while ischar(tline)

                fprintf(ofp,tline);

                tline = fgets(fp);

            end

            fclose(fp);

            fclose(ofp);

        end %---- of MACROMAT

        %% **********************************************************

        % spmodel_genMat (matlab port of spmodel_gen C code)

        % ported by NC - Nov 2011

        %************************************************************

        function spmodel_genMat(infile, outfile)

          % This function copies the middle state (3) from the silence model

          % and makes a sp model out of it by copying it to state 2 of the sp model.

          % From: http://www.voxforge.org/home/dev/acousticmodels/linux/create/htkjulius/tutorial/monophones/step-7

            %  In the last step you created HMM models that did not include an "sp"

            %  (short pause) silence model - which refers to the types of short pauses

            %  that occur between words in normal speech.  However, you did create a

            %  "sil" silence model - sil silence models are typically of longer

            %  duration, and refer to the pauses occur at the end of a sentence.

            %

            % The HTK book says that the sp model needs to have its "emitting state

            % tied to the centre state of the silence model".  What this means is that

            % you need to create a new sp model in your hmmdefs, that it will use the

            % centre state of sil, and then they both need to be 'tied' together.  For

            % a bit of background on HMMs and states, see this example.

            %

            % This can be done by copying the centre state from the sil model in your

            % hmmdefs file and adding it to the sp model, and then running a special

            % tool called HHED to 'tie' the sp model to the sil model so that they

            % share the same centre state.  The HTK book provides some background on

            % what this means, but you need an understanding of the basics of Hidden

            % Markov Modelling before tackling the HTK Book explanations

            ofp = fopen(outfile,'a+'); % we append this time

            fprintf(ofp,'~h "sp"\n');

            fprintf(ofp,'<BEGINHMM>\n<NUMSTATES> 3\n<STATE> 2\n');

            %-- This block gets the hmmdef file to the MODEL

            lNow = [];

            fp = fopen(infile);

            while ~(strcmpi(lNow, '~h "sil"'))

                lNow = fgetl(fp);

            end

            %------------------------------

            %-- This block gets the hmmdef file to the STATE

            lNow = [];

            fp = fopen(infile);

            while ~(strcmpi(lNow, ['<STATE> ' num2str(3)]))

                lNow = fgetl(fp);

            end

            %------------------------------

            %%% This block puts a copy of the hmmdef file in from the

            %%% correct line as found above

            tline = fgetl(fp);

            while ~(strcmpi(tline, ['<STATE> ' num2str(4)]))

                fprintf(ofp,'%s\n', tline);

                tline = fgetl(fp);

            end

            fclose(fp);

            %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

            fprintf(ofp,'<TRANSP> 3\n0.000000e+00 1.000000e+00 0.000000e+00\n');

            fprintf(ofp,'0.000000e+00 5.000000e-01 5.000000e-01\n');

            fprintf(ofp,'0.000000e+00 0.000000e+00 0.000000e+00\n<ENDHMM>\n');

            fclose(ofp);

        end %---- of spmodel_genMat

        %% **********************************************************

        % models_1mixsilMat (matlab port of models_1mixsil C code)

        % ported by NC - Nov 2011

        %************************************************************

        function models_1mixsilMat(infile, outfile)

            % This function takes the hmmdef file and

            % generates a HMM Model file

            ofp = fopen(outfile,'w');

            for ii = 1:11

                %%% This block gets the hmmdef file to the correct line

                lNow = [];

                fp = fopen(infile);

                while ~(strcmpi(lNow, '~h "hmmdef"'))

                    lNow = fgetl(fp);

                end

                %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

                %%% This block puts a numbered header

                if ii==1;  fprintf(ofp,'~h "one"\n'); end

                if ii==2;  fprintf(ofp,'~h "two"\n'); end

                if ii==3;  fprintf(ofp,'~h "three"\n'); end

                if ii==4;  fprintf(ofp,'~h "four"\n'); end

                if ii==5;  fprintf(ofp,'~h "five"\n'); end

                if ii==6;  fprintf(ofp,'~h "six"\n'); end

                if ii==7;  fprintf(ofp,'~h "seven"\n'); end

                if ii==8;  fprintf(ofp,'~h "eight"\n'); end

                if ii==9;  fprintf(ofp,'~h "nine"\n'); end

                if ii==10; fprintf(ofp,'~h "oh"\n'); end

                if ii==11; fprintf(ofp,'~h "zero"\n'); end

                %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

                %%% This block puts a copy of the hmmdef file in from the

                %%% correct line as found above

                tline = fgets(fp);

                while ischar(tline)

                    fprintf(ofp,tline);

                    tline = fgets(fp);

                end

                fclose(fp); %close it as we reached EOF

                %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

            end

            %%% THIS IS THE SILENCE MODEL @ THE END

            fprintf(ofp,'~h "sil"\n');

            fprintf(ofp,'<BEGINHMM>\n<NUMSTATES> 5\n');

            for kk = 2:4

                fprintf(ofp,'<STATE> %d\n<NUMMIXES> 1\n', kk);

                %-- This block gets the hmmdef file to the correct line

                lNow = [];

                fp = fopen(infile);

                while ~(strcmpi(lNow, ['<STATE> ' num2str(kk)]))

                    lNow = fgetl(fp);

                end

                %------------------------------

                %%% This block puts a copy of the hmmdef file in from the

                %%% correct line as found above

                tline = fgetl(fp);

                while ~(strcmpi(tline, ['<STATE> ' num2str(kk+1)]))

                    fprintf(ofp,'%s\n', tline);

                    tline = fgetl(fp);

                end

                fclose(fp);

                %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

            end

            fprintf(ofp,'<TRANSP> 5\n0.000000e+00 1.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00\n');

            fprintf(ofp,'0.000000e+00 6.000000e-01 4.000000e-01 0.000000e+00 0.000000e+00\n');

            fprintf(ofp,'0.000000e+00 0.000000e+00 6.000000e-01 4.000000e-01 0.000000e+00\n');

            fprintf(ofp,'0.000000e+00 0.000000e+00 0.000000e+00 7.000000e-01 3.000000e-01\n');

            fprintf(ofp,'0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00\n<ENDHMM>\n');

            fclose(ofp);

        end %---- of models_1mixsilMat

        %% **********************************************************

        % makeProtoHmm (Make a prototype HMM)

        %************************************************************

        function makeProtoHmm(filename,featureType,numFeatures,numStates)

            % filename is obvious

            % featureType is usually USER_D_A

            % numFeatures is number of features (including differences if used)

            % numStates is usually 18

            ofp = fopen(filename,'w');

            fprintf(ofp,'<BeginHMM>\n');

            fprintf(ofp,' <NumStates> %d <VecSize> %d <%s> <nullD> <diagC>\n',numStates,numFeatures,featureType');

            fprintf(ofp,' <StreamInfo> 1 %d\n',numFeatures);

            for state=2:numStates-1,

                fprintf(ofp,' <State> %d <NumMixes> 1\n',state);

                fprintf(ofp,'  <Stream> 1\n');

                fprintf(ofp,'  <Mixture> 1 1.0\n');

                fprintf(ofp,'    <Mean> %d\n',numFeatures);

                fprintf(ofp,'      ');

                fprintf(ofp,'%1.1f ',zeros(1,numFeatures));

                fprintf(ofp,'\n');

                fprintf(ofp,'   <Variance> %d\n',numFeatures);

                fprintf(ofp,'      ');

                fprintf(ofp,'%1.1f ',ones(1,numFeatures));

                fprintf(ofp,'\n');

            end

            fprintf(ofp,'<TransP> %d\n',numStates);

            transp = zeros(numFeatures);

            transp(1,2)=1;

            for state=2:numStates-2,

                transp(state,state)=0.6;

                transp(state,state+1)=0.4;

            end

            transp(numStates-1,numStates-1)=0.9;

            transp(numStates-1,numStates)=0.1;

            for state=1:numStates,

                fprintf(ofp,'%1.3e  ',transp(state,1:numStates));

                fprintf(ofp,'\n');

            end

            fprintf(ofp,'<EndHMM>\n');

            fclose(ofp);

        end %---- of MAKEPROTOHMM

        %% **********************************************************

        % createMLF - master label file - belongs with hmm class

        %************************************************************

        function createMLF(mapFileFolder)

            mlfFileName = 'labels';

            d = dir(fullfile(mapFileFolder, '*.map'));

            fid = fopen([fullfile(mapFileFolder,mlfFileName) '.mlf'],'w');

            fprintf(fid,'#!MLF!#\n');

            for I = 1:size(d,1)

                fprintf(fid,['"*/' d(I).name(1:end-3) 'lab"\n']);

                fprintf(fid,'sil\n');

                labels = d(I).name(5:end-5);

                for J = 1:length(labels);

                    switch labels(J)

                        case 'O'

                            S = 'oh';

                        case '1'

                            S = 'one';

                        case '2'

                            S = 'two';

                        case '3'

                            S = 'three';

                        case '4'

                            S = 'four';

                        case '5'

                            S = 'five';

                        case '6'

                            S = 'six';

                        case '7'

                            S = 'seven';

                        case '8'

                            S = 'eight';

                        case '9'

                            S = 'nine';

                        case 'Z'

                            S = 'zero';

                    end

                    fprintf(fid,S);

                    fprintf(fid,'\n');

                end;

                fprintf(fid,'sil\n.\n');

            end;

            fclose(fid);

        end % ------ OF CREATEMLF

        %% **********************************************************

        % createSCP - Just a file list

        %************************************************************

        function createSCP(mapFileFolder)

            scpFileName = 'list';

            d = dir(fullfile(mapFileFolder, '*.map'));

            fid = fopen([fullfile(mapFileFolder,scpFileName) '.scp'],'w');

            for I = 1:size(d,1)

                fprintf(fid,[d(I).name(1:end-3) 'map\n']);

            end;

            fclose(fid);

        end% ------ OF createSCP

        %% **********************************************************

        % score - Simple non-dynamic scoring

        %************************************************************

        function score(testMLFpath)

            % First get the test result files and store the file locations

            testIdx = 1;

            temp = dir(testMLFpath);

            for nn = 1:numel(temp)

                if length(temp(nn).name) > 3

                    if strcmp(temp(nn).name(end-3:end), '.mlf')

                        testMLFfiles{testIdx} = temp(nn).name;

                        testIdx = testIdx+1;

                    end

                end

            end

            for fileIdx = 1:numel(testMLFfiles)

                % Extract file name and regonised words as strings -> store in recResults.

                % This code block is Matt Robertson's

                fmlf = fopen(fullfile(testMLFpath , testMLFfiles{fileIdx}),'r');

                MLF = textscan(fmlf,'%s','delimiter','\n');

                replaceDownPattern = '[0-9\-\.\s]';

                for I  = 1:size(MLF{1},1);

                    if strfind(MLF{1}{I},'.rec') >0

                        M{I} = MLF{1}{I};

                    else

                        if strfind(MLF{1}{I},'MLF') >0;

                            M{I} = MLF{1}{I};

                        else

                            M{I} = regexprep(MLF{1}{I},replaceDownPattern,'');

                        end

                    end;

                end;

                fclose(fmlf);

                A = 1;

                for I = 1:size(M,2)

                    if not(strcmpi(M{I},'sp')) && not(strcmpi(M{I},'sil')) && not(strcmpi(M{I},''));

                        M2{A} = M{I};

                        A = A + 1;

                    end;

                end;

                recResults = [];

                X=0;

                for I = 1:length(M2)

                    if strfind(M2{I},'MLF')

                    else

                        if strfind(M2{I},'.rec')

                            X = X+1;

                            % found FileName

                            recResults(X).fileName = M2{I};

                            recResults(X).utterance = {};

                        else

                            recResults(X).utterance = [recResults(X).utterance M2{I}];

                        end;

                    end;

                end;

                % Make separate 2D arrays of input and output digits

                inputDigits  = zeros(numel(recResults),3);

                outputDigits = zeros(numel(recResults),3);

                for nn = 1:numel(recResults)

                    ipStr = recResults(nn).fileName(end-8:end-6);

                    for kk = 1:3

                        inputDigits(nn,kk)  = cHMM.htk_str2num(ipStr(kk)); %see local function @ bottom of script

                        opStr = recResults(nn).utterance(kk);

                        outputDigits(nn,kk) = cHMM.htk_str2num(opStr{:}); %bit of a hack to do cell2str as it were

                    end

                end

                % Now do the scoring (simple whan data in the right format)

                % :D

                scoreArray = (inputDigits == outputDigits);

                pcSent(fileIdx) = 100*sum(all(scoreArray,2)) / numel(recResults);

                pcWord(fileIdx) = 100*sum(scoreArray(:)) / numel(scoreArray);

            end

            %will output table to console if used on *nix

            xlsdataFull = [{'-- File Name --', '%Sent', '%Word'}; testMLFfiles' num2cell(pcSent') num2cell(pcWord')];

            disp(xlsdataFull)

            fid = fopen(fullfile(testMLFpath,['score__' num2str(pcWord(1), '%0.1f') '__.txt']),'w');

            fclose(fid);

        end% ------ OF SCORE

        %% **********************************************************

        % scoreWhole folder - make my life easier @ command line

        %************************************************************

        function scoreWholeFolder(folderToScore, searchString)

            if nargin < 2

                searchString = '*featR*';

            end

            dirInfo = dir(fullfile(folderToScore, searchString));

            numFolders = numel(dirInfo);

            for nn = 1:numFolders;

                currentScoring = fullfile(folderToScore, dirInfo(nn).name);

                disp('');

                disp(currentScoring);

                cHMM.score(currentScoring);

            end

        end

        %**************************************************************************

        % htk_str2num - Convert strings to integers

        %**************************************************************************

        function opNum = htk_str2num(ipString)

            switch ipString

                case {'oh' , 'O'}

                    opNum = 0;

                case {'one' , '1'}

                    opNum = 1;

                case {'two' , '2'}

                    opNum = 2;

                case {'three' , '3'}

                    opNum = 3;

                case {'four' , '4'}

                    opNum = 4;

                case {'five' , '5'}

                    opNum = 5;

                case {'six' , '6'}

                    opNum = 6;

                case {'eight' , '8'}

                    opNum = 8;

                case {'nine' , '9'}

                    opNum = 9;

                otherwise

                    assert(0) % throw error

            end

        end % ------ OF htk_str2num

    end % ------ OF STATIC METHODS

end % ------ OF CLASS