Skip to content

Commit

Permalink
Build sort of gmv parallel wrapper
Browse files Browse the repository at this point in the history
  • Loading branch information
@eaulisa
eaulisa committed Mar 23, 2015
1 parent 7388fe1 commit 53adbea
Show file tree
Hide file tree
Showing 2 changed files with 75 additions and 63 deletions.
134 changes: 73 additions & 61 deletions src/solution/GMVWriter.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -24,7 +24,8 @@
#include <iostream>
#include <sstream>
#include <algorithm>
#include <cstring>
#include <cstring>
#include "Files.hpp"


namespace femus {
Expand Down Expand Up @@ -331,6 +332,11 @@ void GMVWriter::write(const std::string output_path, const char order[], const s

void GMVWriter::Pwrite(const std::string output_path, const char order[], const std::vector<std::string>& vars, const unsigned time_step) const {

if(_nprocs == 1){
write( output_path, order, vars, time_step);
return;
}

unsigned gridn = _gridn; // aggiunta da me

if ( gridn == 0 ) gridn = _gridn;
Expand All @@ -340,20 +346,21 @@ void GMVWriter::Pwrite(const std::string output_path, const char order[], const
// ********** linear -> index==0 *** quadratic -> index==1 **********
unsigned index=( strcmp(order, "linear") ) ? 1 : 0;

std::string dirnamePGMV = "GMVParallelFiles/";
Files files;
files.CheckDir(output_path,dirnamePGMV);

std::string filename_prefix;
if( _ml_sol != NULL ) filename_prefix = "sol";
else filename_prefix = "mesh";

std::ostringstream filename;
filename << output_path << "/" << filename_prefix << ".level" << _gridn << "." <<_iproc<<"."<< time_step << "." << order << ".gmv";
filename << output_path << "/" << dirnamePGMV << filename_prefix << ".level" << _gridn << "." <<_iproc<<"."<< time_step << "." << order << ".gmv";

std::ofstream fout;

fout.open(filename.str().c_str());
if (fout.is_open()) {
std::cout << std::endl << " The output is printed to file " << filename.str() << " in GMV format" << std::endl;
}
else {
if (!fout.is_open()) {
std::cout << std::endl << " The output file "<< filename.str() <<" cannot be opened.\n";
abort();
}
Expand Down Expand Up @@ -700,61 +707,66 @@ void GMVWriter::Pwrite(const std::string output_path, const char order[], const
}
delete [] det;

// // print xml wrapper
//
// // *********** open pvtu file *************
// std::ofstream Pfout;
// if(_iproc!=0) {
// Pfout.rdbuf(); //redirect to dev_null
// }
// else {
// std::ostringstream Pfilename;
// Pfilename << output_path << "/" << filename_prefix << "_level" << _gridn <<"_"<< time_step << "_" << order << "_gmv.xml";
// Pfout.open(Pfilename.str().c_str());
// if (Pfout.is_open()) {
// std::cout << std::endl << " The output is printed to file " << Pfilename.str() << " in parallel GMV-XML (64-based) format" << std::endl;
// }
// else {
// std::cout << std::endl << " The output file "<< Pfilename.str() <<" cannot be opened.\n";
// abort();
// }
// }
//
// // *********** write pvtu header ***********
// Pfout<< "<?xml version=\"1.0\"?>" << std::endl;
// Pfout<< "<pqevents>\n";
// for(int jproc=0;jproc<_nprocs;jproc++){
// Pfout<< " <pqevent object=\"pqClientMainWindow/MainControlsToolbar/1QToolButton0\" command=\"activate\" arguments=\"\" />\n";
// Pfout<< " <pqevent object=\"pqClientMainWindow/FileOpenDialog\" command=\"filesSelected\" arguments=\"";
// Pfout<< filename_prefix << ".level" << _gridn << "." <<jproc<<"."<< time_step << "." << order << ".gmv\" />\n";
// }
//
// Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mousePress\" arguments=\"1,1,0,50,0,/0:0\" />\n";
// Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mouseRelease\" arguments=\"1,0,0,50,0,/0:0\" />\n";
//
// Pfout<< " <pqevent object=\"pqClientMainWindow/propertiesDock/propertiesPanel/Accept\" command=\"activate\" arguments=\"\" />\n";
//
// Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mousePress\" arguments=\"1,1,33554432,50,0,/0:0/"<<_nprocs-1<<":0\" />\n";
// Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mouseRelease\" arguments=\"1,0,33554432,50,0,/0:0/"<<_nprocs-1<<":0\" />\n";
// Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"keyEvent\" arguments=\"7,16777248,0,,0,1\" />\n";
// Pfout<< " <pqevent object=\"pqClientMainWindow/menubar\" command=\"activate\" arguments=\"menuFilters\" />\n";
// Pfout<< " <pqevent object=\"pqClientMainWindow/menubar/menuFilters/Alphabetical\" command=\"activate\" arguments=\"GroupDataSets\" />\n";
//
// Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mousePress\" arguments=\"1,1,0,50,0,/0:0\" />\n";
// Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mouseRelease\" arguments=\"1,0,0,50,0,/0:0\" />\n";
// Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mousePress\" arguments=\"1,1,33554432,50,0,/0:0/"<<_nprocs<<":0\" />\n";
// Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mouseRelease\" arguments=\"1,0,33554432,50,0,/0:0/"<<_nprocs<<":0\" />\n";
// Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"keyEvent\" arguments=\"7,16777248,0,,0,1\" />\n";
// Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mousePress\" arguments=\"1,1,67108864,0,0,/0:0/0:1\" />\n";
// Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mouseRelease\" arguments=\"1,0,67108864,0,0,/0:0/0:1\" />\n";
//
// Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mousePress\" arguments=\"1,1,0,50,0,/0:0/"<<_nprocs<<":0\" />\n";
// Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mouseRelease\" arguments=\"1,0,0,50,0,/0:0/"<<_nprocs<<":0\" />\n";
// Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mousePress\" arguments=\"1,1,0,0,0,/0:0/"<<_nprocs<<":1\" />\n";
// Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mouseRelease\" arguments=\"1,0,0,0,0,/0:0/"<<_nprocs<<":1\" />\n";
//
// Pfout<< "</pqevents>\n";
// Pfout.close();
// print xml wrapper

// *********** open pvtu file *************
std::ofstream Pfout;
if(_iproc!=0) {
Pfout.rdbuf(); //redirect to dev_null
}
else {
std::ostringstream Pfilename;
Pfilename << output_path << "/" << filename_prefix << "_level" << _gridn <<"_"<< time_step << "_" << order << "_gmv.xml";
Pfout.open(Pfilename.str().c_str());
if (Pfout.is_open()) {
std::cout << std::endl << " The output is printed to file " << Pfilename.str() << " in parallel GMV-XML (64-based) format" << std::endl;
std::cout << " Laod it in ParaView -> Tools -> Play Test ... \n";
}
else {
std::cout << std::endl << " The output file "<< Pfilename.str() <<" cannot be opened.\n";
abort();
}
}

// *********** write pvtu header ***********
Pfout<< "<?xml version=\"1.0\"?>" << std::endl;
Pfout<< "<pqevents>\n";
for(int jproc=0;jproc<_nprocs;jproc++){
Pfout<< " <pqevent object=\"pqClientMainWindow/MainControlsToolbar/1QToolButton0\" command=\"activate\" arguments=\"\" />\n";
Pfout<< " <pqevent object=\"pqClientMainWindow/FileOpenDialog\" command=\"filesSelected\" arguments=\"";
Pfout<< dirnamePGMV << filename_prefix << ".level" << _gridn << "." <<jproc<<"."<< time_step << "." << order << ".gmv\" />\n";
}

Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mousePress\" arguments=\"1,1,0,50,5,/0:0\" />\n";
Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mouseRelease\" arguments=\"1,0,0,50,5,/0:0\" />\n";

Pfout<< " <pqevent object=\"pqClientMainWindow/propertiesDock/propertiesPanel/Accept\" command=\"activate\" arguments=\"\" />\n";

Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mousePress\" arguments=\"1,1,33554432,50,5,/0:0/"<<_nprocs-1<<":0\" />\n";
Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mouseRelease\" arguments=\"1,0,33554432,50,5,/0:0/"<<_nprocs-1<<":0\" />\n";
Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"keyEvent\" arguments=\"7,16777248,0,,0,1\" />\n";
Pfout<< " <pqevent object=\"pqClientMainWindow/menubar\" command=\"activate\" arguments=\"menuFilters\" />\n";
Pfout<< " <pqevent object=\"pqClientMainWindow/menubar/menuFilters/Alphabetical\" command=\"activate\" arguments=\"GroupDataSets\" />\n";

Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mousePress\" arguments=\"1,1,0,10,5,/0:0/"<<_nprocs<<":1\" />\n";
Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mouseRelease\" arguments=\"1,0,0,10,5,/0:0/"<<_nprocs<<":1\" />\n";

Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mousePress\" arguments=\"1,1,33554432,25,5,/0:0\" />\n";
Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mouseRelease\" arguments=\"1,0,33554432,25,5,/0:0\" />\n";
Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mousePress\" arguments=\"1,1,33554432,50,5,/0:0/"<<_nprocs<<":0\" />\n";
Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mouseRelease\" arguments=\"1,0,33554432,50,5,/0:0/"<<_nprocs<<":0\" />\n";
Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"keyEvent\" arguments=\"7,16777248,0,,0,1\" />\n";
Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mousePress\" arguments=\"1,1,0,10,5,/0:0/"<<_nprocs<<":1\" />\n";
Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mouseRelease\" arguments=\"1,0,0,10,5,/0:0/"<<_nprocs<<":1\" />\n";

Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mousePress\" arguments=\"1,1,0,50,5,/0:0/"<<_nprocs<<":1\" />\n";
Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mouseRelease\" arguments=\"1,0,0,50,5,/0:0/"<<_nprocs<<":1\" />\n";

Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mousePress\" arguments=\"1,1,0,10,5,/0:0/"<<_nprocs<<":1\" />\n";
Pfout<< " <pqevent object=\"pqClientMainWindow/pipelineBrowserDock/pipelineBrowser\" command=\"mouseRelease\" arguments=\"1,0,0,10,5,/0:0/"<<_nprocs<<":1\" />\n";

Pfout<< "</pqevents>\n";
Pfout.close();



Expand Down
4 changes: 2 additions & 2 deletions src/utils/Files.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -216,13 +216,13 @@ void Files::CheckDir(const std::string& dir_name_in, const std::string& my_name_
std::ifstream in; in.open(dirname.str().c_str());

if (!in.is_open()) {
std::cout << " No " << dirname.str() << " directory: I'll create it. " << std::endl;
std::cout << std::endl << " No " << dirname.str() << " directory: I'll create it. " << std::endl;

int status = mkdir(dirname.str().c_str(), S_IRWXU | S_IRGRP | S_IXGRP | S_IROTH | S_IXOTH);
if (status != 0) {std::cout << "MKDIR error: " << status << std::endl; abort();}

}
else {std::cout<< "That's alright, " << my_name_in << " is already there." << std::endl;}
else {std::cout << std::endl <<" That's alright, " << my_name_in << " is already there." << std::endl;}

in.close();

Expand Down

0 comments on commit 53adbea

Please sign in to comment.