alternaria_moves.h

Go to the documentation of this file.

/*
 * This work is licensed under a Creative Commons 
 * Attribution-Noncommercial-Share Alike 3.0 United States License.
 * 
 *    http://creativecommons.org/licenses/by-nc-sa/3.0/us/
 * 
 * You are free:
 * 
 *    to Share - to copy, distribute, display, and perform the work
 *    to Remix - to make derivative works
 * 
 * Under the following conditions:
 * 
 *    Attribution. You must attribute the work in the manner specified by the
 *    author or licensor (but not in any way that suggests that they endorse you
 *    or your use of the work).
 * 
 *    Noncommercial. You may not use this work for commercial purposes.
 * 
 *    Share Alike. If you alter, transform, or build upon this work, you may
 *    distribute the resulting work only under the same or similar license to
 *    this one.
 * 
 * For any reuse or distribution, you must make clear to others the license
 * terms of this work. The best way to do this is by including this header.
 * 
 * Any of the above conditions can be waived if you get permission from the
 * copyright holder.
 * 
 * Apart from the remix rights granted under this license, nothing in this
 * license impairs or restricts the author's moral rights.
 */


#ifndef ALTERNARIA_MOVES_H
#define ALTERNARIA_MOVES_H


#include <config.h>

#include <inttypes.h>

#include <jwsc++/base/exception.h>

#include "sampler.h"


class Apical_hypha_birth_death_move : public Sampler_jump_move
{
    public:

        Apical_hypha_birth_death_move
        (
            const char* birth_name, 
            const char* death_name,
            double      birth_prob=1.0,
            double      death_prob=1.0
        ) 
        throw (jwscxx::base::Arg_error) 
        : Sampler_jump_move(birth_name, death_name, birth_prob, death_prob) { }


        bool run_1(Sampler_move_parameters* params)
            throw (jwscxx::base::Exception);


        bool run_2(Sampler_move_parameters* params)
            throw (jwscxx::base::Exception);
};


class Lateral_hypha_birth_death_move : public Sampler_jump_move
{
    public:

        Lateral_hypha_birth_death_move
        (
            const char* birth_name, 
            const char* death_name,
            double      birth_prob=1.0,
            double      death_prob=1.0
        ) 
        throw (jwscxx::base::Arg_error) 
        : Sampler_jump_move(birth_name, death_name, birth_prob, death_prob) { }


        bool run_1(Sampler_move_parameters* params)
            throw (jwscxx::base::Exception);


        bool run_2(Sampler_move_parameters* params)
            throw (jwscxx::base::Exception);
};


class Spore_birth_death_move : public Sampler_jump_move
{
    public:

        Spore_birth_death_move
        (
            const char* birth_name, 
            const char* death_name,
            double      birth_prob=1.0,
            double      death_prob=1.0
        ) 
        throw (jwscxx::base::Arg_error) 
        : Sampler_jump_move(birth_name, death_name, birth_prob, death_prob) { }


        bool run_1(Sampler_move_parameters* params)
            throw (jwscxx::base::Exception);


        bool run_2(Sampler_move_parameters* params)
            throw (jwscxx::base::Exception);
};


class Spore_transition_apical_move : public Sampler_jump_move
{
    public:

        Spore_transition_apical_move
        (
            const char* apical_to_spore_name, 
            const char* spore_to_apical_name,
            double      apical_to_spore_prob=1.0,
            double      spore_to_apical_prob=1.0
        ) 
        throw (jwscxx::base::Arg_error) 
        : Sampler_jump_move(apical_to_spore_name, spore_to_apical_name, 
                apical_to_spore_prob, spore_to_apical_prob) { }


        bool run_1(Sampler_move_parameters* params)
            throw (jwscxx::base::Exception);


        bool run_2(Sampler_move_parameters* params)
            throw (jwscxx::base::Exception);
};


class Spore_split_merge_2_apical_move : public Sampler_jump_move
{
    public:

        Spore_split_merge_2_apical_move
        (
            const char* split_name, 
            const char* merge_name,
            double      split_prob=1.0,
            double      merge_prob=1.0
        ) 
        throw (jwscxx::base::Arg_error) 
        : Sampler_jump_move(split_name, merge_name, split_prob, merge_prob) { }


        bool run_1(Sampler_move_parameters* params)
            throw (jwscxx::base::Exception);


        bool run_2(Sampler_move_parameters* params)
            throw (jwscxx::base::Exception);


    protected:

        static void propose_1st_split
        (
            Apical_hypha**              proposal_out,
            double*                     log_prob_out,
            const Apical_hypha_density* density,
            const Spore*                spore
        )
        throw (jwscxx::base::Arg_error);


        static void propose_2nd_split
        (
            Apical_hypha**      proposal_out,
            double*             log_prob_out,
            const Apical_hypha* split_1
        )
        throw (jwscxx::base::Arg_error);

        static void propose_merge
        (
            Spore**              proposal_out,
            double*              log_prob_out,
            const Structure*     parent,
            const Spore_density* density,
            const Apical_hypha*  hypha_1,
            const Apical_hypha*  hypha_2
        )
        throw (jwscxx::base::Arg_error);
};




class Spore_split_merge_spore_move : public Sampler_jump_move
{
    public:

        Spore_split_merge_spore_move
        (
            const char* split_name, 
            const char* merge_name,
            double      split_prob=1.0,
            double      merge_prob=1.0
        ) 
        throw (jwscxx::base::Arg_error) 
        : Sampler_jump_move(split_name, merge_name, split_prob, merge_prob) { }


        bool run_1(Sampler_move_parameters* params)
            throw (jwscxx::base::Exception);


        bool run_2(Sampler_move_parameters* params)
            throw (jwscxx::base::Exception);


    protected:

        static void propose_1st_split
        (
            Spore**      proposal_out,
            double*      log_prob_out,
            const Spore* spore
        )
        throw (jwscxx::base::Arg_error);


        static void propose_2nd_split
        (
            Spore**      proposal_out,
            double*      log_prob_out,
            const Spore* split_1
        )
        throw (jwscxx::base::Arg_error);

        static void propose_merge
        (
            Spore**          proposal_out,
            double*          log_prob_out,
            const Structure* parent,
            const Spore*     hypha_1,
            const Spore*     hypha_2
        )
        throw (jwscxx::base::Arg_error);
};




// TODO Spore_split_merge_apical_move


// TODO Apical_hypha_split_merge_spore_move


class Apical_hypha_split_merge_apical_move : public Sampler_jump_move
{
    public:

        Apical_hypha_split_merge_apical_move
        (
            const char* split_name, 
            const char* merge_name,
            double      split_prob=1.0,
            double      merge_prob=1.0
        ) 
        throw (jwscxx::base::Arg_error) 
        : Sampler_jump_move(split_name, merge_name, split_prob, merge_prob) { }


        bool run_1(Sampler_move_parameters* params)
            throw (jwscxx::base::Exception);


        bool run_2(Sampler_move_parameters* params)
            throw (jwscxx::base::Exception);


    protected:

        static void propose_1st_split
        (
            Apical_hypha**      proposal_out,
            double*             log_prob_out,
            const Apical_hypha* hypha
        )
        throw (jwscxx::base::Arg_error);


        static void propose_2nd_split
        (
            Apical_hypha**      proposal_out,
            double*             log_prob_out,
            const Apical_hypha* split_1
        )
        throw (jwscxx::base::Arg_error);

        static void propose_merge
        (
            Apical_hypha**      proposal_out,
            double*             log_prob_out,
            const Structure*    parent,
            const Apical_hypha* hypha_1,
            const Apical_hypha* hypha_2
        )
        throw (jwscxx::base::Arg_error);
};


class Apical_hypha_split_merge_lateral_move : public Sampler_jump_move
{
    public:

        Apical_hypha_split_merge_lateral_move
        (
            const char* split_name, 
            const char* merge_name,
            double      split_prob=1.0,
            double      merge_prob=1.0
        ) 
        throw (jwscxx::base::Arg_error) 
        : Sampler_jump_move(split_name, merge_name, split_prob, merge_prob) { }


        bool run_1(Sampler_move_parameters* params)
            throw (jwscxx::base::Exception);


        bool run_2(Sampler_move_parameters* params)
            throw (jwscxx::base::Exception);


    protected:

        static void propose_1st_split
        (
            Apical_hypha**      proposal_out,
            double*             log_prob_out,
            const Apical_hypha* hypha
        )
        throw (jwscxx::base::Arg_error);


        static void propose_2nd_split
        (
            Lateral_hypha**              proposal_out,
            double*                      log_prob_out,
            const Apical_hypha*          split_1,
            const Lateral_hypha_density* density
        );

        static void propose_merge
        (
            Apical_hypha**       proposal_out,
            double*              log_prob_out,
            const Structure*     parent,
            const Apical_hypha*  hypha_1,
            const Lateral_hypha* hypha_2
        )
        throw (jwscxx::base::Arg_error);
};


class Structure_resize_move : public Sampler_diffusion_move
{
    public:

        Structure_resize_move(const char* name, double prob=1.0) 
            throw (jwscxx::base::Arg_error) 
            : Sampler_diffusion_move(name, prob) { }


        bool run(Sampler_move_parameters* params) 
            throw (jwscxx::base::Exception);


        const std::ostringstream& get_info(Sampler_move_parameters* params);


        const std::ostringstream& get_proposal_info
        (
            Sampler_move_parameters* params
        );
};


class Structure_rotate_move : public Sampler_diffusion_move
{
    public:

        Structure_rotate_move(const char* name, double prob=1.0) 
            throw (jwscxx::base::Arg_error) 
            : Sampler_diffusion_move(name, prob) { }


        bool run(Sampler_move_parameters* params) 
            throw (jwscxx::base::Exception);


        const std::ostringstream& get_info(Sampler_move_parameters* params);


        const std::ostringstream& get_proposal_info
        (
            Sampler_move_parameters* params
        );
};


class Structure_opacity_move : public Sampler_diffusion_move
{
    public:

        Structure_opacity_move(const char* name, double prob=1.0) 
            throw (jwscxx::base::Arg_error) 
            : Sampler_diffusion_move(name, prob) { }


        bool run(Sampler_move_parameters* params) 
            throw (jwscxx::base::Exception);


        const std::ostringstream& get_info(Sampler_move_parameters* params);


        const std::ostringstream& get_proposal_info
        (
            Sampler_move_parameters* params
        );
};


class Apical_structure_position_move : public Sampler_diffusion_move
{
    public:

        Apical_structure_position_move(const char* name, double prob=1.0) 
            throw (jwscxx::base::Arg_error) 
            : Sampler_diffusion_move(name, prob) { }


        bool run(Sampler_move_parameters* params) 
            throw (jwscxx::base::Exception);


        const std::ostringstream& get_info(Sampler_move_parameters* params);


        const std::ostringstream& get_proposal_info
        (
            Sampler_move_parameters* params
        );
};


class Lateral_structure_lat_dist_move : public Sampler_diffusion_move
{
    public:

        Lateral_structure_lat_dist_move(const char* name, double prob=1.0) 
            throw (jwscxx::base::Arg_error) 
            : Sampler_diffusion_move(name, prob) { }


        bool run(Sampler_move_parameters* params) 
            throw (jwscxx::base::Exception);


        const std::ostringstream& get_info(Sampler_move_parameters* params);


        const std::ostringstream& get_proposal_info
        (
            Sampler_move_parameters* params
        );
};


#endif