Turn on suggestions

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

Showing results for

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

pe_rasouliyahoo_com

Beginner

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

08-23-2012
12:29 PM

100 Views

1D Poisson example

d_Poisson_2D_f.f90 error:

=======================================================================

The number of mesh intervals in x-direction is nx=6

The number of mesh intervals in y-direction is ny=1

In the mesh point (0.167,0.000) the error between the computed and the true solution is equal to -0.272E+00

In the mesh point (0.167,1.000) the error between the computed and the true solution is equal to 0.272E+00

The computed solution seems to be inaccurate.

Double precision 2D Poisson example FAILED to compute the solution...

and d_Poisson_3D_f.f90 error:

------------------------------------------------------------------------------

MKL TRIG TRANSFORMS ERROR:

The dimension of the trigonometric transform 1 should be

an integer number greater or equal to 2.

-------------------------------------------------------------------------------

MKL POISSON LIBRARY ERROR:

Fatal error: Trigonometric Transform commit step has failed to complete.

Error code = 1. Computations has been stopped...

Double precision 3D Poisson example FAILED to compute the solution...

Is there any quick resolution for it or is there a 1D Poisson fortran example out there?

thanks

Link Copied

1 Reply

Alexander_K_Intel2

Employee

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

08-23-2012
07:26 PM

100 Views

If you want to get solution of 1D Poisson solver of size N and rhs with element f(k) you need to made next steps:

1. Set Nz = N, Nx and Ny any number greater or equal to 2.

2. Set zero Neumann condition in both x and y dimension

3. Set rhs(i,j,k) = f(k)

In such case you get Nx*Ny solution of 1D Poisson equation (x(k) = solution(i,j,k) independently of i and j).

But I'm strongly recommend you to use LAPACK functionality to solve such equation

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.