@@ -475,27 +475,18 @@ int main(int ac, char* av[])
475
475
double primetime = 3 *tau_hom;
476
476
477
477
478
- // BEGIN Global stuff
479
- mpi::environment env (ac, av);
480
- mpi::communicator world;
481
- communicator = &world;
478
+ sprintf (strbuf, " %s/%s_e%.2et%.2f%s" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str ());
479
+ std::string logfile_prefix = strbuf;
482
480
483
- sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.log" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), world.rank ());
484
- std::string logfile = strbuf;
481
+ auryn_init (ac, av, dir=" ." , " " , logfile_prefix);
485
482
486
- LogMessageType log_level_file = PROGRESS;
487
- if ( verbose ) log_level_file = EVERYTHING;
488
- logger = new Logger (logfile,world.rank (),PROGRESS,log_level_file);
483
+ logger->set_logfile_loglevel ( PROGRESS );
484
+ if ( verbose ) logger->set_logfile_loglevel ( EVERYTHING );
489
485
490
- sys = new System (&world);
491
- // boost::filesystem::path p = av[0];
492
- // string binaryname = p.stem().string();
493
- // sys->set_simulation_name(binaryname);
494
- // END Global stuff
495
486
496
487
if (!infilename.empty ()) {
497
488
std::stringstream iss;
498
- iss << infilename << " ." << world. rank ();
489
+ iss << infilename << " ." << sys-> mpi_rank ();
499
490
infilename = iss.str ();
500
491
}
501
492
@@ -615,12 +606,12 @@ int main(int ac, char* av[])
615
606
logger->msg (msg,PROGRESS,true );
616
607
617
608
if (wmatdump) {
618
- sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.weight" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), world. rank ());
609
+ sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.weight" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), sys-> mpi_rank ());
619
610
WeightMatrixMonitor * wmatmon = new WeightMatrixMonitor ( con_ee, strbuf , wmat_interval );
620
611
}
621
612
622
613
if ( !fast ) {
623
- sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.syn" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), world. rank ());
614
+ sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.syn" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), sys-> mpi_rank ());
624
615
WeightMonitor * wmon = new WeightMonitor ( con_ee, strbuf, 10 );
625
616
for ( int i = 0 ; i < 5 ; ++i ) {
626
617
for ( int j = 0 ; j < 5 ; ++j ) {
@@ -630,28 +621,28 @@ int main(int ac, char* av[])
630
621
}
631
622
}
632
623
633
- sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.%c.ras" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), world. rank (), ' e' );
624
+ sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.%c.ras" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), sys-> mpi_rank (), ' e' );
634
625
SpikeMonitor * smon_e = new SpikeMonitor ( neurons_e, strbuf , 2500 );
635
626
636
- // sprintf(strbuf, "%s/%s_e%.2et%.2f%s.%d.%c.ras", dir.c_str(), file_prefix, eta, tau_hom, label.c_str(), world.rank (), 'p');
627
+ // sprintf(strbuf, "%s/%s_e%.2et%.2f%s.%d.%c.ras", dir.c_str(), file_prefix, eta, tau_hom, label.c_str(), sys->mpi_rank (), 'p');
637
628
// SpikeMonitor * smon_p = new SpikeMonitor( poisson, strbuf , 2500);
638
629
}
639
630
640
- sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.%c.prate" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), world. rank (), ' e' );
631
+ sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.%c.prate" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), sys-> mpi_rank (), ' e' );
641
632
PopulationRateMonitor * pmon_e = new PopulationRateMonitor ( neurons_e, strbuf, 1.0 );
642
633
643
634
if ( wall ) {
644
- sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.rt" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), world. rank ());
635
+ sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.rt" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), sys-> mpi_rank ());
645
636
RealTimeMonitor * rtmon = new RealTimeMonitor ( strbuf );
646
637
}
647
638
648
- // sprintf(strbuf, "%s/%s_e%.2et%.2f%s.%d.%c.prate", dir.c_str(), file_prefix, eta, tau_hom, label.c_str(), world.rank (), 'i');
639
+ // sprintf(strbuf, "%s/%s_e%.2et%.2f%s.%d.%c.prate", dir.c_str(), file_prefix, eta, tau_hom, label.c_str(), sys->mpi_rank (), 'i');
649
640
// PopulationRateMonitor * pmon_i = new PopulationRateMonitor( neurons_i, strbuf, 1.0 );
650
641
651
- // sprintf(strbuf, "%s/%s_e%.2et%.2f%s.%d.%c.mem", dir.c_str(), file_prefix, eta, tau_hom, label.c_str(), world.rank (), 'e');
642
+ // sprintf(strbuf, "%s/%s_e%.2et%.2f%s.%d.%c.mem", dir.c_str(), file_prefix, eta, tau_hom, label.c_str(), sys->mpi_rank (), 'e');
652
643
// VoltageMonitor * vmon_e = new VoltageMonitor( neurons_e, 33, strbuf );
653
644
654
- // sprintf(strbuf, "%s/%s_e%.2et%.2f%s.%d.%c.ras", dir.c_str(), file_prefix, eta, tau_hom, label.c_str(), world.rank (), 'i');
645
+ // sprintf(strbuf, "%s/%s_e%.2et%.2f%s.%d.%c.ras", dir.c_str(), file_prefix, eta, tau_hom, label.c_str(), sys->mpi_rank (), 'i');
655
646
// SpikeMonitor * smon_i = new SpikeMonitor( neurons_i, strbuf , 500);
656
647
657
648
RateChecker * chk = new RateChecker ( neurons_e , 0.1 , 20 .*kappa , tau_chk);
@@ -669,19 +660,19 @@ int main(int ac, char* av[])
669
660
oss << " Changing cell input ... " ;
670
661
logger->msg (oss.str (),PROGRESS,true );
671
662
672
- sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.pat" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), world. rank ());
663
+ sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.pat" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), sys-> mpi_rank ());
673
664
PatternMonitor * patmon = new PatternMonitor (neurons_e,strbuf,patfile.c_str (), 10 , 1 .);
674
665
PatternStimulator * patstim = new PatternStimulator (neurons_e,currentfile.c_str (),patfile.c_str ());
675
666
676
- sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.scal" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), world. rank ());
667
+ sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.scal" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), sys-> mpi_rank ());
677
668
WeightMonitor * wmon_scal = new WeightMonitor ( con_ee, 0 , 0 , strbuf, 1 , ELEMENTLIST);
678
669
wmon_scal->add_to_list ( con_ee->get_pre_partners (0 ) );
679
670
wmon_scal->add_to_list ( con_ee->get_pre_partners (10 ) );
680
671
681
672
}
682
673
683
674
// if ( corr ) {
684
- // sprintf(strbuf, "%s/%s_e%.2et%.2f%s.%d.%c.pat", dir.c_str(), file_prefix, eta, tau_hom, label.c_str(), world.rank (), 'c');
675
+ // sprintf(strbuf, "%s/%s_e%.2et%.2f%s.%d.%c.pat", dir.c_str(), file_prefix, eta, tau_hom, label.c_str(), sys->mpi_rank (), 'c');
685
676
// PatternMonitor * pmon_c = new PatternMonitor( neurons_e, strbuf,corr_pat_file.c_str() );
686
677
// }
687
678
@@ -719,7 +710,7 @@ int main(int ac, char* av[])
719
710
SparseConnection * con_corr_i = new SparseConnection (corr_e,neurons_i,w*strength,sparseness*plen*psize/ni,GLUT);
720
711
721
712
// set up Weight monitor
722
- sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.%c.ras" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), world. rank (), ' c' );
713
+ sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.%c.ras" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), sys-> mpi_rank (), ' c' );
723
714
SpikeMonitor * smon_c = new SpikeMonitor ( corr_e, strbuf , psize );
724
715
}
725
716
@@ -762,7 +753,7 @@ int main(int ac, char* av[])
762
753
// SparseConnection * con_corr_i = new SparseConnection(corr_e,neurons_i,w*strength,sparseness*plen*psize/ni,GLUT);
763
754
764
755
// set up Weight monitor
765
- sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.%c.ras" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), world. rank (), ' c' );
756
+ sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.%c.ras" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), sys-> mpi_rank (), ' c' );
766
757
SpikeMonitor * smon_c = new SpikeMonitor ( corr_e, strbuf , psize );
767
758
}
768
759
@@ -796,7 +787,7 @@ if ( prefile != "" ) {
796
787
// adding patterns
797
788
if ( patfile != " " ) {
798
789
logger->msg (" Preparing stimulus ..." ,PROGRESS,true );
799
- sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.stim" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), world. rank ());
790
+ sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.stim" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), sys-> mpi_rank ());
800
791
StimulusGroup * stimgroup = new StimulusGroup (ne,patfile,strbuf);
801
792
stimgroup->set_mean_on_period (onperiod);
802
793
stimgroup->set_mean_off_period (offperiod);
@@ -829,33 +820,33 @@ if ( patfile != "" ) {
829
820
830
821
if (!fast) {
831
822
logger->msg (" Saving neurons state ..." ,PROGRESS,true );
832
- sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.e.nstate" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), world. rank ());
823
+ sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.e.nstate" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), sys-> mpi_rank ());
833
824
neurons_e->write_to_file (strbuf);
834
- sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.i.nstate" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), world. rank ());
825
+ sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.i.nstate" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), sys-> mpi_rank ());
835
826
neurons_i->write_to_file (strbuf);
836
827
837
828
logger->msg (" Saving weight matrix ..." ,PROGRESS,true );
838
- sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.wmat" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), world. rank ());
829
+ sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.wmat" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), sys-> mpi_rank ());
839
830
con_ee->write_to_file (strbuf);
840
831
841
- sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.ei.wmat" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), world. rank ());
832
+ sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.ei.wmat" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), sys-> mpi_rank ());
842
833
con_ei->write_to_file (strbuf);
843
834
844
835
for ( int i = 0 ; i < corr_connections.size () ; ++i ) {
845
- sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.%d.wmat" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), i, world. rank ());
836
+ sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.%d.wmat" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), i, sys-> mpi_rank ());
846
837
corr_connections[i]->write_to_file (strbuf);
847
838
}
848
839
}
849
840
850
841
// save lifetime
851
- sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.lifetime" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), world. rank ());
842
+ sprintf (strbuf, " %s/%s_e%.2et%.2f%s.%d.lifetime" , dir.c_str (), file_prefix, eta, tau_hom, label.c_str (), sys-> mpi_rank ());
852
843
std::ofstream killfile;
853
844
killfile.open (strbuf);
854
845
killfile << sys->get_time ()-primetime << std::endl;
855
846
killfile.close ();
856
847
857
848
858
- if ( world. rank () == 0 ) {
849
+ if ( sys-> mpi_rank () == 0 ) {
859
850
logger->msg (" Saving elapsed time ..." ,PROGRESS,true );
860
851
char filenamebuf [255 ];
861
852
sprintf (filenamebuf, " %s/elapsed.dat" , dir.c_str ());
@@ -869,7 +860,7 @@ if ( patfile != "" ) {
869
860
delete sys;
870
861
871
862
if (errcode) {
872
- env. abort (errcode);
863
+ mpienv-> abort (errcode);
873
864
}
874
865
875
866
return errcode;
0 commit comments