kinunit.cpp

Go to the documentation of this file.

#include <../../nmodlconf.h>
#include "model.h"
#include "symbol.h"
#include "units.h"
#include "parse1.hpp"
extern char *indepunits;
static List *reactnames;

static void set_flux_units(unit*);
static void react_unit_err(char*, unit*);

void kinunits(Item* type, int pass)
{
    struct unit ux1, ux2, ur1, ur2, uflux;
    Item *q;
    char *s;
    
    if (type->itemsubtype == REACT1) {
        ucopypop(&ux2);
    }
    ucopypop(&ux1);
    if (type->itemsubtype == REACT1) {
        ucopypop(&ur2);
    }
    ucopypop(&ur1);
    if (pass == 0) {
        return;
    }
    
    /* case reaction  with no numbers */
    /* should check that all states have same units */
    Unit_push((char *)0); /* this will go to any units */
    ITERATE(q, reactnames) {
        unit_push(ITM(q));
        s = (char *)(ITMA(SYM(ITM(q))->info)[6]);
        if (s) {
            Unit_push(s);
            unit_mul();
        }
        if (unit_diff()) {
Fprintf(stderr, "REACTION quantity units for %s is: %s\n", SYM(ITM(q))->name, unit_str());
            unit_pop();
Fprintf(stderr, "but the quantity units of the first term is: %s\n", unit_str());
            diag("Inconsistent material quantity units\n",
              "Need a correct COMPARTMENT statement");
        }
    }
        
    Unit_push(indepunits);
    unit_div();
    ucopypop(&uflux);
    set_flux_units(&uflux);
    
    if (type->itemsubtype == REACT1 || type->itemsubtype == '-') {
        ucopypush(&ux1);
        ucopypush(&uflux);
        ucopypush(&ur1);
        unit_div();
        if (unit_diff()) {
            react_unit_err("forward", &uflux);
        }
        unit_pop();
    }

    if (type->itemsubtype == REACT1) {
        ucopypush(&ux2);
        ucopypush(&uflux);
        ucopypush(&ur2);
        unit_div();
        if (unit_diff()) {
            react_unit_err("backward", &uflux);
        }
        unit_pop();
    }

    if (type->itemsubtype == LT) {
        ucopypush(&ux1);
        ucopypush(&uflux);
        if (unit_diff()) {
Fprintf(stderr, "Flux units are: %s\n", Unit_str(&uflux));
Fprintf(stderr, "But users << flux units are:%s\n", Unit_str(&ux1));
diag("Inconsistent flux units", (char *)0);
        }
        unit_pop();
    }
    
    freelist(&reactnames);
}

static void set_flux_units(unit* up)
{
    Symbol *s;

    Sprintf(buf, "%s", Unit_str(up));
    if ((s = lookup("f_flux")) == SYM0) {
        s = install("f_flux", NAME);
    }
    s->u.str = stralloc(buf, (char *)0);
    if ((s = lookup("b_flux")) == SYM0) {
        s = install("b_flux", NAME);
    }
    s->u.str = stralloc(buf, (char *)0);
    
}

static void react_unit_err(char* s, unit* up)
{
    
    Fprintf(stderr, "Flux units for this reaction: %s\n", Unit_str(up));
    ucopypop(up);
    Fprintf(stderr, "This implies %s rate units: %s\n", s, Unit_str(up));
    ucopypop(up);
Fprintf(stderr, "But the users %s rate units are: %s\n", s, Unit_str(up));
    diag("inconsistent reaction units", (char *)0);
}

void clear_compartlist() {
    SYMITER_STAT {
        ITMA(s->info)[6] = ITEM0;
    }}
}

void unit_compartlist(Item* q)
{
    char *ustr;
    
    ustr = (char *)(ITMA(SYM(q)->info)[6]);
    if (ustr) {
        diag(SYM(q)->name, "already in previous COMPARTMENT");
    }
    ITMA(SYM(q)->info)[6] = (Item *) stralloc(unit_str(), (char *)0);
}

void unit_ldifuslist(Item* q, int flag)
{
    char *ustr;
    unitonflag = flag;
    ustr = (char *)(ITMA(SYM(q)->info)[6]);
    if (!ustr) {
        diag(SYM(q)->name, "not declared in previous COMPARTMENT");
    }
    Unit_push("micron4/ms");
    if (!unit_cmp_exact()) {
        unit_pop();
diag(unit_str(), ": relevant area * diffusion constant must\n   be micron2 micron2/ms (1-21 m4/s)");
    }
    unit_pop();
    Unit_push("micron2");
    Unit_push(ustr);
    if (!unit_cmp_exact()) {
        diag(ustr, ": With LONGDITUDINAL_DIFFUSION the compartment \
volume\nmust be measured in micron3/micron (1-12 m2)");
    }
    unit_pop();
    unit_pop();
    unitonflag = 0;
}

void consreact_push(Item* q)
{
    char* ustr;
    unit_push(q);
    if (SYM(q)->info) {
        ustr = (char *)(ITMA(SYM(q)->info)[6]);
        if (ustr) {
            Unit_push(ustr);
            unit_mul();
        }
    }
}

void ureactadd(Item* q)
{
    if (!reactnames) {
        reactnames = newlist();
    }
    Lappenditem(reactnames, q);
}