Intel® Moderncode for Parallel Architectures
Support for developing parallel programming applications on Intel® Architecture.
Announcements
Welcome to the Intel Community. If you get an answer you like, please mark it as an Accepted Solution to help others. Thank you!
For the latest information on Intel’s response to the Log4j/Log4Shell vulnerability, please see Intel-SA-00646
1687 Discussions

## Parallel Gauss-Seidel with relaxation iterative algorithm...

Beginner
127 Views

Hello,

Description:

The Parallel iterative with relaxation method that i programmed here is
designed to be used to solve large sparse systems of linear equations
where the direct methods can exceed available machine memory and/or
be extremely time-consuming. for example the direct method of the
Gauss algorithm takes O(n^2) in the forward elimination process and
is dominated by the O(n^3) back substitution process, that means, if
for example an operation takes 10^9 second and we have 1000 equations ,
the elimination process in the Gauss algorithm will takes 0.7 second, but
if we have 10000 equations in the system , the elimination process in the
Gauss algorithm will take 11 minutes !. This is why i have develloped for
you the Parallel Gauss-Seidel with relaxation iterative algorithm in Object Pascal,
that is very fast.

And please take a look at my article on my Parallel Gauss-Seidel
with relaxation algorithm:

http://pages.videotron.com/aminer/ParallelGaussSeidel/gsrp.htm

The benchmarks here:

Please look at my parallel program gsp.pas inside the zip file , compile and execute it ... -

Language: FPC Pascal v2.2.0+ / Delphi 7+: http://www.freepascal.org/

Operating Systems: Win , Linux and Mac (x86).

Required FPC switches: -O3 -Sd -dFPC -dWin32 -dFreePascal

-Sd for delphi mode....

Required Delphi switches: -DMSWINDOWS -\$H+

For Delphi 5,6,7 use -DDelphi

For Delphi 2005,2006,2007,2009,2010+ use the switch -DDELPHI2005

Regards,
Amine Moulay Ramdane.
http://pages.videotron.com/aminer/

4 Replies
Beginner
127 Views

I wrote:
[...]
>for example an operation takes 10^9 second

I correct , i mean 10^-9 second.

Regards,
Amine Moulay Ramdane.

Beginner
127 Views

Hello,

I correct a typo:

It's the forward elimination that takes O(n^3) in the Gauss algorithm...

The Parallel iterative with relaxation method that i programmed here is
designed to be used to solve large sparse systems of linear equations
where the direct methods can exceed available machine memory and/or
be extremely time-consuming. for example the direct method of the
Gauss algorithm takes O(n^2) in the back substitution process and is
dominated by the O(n^3) forward elimination process, that means, if for
example an operation takes 10^-9 second and we have 1000 equations ,
the elimination process in the Gauss algorithm will takes 0.7 second, but
if we have 10000 equations in the system , the elimination process in the
Gauss algorithm will take 11 minutes !. This is why i have develloped for
you the Parallel Gauss-Seidel with relaxation iterative algorithm in Object Pascal,
that is very fast.

And welcome: http://pages.videotron.com/aminer/

Regards,
Amine Moulay Ramdane.

Beginner
127 Views

Hello,

I have changed my notation to something like this:

pderivative(u)/pderivative(x1)

That means partial derivatives of function u with respect to x1....

Welcome:

http://pages.videotron.com/aminer/ParallelGaussSeidel/gsrp.htm

http://pages.videotron.com/aminer/

Regards,
Amine Moulay Ramdane.

Beginner
127 Views

Phil wrote:
>the beginning of the uses statement in gsr.pas

I have updated gsp.zip with that...

Phil wrote:
>Now it runs:
>The system may or may not converge...
>The system converge...
>The system solved...
>-2.1999999779999871E+0009 -2.0910889788668756E+0009 ->2.0490340859696247E+0009

That means that Parallel Gauss-Seidel with relaxation
algorithm is working on Windows and Mac OSX also :)