plot3d.C File Reference

#include "roccom.h"
#include <cstdio>
#include <iostream>
#include <fstream>
#include <vector>
#include <algorithm>
#include <cstring>
#include <string>
#include <cstdlib>
#include <cmath>
#include <cassert>
#include "commpi.h"

Include dependency graph for plot3d.C:

Go to the source code of this file.

Classes
class	Point_3< T >

Functions
	COM_EXTERN_MODULE (Rocface)
void	read_plot3d (istream &is, vector< vector< Point_3 > > &coors, vector< pair< int, int > > &dims)
int	main (int argc, char *argv[])

Function Documentation

COM_EXTERN_MODULE

(

Rocface

)

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 77 of file plot3d.C.

References COM_call_function(), COM_finalize(), COM_get_attribute_handle(), COM_get_function_handle(), COM_init(), COM_LOAD_MODULE_STATIC_DYNAMIC, COM_new_window(), COM_print_profile(), COM_set_array(), COM_set_profiling(), COM_set_verbose(), COM_window_init_done(), i, k, MPI_COMM_WORLD, ni, and read_plot3d().

                                 {
  MPI_Init( &argc, &argv);
  COM_init( &argc, &argv);

  if ( argc < 3) {
    std::cout << "Usage: " << argv[0] 
              << " file1 file2 [<out_path>] " 
              << " [<scale_1>] [<scale_2>]" << std::endl;
    exit(-1);
  }

  COM_set_verbose( 1);
  COM_set_profiling( 1);

  const string wnames[2] = {string(argv[1], string( argv[1]).find(".")), 
                            string(argv[2], string( argv[2]).find("."))};

  vector< vector<Point_3> >    coors[2];

  COM_LOAD_MODULE_STATIC_DYNAMIC( Rocface, "RFC");

  for ( int k=0; k<2; ++k) {
    COM_new_window( wnames[k].c_str());

    ifstream fi( argv[k+1]); 
    vector<pair< int, int> >  dims;
    read_plot3d( fi, coors[k], dims);

    for ( int i=0, ni=dims.size(); i<ni; ++i) {
      COM_set_array( (wnames[k]+".:st2:").c_str(), i+1, &dims[i].first);
      COM_set_array( (wnames[k]+".nc").c_str(), i+1, &coors[k][i][0]);
    }

    COM_window_init_done( wnames[k].c_str());
  }

  MPI_Comm comm = MPI_COMM_WORLD;

  int RFC_overlay = COM_get_function_handle( "RFC.overlay");
  int RFC_write = COM_get_function_handle( "RFC.write_overlay");

  int tri1_mesh = COM_get_attribute_handle( (wnames[0]+".mesh").c_str());
  int tri2_mesh = COM_get_attribute_handle( (wnames[1]+".mesh").c_str());

  COM_call_function( RFC_overlay, &tri1_mesh, &tri2_mesh, &comm, "./");
  COM_call_function( RFC_write, &tri1_mesh, &tri2_mesh);

  COM_print_profile( "", "");

  COM_finalize();
  MPI_Finalize();  
}

Here is the call graph for this function:

void read_plot3d	(	istream &	is,
		vector< vector< Point_3 > > &	coors,
		vector< pair< int, int > > &	dims
	)

Definition at line 55 of file plot3d.C.

References i, j, and nj.

Referenced by main().

                                                 {
  int ngrid;
  is >> ngrid;

  dims.resize(ngrid);
  for ( int i=0; i<ngrid; ++i) {
    int t;
    is >> dims[i].first >> dims[i].second >> t;
  }

  coors.resize( ngrid);
  for ( int i=0; i<ngrid; ++i) {
    coors[i].resize( dims[i].first*dims[i].second);

    for ( int l=0; l<3; ++l) {
      for ( int j=0, nj=coors[i].size(); j<nj; ++j)
        is >> coors[i][j][l];
    }
  }
}

Here is the caller graph for this function: