Function main

SRC/tcl/mpiMainTest.cpp:52–134

#include <Matrix.h>

int
main(int argc, char **argv)
{
  FEM_ObjectBroker theBroker;
  MPI_MachineBroker theMachine(&theBroker, argc, argv);
  theMachineBroker = &theMachine;

  int rank = theMachine.getPID();
  int np = theMachine.getNP();

  //
  // depending on rank we do something
  //
  if (rank != 0) {

    //
    // on secondary processes we spin waiting to create & run actors
    //
    fprintf(stderr, "Secondary Process Running\n");
    theMachine.runActors();

  } else {

    //
    // on process 0 we create some ShadowSubdomains & then start the OpenSees interpreter
    //
    fprintf(stderr, "Primary Process Running OpenSees Interpreter\n");   

    //
    // set some global parameters
    //
    OPS_OBJECT_BROKER = &theBroker;
    OPS_MACHINE = &theMachine;
    OPS_PARALLEL_PROCESSING = np;
    OPS_NUM_SUBDOMAINS = np - 1;
    OPS_PARTITIONED    = false;
    Channel  **OPS_theChannels = 0;
    OPS_theChannels = new Channel *[OPS_NUM_SUBDOMAINS];
    Subdomain **OPS_theSubdomains = 0;
    OPS_theSubdomains = new Subdomain *[OPS_NUM_SUBDOMAINS];

    // start the remote actors
    for (int i=1; i<=OPS_NUM_SUBDOMAINS; i++) {
      ShadowSubdomain *theSubdomain = new ShadowSubdomain(i, 
							  *OPS_MACHINE, 
							  *OPS_OBJECT_BROKER);
      theDomain.addSubdomain(theSubdomain);
      OPS_theChannels[i-1] = theSubdomain->getChannelPtr();
      OPS_theSubdomains[i-1] = theSubdomain;
    }

    Subdomain *theSubdomain = OPS_theSubdomains[0];
    Node *node1 = new Node(1, 2,   0.0,  0.0);
    Node *node2 = new Node(2, 2, 144.0,  0.0);
    Node *node3 = new Node(3, 2, 168.0,  0.0);    
    Node *node4 = new Node(4, 2,  72.0, 96.0);        
    Matrix mass(2,2); mass(0,0)=1.0; mass(1,1)=2.0;
    node1->setMass(mass);