Use test_code3.cpp to get a working reaction divider