[Trilinos-Users] scaling with trilinos tpetra and anasazi solvers
Mike Atambo
mikeat4999 at gmail.com
Thu Jan 29 09:15:24 MST 2015
Mike (Heroux),
Thanks i will do this, and see what kind of information i get, i will
post back what i find to
make it easier to for one of you to help (or explain what i may be
missing).
Mike
On Thu, Jan 29, 2015 at 4:29 PM, Heroux, Mike <MHeroux at csbsju.edu> wrote:
> Mike,
>
> Just to follow up a bit more: The reason I suggest this manual approach
> is to determine if your MPI environment can scale properly. If your cores
> are on a single workstation, you will not get full speed up on this kind of
> solver, but I am surprised you are not seeing any. This test would help
> diagnose the issue.
>
> Mike
>
> From: <Heroux>, "maherou at sandia.gov" <maherou at sandia.gov>
> Date: Thursday, January 29, 2015 at 1:03 PM
> To: Mike Atambo <mikeat4999 at gmail.com>, "
> Trilinos-Users at software.sandia.gov" <Trilinos-Users at software.sandia.gov>
> Subject: Re: [Trilinos-Users] scaling with trilinos tpetra and anasazi
> solvers
>
> Mike,
>
>
> It is hard to say for sure why you are seeing the results you have. One
> thing you could do is write a simple MPI code that does a vector operation
> such as:
>
>
> int localn = 50000000/size; // size is the number of MPI ranks
>
>
> // define arrays x, y, w appropriately
>
> // define ntrials so that the job runs a significant amount of time.
>
>
> for (int j=0; j< ntrials; ++j)
>
> for (int i=0; i< localn; ++i) w[i] += alpha*x[i] + y[i];
>
>
> With a tridiagonal system and the power method, your computation and
> data access should be similar to the above operation. This should give you
> some sense of what to expect.
>
>
> Mike
>
>
> ------------------------------
> *From:* trilinos-users-bounces at software.sandia.gov <
> trilinos-users-bounces at software.sandia.gov> on behalf of Mike Atambo <
> mikeat4999 at gmail.com>
> *Sent:* Monday, January 26, 2015 7:29 AM
> *To:* Trilinos-Users at software.sandia.gov
> *Subject:* [EXTERNAL] [Trilinos-Users] scaling with trilinos tpetra and
> anasazi solvers
>
> My apology, i send an earlier email to an older thread, this is a
> correction.
>
> Im trying to solve a system using tpetra and anasazi's krylov-shur,
> but there seems to
> be no scaling with number of processors, i tried to run some of the
> examples, and
> they reflect the same things i saw in my code,
>
> attached (If list accepts attachments... hoping) is the lesson_03
> power method, changed to create a tridiagonal sparse matrix,
> with about 50 000 000 global indices, the run time seems the same
> for 2 to 16 processors,
> is there something im missing?
>
>
>
> $ time mpirun -np 2 ./powerm.x
> real 3m16.014s
> user 3m5.610s
> sys 0m4.452s
>
> $ time mpirun -np 8 ./powerm.x | grep real
> real 3m14.229s
> user 3m7.938s
> sys 0m4.939s
>
>
> $ time mpirun -np 16 ./powerm.x
> real 3m18.153s
> user 3m10.202s
> sys 0m6.108s
>
> Mike.
>
> matambo at ictp.it
> Ext .139
>
>
>
--
M. O. Atambo
mikeat4999 at gmail.com
matambo at ictp.it
Ext .139
Room 209.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://software.sandia.gov/pipermail/trilinos-users/attachments/20150129/72579980/attachment.html>
More information about the Trilinos-Users
mailing list