#ifndef lint
static char vcid[] = "$Id: $";
#endif

#include "sys.h"
#include "options.h"
#include "sysio.h"
#include "mat.h"

/*@C
      MatCreateInitialMatrix - Reads from command line to determine 
           what type of matrix to create. Also uses MPI matrices if 
           number processors in MPI_COMM_WORLD greater then one.

  Input Parameters:
.   m,n - matrix dimensions
 
  Output Parameter:
.   V - location to stash resulting matrix.
@*/
int MatCreateInitialMatrix(int m,int n,Mat *V)
{
  int numtid;
  MPI_Comm_size(MPI_COMM_WORLD,&numtid);
  if (OptionsHasName(0,0,"-dense_mat")) {
    return MatCreateSequentialDense(m,n,V);
  }
  if (OptionsHasName(0,0,"-mpi_mats") || numtid >1) {
    return MatCreateMPIAIJ(MPI_COMM_WORLD,-1,-1,m,n,5,0,0,0,V);
  }
  return MatCreateSequentialAIJ(m,n,10,0,V);
}