double All_hadrons_triple_integral::get_triple_integral_value(double time)
{
value = 0.0;
Hadron_triple_integral hti(params);
for(std::map<Hadron, unsigned, Hadron_comparator>::iterator it = map_of_different_hadrons.begin() ; it != map_of_different_hadrons.end() ; ++it )
{
double val = hti.get_triple_integral_value(time, it->first );
std::cout << it->second << " " << val << "\n";
value += it->second * val;
}
return value;
}
Z void growSymHashTable(void)
{
I i;
S s0,s1,*sp;
HT newtable=hti(4*SymHashTable->nb);
for(i=0;i<SymHashTable->nb;++i)
for(s0=(S)SymHashTable->b[i];s0;s0=s1)
{
s1=s0->s;
sp=(S *)(newtable->b+((newtable->nb-1)&HA(s0->n)));
s0->s=*sp;
*sp=s0;
++newtable->ni;
}
bfree((char *)SymHashTable);
SymHashTable=newtable;
}
void symhti(void){SymHashTable=hti(STARTHASHSIZE);}
