newmen/versatile-diamond

#include "two_bridges_to_high_bridge.h"

const char TwoBridgesToHighBridge::__name[] = "two bridges to high bridge";

double TwoBridgesToHighBridge::RATE()
{
    static double value = getRate("TWO_BRIDGES_TO_HIGH_BRIDGE");
    return value;
}

void TwoBridgesToHighBridge::find(TwoBridgesCTRiCBRs *target)
{
    create<TwoBridgesToHighBridge>(target);
}

void TwoBridgesToHighBridge::doIt()
{
    assert(target()->type() == TwoBridgesCTRiCBRs::ID);

    Atom *atoms[3] = { target()->atom(3), target()->atom(4), target()->atom(5) };
    Atom *a = atoms[0], *b = atoms[1], *c = atoms[2];

    assert(a->is(0));
    assert(b->is(5));
    assert(c->is(24));

    // erase from crystal should be before bond-unbond atoms
    a->eraseFromCrystal();
    Handbook::amorph().insert(a);

    a->unbondFrom(c);
    a->bondWith(b);

    if (a->is(2)) a->changeType(17);
    else if (a->is(1)) a->changeType(16);
    else a->changeType(15);

    b->changeType(19);
    c->changeType(5);

    Finder::findAll(atoms, 3);
}