Ejemplo n.º 1
0
void Neuron::calcError() {
    if (!error_recalc) return;
    error_recalc = false;
    calcEnergy();
    double sum = 0;
    for (Link * i: outputs) {
        sum += i->getError();
    }
    error = ::activation_function_derivative(energy) * sum;
}
Ejemplo n.º 2
0
void outputData(Tsph & sph) {
    if(RP::Time >= RP::TimeAscii) {
        char filename[64];
        sprintf(filename, "snap/sph_t%04d.dat", RP::NumberOfAscii);
        sph.writeParticleAscii(filename);
        RP::TimeAscii += RP::TimestepAscii;
        RP::NumberOfAscii++;        
    }
    PS::F64    etot = calcEnergy(sph);
    PS::F64vec vtot = calcMomentum(sph);
    if(PS::Comm::getRank() == 0) {
        fprintf(RP::FilePointerForLog, "time: %.10f %+e %+e %+e %+e\n",
                RP::Time, RP::Timestep, etot, vtot[0], WT::getTimeTotal());
        WT::dump(RP::Time, RP::FilePointerForTime);
        fflush(RP::FilePointerForLog);
    }
}
Ejemplo n.º 3
0
void fill_stdhep(int i)
{
    int j, k, num;

    num = i * 20;
    hepevt_.nevhep = i;
    hepevt_.nhep = num;
    for (j = 0; j < num; j++) {
        hepevt_.idhep[j] = 211;
        hepevt_.isthep[j] = 1;
        for (k = 0; k < 2; k++) {
            hepevt_.jmohep[j][k] = j - 1;
            hepevt_.jdahep[j][k] = j + 1;
            }
        for (k = 0; k < 3; k++)
            hepevt_.phep[j][k] = 10.2 * k + 0.31 * i;
	hepevt_.phep[j][4] = .511 + i * 0.04;
	hepevt_.phep[j][3] = calcEnergy(j);
        for (k = 0; k < 4; k++)
            hepevt_.vhep[j][k] = .0001;
        }
}
Ejemplo n.º 4
0
int main(int argc, char *argv[]) {
    PS::F32 time  = 0.0;
    PS::F32 tend  = 10.0;
    PS::F32 dtime = 1.0 / 128.0;
    PS::F32 dtout = 1.0 / 8.0;
    PS::S64 ntot  = 1024;
    
    PS::Initialize(argc, argv);
    
    PS::DomainInfo dinfo;
    dinfo.initialize();

    PS::ParticleSystem<FPGrav> system_grav;
    system_grav.initialize();

    PS::TreeForForceLong<FPGrav, FPGrav, FPGrav>::
        Monopole tree_grav;
    tree_grav.initialize(ntot);

    setParticleColdUniformSphere(system_grav, ntot);

    dinfo.decomposeDomainAll(system_grav);

    system_grav.exchangeParticle(dinfo);
    
    tree_grav.calcForceAllAndWriteBack
        (CalcGravity<FPGrav>(),
         CalcGravity<PS::SPJMonopole>(),
         system_grav,
         dinfo);
    
    PS::F64 etot0 = calcEnergy(system_grav);
    if(PS::Comm::getRank() == 0) {
        fprintf(stderr, "time: %10.7f energy: %+e energy error: %+e\n",
                time, etot0, (etot0 - etot0) / etot0);
    }

    while(time < tend) {

        predict(system_grav, dtime);        
        dinfo.decomposeDomainAll(system_grav);
        system_grav.exchangeParticle(dinfo);        
        tree_grav.calcForceAllAndWriteBack
            (CalcGravity<FPGrav>(),
             CalcGravity<PS::SPJMonopole>(),
             system_grav,
             dinfo);
        correct(system_grav, dtime);
        
        time += dtime;
        PS::F64 etot1 = calcEnergy(system_grav);
        if(fmod(time, dtout) == 0.0 &&
           PS::Comm::getRank() == 0) {
                fprintf(stderr, "time: %10.7f energy: %+e energy error: %+e\n",
                        time, etot1, (etot1 - etot0) / etot0);
        }

    }

    PS::Finalize();

    return 0;
}
Ejemplo n.º 5
0
double Neuron::getSignal() {
    calcEnergy();
    return energy;
}