scalapack+C+mpi beginner

Open discussion regarding features, bugs, issues, vendors, etc.

scalapack+C+mpi beginner

Postby Raimondas » Thu Dec 02, 2010 5:44 am

Hello,

1. I have some knowledge in C/C++ programming, because i am electronic engineering student.
2. I have created Mpich 2 cluster and installed scalapack.
3. And most important I need to calculate Ax=b, where A is N by N matrix and b is N row.
4. From google i found some examples, but they are without comments and difficult to understood.
5. Please help me !!

my simplest code is bellow:
Code: Select all
#include <stdio.h>
#include <stdlib.h>
#include "mpi.h"

int main(int argc, char **argv)
{

   int rank, size;

   MPI_Init( &argc, &argv);
      MPI_Comm_rank(MPI_COMM_WORLD, &rank);
      MPI_Comm_size(MPI_COMM_WORLD, &size);
   
   int N;
   N=9;

   //2D matrix converted to 1D by row
   double *A=(double*)malloc(N*N*sizeof(double));

   A[0]=1;
   A[1]=-0.5;
   A[2]=0;
   A[3]=-0.5;
   A[4]=0;
   A[5]=0;
   A[6]=0;
   A[7]=0;
   A[8]=0;
   A[9]=-0.5;
   A[10]=2;
   A[11]=-0.5;
   A[12]=0;
   A[13]=-1;
   A[14]=0;
   A[15]=0;
   A[16]=0;
   A[17]=0;
   A[18]=0;
   A[19]=-0.5;
   A[20]=1;
   A[21]=0;
   A[22]=0;
   A[23]=-0.5;
   A[24]=0;
   A[25]=0;
   A[26]=0;
   A[27]=-0.5;
   A[28]=0;
   A[29]=0;
   A[30]=2;
   A[31]=-1;
   A[32]=0;
   A[33]=-0.5;
   A[34]=0;
   A[35]=0;
   A[36]=0;
   A[37]=-1;
   A[38]=0;
   A[39]=-1;
   A[40]=4;
   A[41]=-1;
   A[42]=0;
   A[43]=-1;
   A[44]=0;
   A[45]=0;
   A[46]=0;
   A[47]=-0.5;
   A[48]=0;
   A[49]=-1;
   A[50]=2;
   A[51]=0;
   A[52]=0;
   A[53]=-0.5;
   A[54]=0;
   A[55]=0;
   A[56]=0;
   A[57]=-0.5;
   A[58]=0;
   A[59]=0;
   A[60]=1;
   A[61]=-0.5;
   A[62]=0;
   A[63]=0;
   A[64]=0;
   A[65]=0;
   A[66]=0;
   A[67]=-1;
   A[68]=0;
   A[69]=-0.5;
   A[70]=2;
   A[71]=-0.5;
   A[72]=0;
   A[73]=0;
   A[74]=0;
   A[75]=0;
   A[76]=0;
   A[77]=-0.5;
   A[78]=0;
   A[79]=-0.5;
   A[80]=1;
   
   //RHS
   double *B=(double*)malloc(N*sizeof(double));
   B[0]=0;
   B[1]=0;
   B[2]=0;
   B[3]=0;
   B[4]=1;
   B[5]=0;
   B[6]=0;
   B[7]=0;
   B[8]=0;

   //=============RESULT MUST BE=============//
   //0.8      0.9      0.8
   //0.7      1      0.7
   //0         0      0
   //========================================//

   MPI_Finalize();
      exit(0);
}
Raimondas
 
Posts: 1
Joined: Wed Dec 01, 2010 12:04 pm

Re: scalapack+C+mpi beginner

Postby admin » Tue Feb 08, 2011 9:20 am

Hi,
we have a User's guide available online.
See http://www.netlib.org/scalapack/slug/index.html
Julie
admin
Site Admin
 
Posts: 504
Joined: Wed Dec 08, 2004 7:07 pm


Return to User Discussion

Who is online

Users browsing this forum: No registered users and 2 guests