CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

Linear Solver for Particle Methods

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   November 30, 2013, 15:50
Default Linear Solver for Particle Methods
Senior Member
Prapanch Nair
Join Date: Mar 2009
Location: Bangalore, India
Posts: 105
Rep Power: 8
prapanj is on a distinguished road
Dear All,

I am using a particle method to solve a fluid mechanics problem. This involves solving a poisson equation every time step. I have my BiCGSTAB code to solve the linear system. But since the particle configuration changes with time, I have to build a new coefficient matrix for every time step. So, I believe there must be a linear solver library out there that can help me do this efficiently.

I started with LIS (since I wish to use a library that is simple with less dependencies). But the problem is that, since the coeff matrix changes size and sparsity every time step, the sparse matrix build involves lot of reallocation of memory and is costing a lot of computation time. I am unable to use OMP parallelization effectively, for the same reason.

I was told that PetsC has a way to allocate the sparse matrix based on a rough estimate of sparsity, to begin with. Will this help? Are there libraries that can solve my issue and get help me solve the poisson equation reasonably fast?

My code is in C, I use OpenMP. I have plans to write a GPU version too.
prapanj is offline   Reply With Quote

Old   December 2, 2013, 21:03
New Member
Join Date: Mar 2011
Posts: 19
Rep Power: 6
nw_ds is on a distinguished road
How about a multigrid method (Vcycle or FMG)?? after all it's build to solve poisson equation
and it's easily paralizable whether in omp or gpu
nw_ds is offline   Reply With Quote


linear solver, particle methods, petsc, reallocation, sparsity

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Edit linear solver luckycfd OpenFOAM Programming & Development 32 January 24, 2014 14:28
Use "bounded" in scheme or not to use? immortality OpenFOAM Pre-Processing 0 June 11, 2013 16:20
suitable boundary condition for scavenging process? immortality OpenFOAM Running, Solving & CFD 3 January 25, 2013 20:10
compressible two phase flow in CFX4.4 youngan CFX 0 July 1, 2003 23:32
solver for linear system with large sparse matrix Yangang Bao Main CFD Forum 1 October 25, 1999 04:22

All times are GMT -4. The time now is 18:26.