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
- Problems with 2D poisson solver and periodic conditions

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

Camilo_P_

Beginner

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

04-12-2016
12:30 PM

112 Views

Problems with 2D poisson solver and periodic conditions

Hello,

I am coding a fluid solver in the Vorticity-Potential formulation. I need to solve a poisson equation for the stream function. I was trying to use MKL's. In my problem, I have Dirichlet BCs in the y walls and periodic in the x walls. When I solve the poisson equation with the MKL solver I get wrong results at the corners, but when I solve it using SOR I get the correct results.

I know that I am using the solver appropiately because I have tested it with other functions and I have gotten correct results.

I can see from the results that the periodicity is not being enforced when MKL is used.

I am using Intel Visual Fortran XE 2011 Update 10. I have attached the code, I use the variable *useMKL* to switch between MKL solver and SOR. It would be great if I can use MKL's as I need to move to a 3D implementation and the SOR method will be very slow then.

If more information is needed let me know. Thanks in advance,

Camilo

Link Copied

3 Replies

Gennady_F_Intel

Moderator

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

04-12-2016
09:40 PM

112 Views

Camilo_P_

Beginner

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

04-13-2016
07:23 AM

112 Views

I found my error, in line 205 it read:

do ix=2,nx

it should have read

do ix=1,nx+1

It works perfect now.

Chatterjee__Swagnik

Beginner

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

07-12-2018
11:39 PM

112 Views

Hello Sir,

Can you please explain me how to solve a diffusion equation using the mkl poisson solver with periodic in X and Y and neumann in Z ? i am very new in this field and the exact algorithm is not clear to me.

Thanks,

Swagnik

- 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.