Turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Intel Community
- Software
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library
- Use OpenMP when invoking MKL poisson library

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

Chang_Lei

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

08-02-2010
01:43 AM

85 Views

Use OpenMP when invoking MKL poisson library

It seems that the accelerating speed is tiny.

Is this ordinary to the MKL poisson library(or where could I find the PARALLEL PERFORMANCE of MKL Poisson solver) ?

Thanks.

1 Solution

Alexander_K_Intel2

Employee

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

08-02-2010
07:38 PM

85 Views

Try to rewrite your code in next way:

call d_init_Helmholtz_3D(ax,bx,...,...,BCtype,q,ipar,dpar,stat)

call d_commit_Helmholtz_3D(ff,bd_ax,bd_bx,...,...,dpar,stat)

call d_Helmholtz_3D(ff,bd_ax,,...,...,ipar,dpar,stat)

call free_Helmholtz_3D(xhandle,yhandle,ipar,stat)

With best regards,

Alexander Kalinkin

Link Copied

5 Replies

TimP

Black Belt

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

08-02-2010
05:37 AM

85 Views

The goal of OpenMP parallelism is not to reduce CPU time, but to reduce elapsed time. You are lucky if the total CPU time is reduced; more often it would increase somewhat.

Chang_Lei

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

08-02-2010
06:45 PM

85 Views

The exact command is as follows:*ifort -openmp -fast mkl_dfti.f90 mkl_poisson.f90 MyProgram.f90 -L/export/mathlib/cmkl90/lib/em64t/ -lmkl_em64t -lguide -lpthread -lm*And the code where I invoked Poisson library isin MyProgram.f90:

call d_init_Helmholtz_3D(ax,bx,...,...,BCtype,q,ipar,dpar,stat)

call d_commit_Helmholtz_3D(ff,bd_ax,bd_bx,...,...,dpar,stat)

ipar(23)=4 !OpenMP number of threads

call d_Helmholtz_3D(ff,bd_ax,,...,...,ipar,dpar,stat)

call free_Helmholtz_3D(xhandle,yhandle,ipar,stat)

call cpu_time(t1)

The program is submitted and run on cluster, and 4 threads are used(setin *.sh file as"nodes=1:ppn=4").

So, if I'm wrong. How can I get the right usage? You didn't mean the CMKL9.0 doesn't support OpenMP, did you?

Anywhere can I find the usage of Poisson library with OpenMP?

With best regards.

Alexander_K_Intel2

Employee

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

08-02-2010
07:38 PM

86 Views

Try to rewrite your code in next way:

call d_init_Helmholtz_3D(ax,bx,...,...,BCtype,q,ipar,dpar,stat)

call d_commit_Helmholtz_3D(ff,bd_ax,bd_bx,...,...,dpar,stat)

call d_Helmholtz_3D(ff,bd_ax,,...,...,ipar,dpar,stat)

call free_Helmholtz_3D(xhandle,yhandle,ipar,stat)

With best regards,

Alexander Kalinkin

Chang_Lei

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

08-03-2010
02:19 AM

85 Views

I have rewrote as you said anddecreased the elapsed time from about 0.17 to 0.06.

That seems to be very good accelerating speed.

Here I have another question. Can MKL poisson library be used in cylindrical coordinate (r,z,thelta)? Especially for the axisymmetrical coordinate (r,z). (I'm holding my breath...)

Best regards.

Alexander_K_Intel2

Employee

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

08-03-2010
04:02 AM

85 Views

Current version of MKL support only cartesian and spherical coordinate system.

With best regards,

Alexander Kalinkin

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

For more complete information about compiler optimizations, see our Optimization Notice.