The study of solving computational problems by exploiting the use of multiple processors simultaneously.

learn more… | top users | synonyms

1
vote
0answers
69 views

Which topics are hot in computational science? [closed]

I am accepted to computational science master program. The master program just started to accept students in our university and today each group of 3 students was assigned to an advisor teacher (we ...
0
votes
0answers
43 views

Parallel Boost Graph Library (PBGL) and ParMETIS

I want to read a graph from a text file (in parallel) and use the adjacency list data structure to do some computations on it. By computations, I mean doing parallel Jacobi relaxation method. I wanted ...
5
votes
2answers
97 views

Intel Knights Landing work loads vs NVIDIA GeForce

There are lot of articles being written about how the newly launched Intel Xeon Phi will steal the HPC\Super Computer market share from the competitors. Intel Knights is equipped with 72 cores and 4 ...
4
votes
2answers
88 views

Advice on handling many “small” matrices in parallel

I'm working on a small fun project on the side for a numerical method I've been working a bit with. Roughly, the computational problem I have to solve is the following: Assume you have a collection of ...
2
votes
1answer
75 views

Templated Numerical Linear Algebra in Parallel

I have to invert large, but densely populated matrices with higher precision arithmetic. Therefore I am looking for something like the PLASMA library, which can do Cholesky or LU factorization in ...
3
votes
3answers
235 views

Does OOP work easily and efficiently in parallel computations? [closed]

I am working on a task for which the object oriented (OO) approach fits well. I am doing it in MATLAB, since I am in the prototyping phase. However, I know that later I will surely have to perform ...
2
votes
1answer
82 views

Parallel efficiency

I would like to calculate efficiency of parallel alghoritm, using the number of computations instead of time computations. In materials from my studies I have a formula like below: $$ \eta(n,p) = ...
3
votes
0answers
100 views

What are the career options for computational scientists other than being in software Industry?

I would like to know the career options for a computational scientist other than getting into software industry. I wish to use the my skills of Numerical Linear Algebra,Parallel Programming, ...
2
votes
2answers
193 views

GPU vs CPU calculation

I've been working on calculating large factorials ($N>10^9$) and I was wondering if it wasn't faster to use the GPU to run the calculations on something like openCL. What I realized however was ...
1
vote
1answer
96 views

Is there a single C++ IDE packed with necessary scientific libraries for computing purpose?

I am basically looking for a good open source C++ IDE with scientific and visualization libraries built in. Can any one suggest few projects?
2
votes
0answers
67 views

preconditioned Uzawa method with Petsc

I am trying to improve the resolution of a Stokes problem (P2/P1 on unstructured mesh) defined by the matrix $M$: $M= \begin{pmatrix} A_u & 0 & B_u \\ 0 & A_v & B_v\\ B_u^T & ...
5
votes
2answers
235 views

How to Run MPI-3.0 in shared memory mode like OpenMP

I am parallelizing code to numerically solve a 5 Dimensional population balance model. Currently I have a very good MPICH2 parallelized code in FORTRAN but as we increase parameter values the arrays ...
4
votes
0answers
73 views

Computing entries of a matrix in parallel using GPU

I am trying to solve $Ax=b$ a system of linear equation. The matrix $A=(a_{ij})$ is constructed as $$a_{ij}=\int_0^1\int_0^1 K(x,y)h_i(x)h_j(y) dxdy.$$ The evaluation of integral takes a lot of time, ...
3
votes
1answer
123 views

Challenges in implementing Algebraic Multigrid on millions of processors

I just implemented an Algebraic Multigrid solver for a Mixed Dirichlet-Neumann Boundary Value problem and was surprised to see the speed-up as compared to a simple iterative solver for a large problem ...
0
votes
1answer
115 views

Parallel slower than serial

Issue Hello everyone, I have got a program (from the net) that I intend to speed up by converting it into its parallel version thru the use of pthreads ...
1
vote
1answer
62 views

C++ - Starting paraller instances of an executable from a C++ code

Can I start parallel run of LAMMPS from inside a C++ code. Is it possible to compile my code and then do mpirun with my compiled code. I also have one system command inside the code, which starts ...
2
votes
2answers
209 views

GPU-enabled Lattice Boltzmann solvers?

Is anybody aware of any GPU-enabled Lattice Boltzmann solvers (preferably on C++/OpenCL and open-source) that would be recommended? I have found Advanced Simulation Library, but it seems to be very ...
1
vote
2answers
389 views

How to use the basic Sparse matrix operations (multiplication, .etc) in PyCUDA

I try to use sparse matrix operations in GPU in Python and now try to use PyCUDA with theano. But I can't find how to do sparse matrix and vector multiplication. I only got an example showing how to ...
2
votes
1answer
122 views

Why is MPI/PETSc only showing 1 processor even though it works for helloworld?

I am a beginner with Unix and MPI. For the serial program, I run the program by typing make in /directory/buildagain and then ...
2
votes
1answer
52 views

Performance based on the Roofline model

To somewhat follow up on the question asked here, I have been told that the Roofline model is one way of assessing the performance of any scientific code. Basically I compute the Arithmetic Intensity ...
0
votes
0answers
58 views

How come in distributed computing 10X 32GB systems doesn't equal 1 320GB system

is there a mathematical formula to relate how much memory per node is lost when you use parallel computing? for example, when running COMSOL (FEA, MUMPS) in parallel, a problem that would require ...
1
vote
2answers
140 views

Parallelization of element-wise matrix multiplication

I use Armadillo as an interface to OpenBLAS. In my current program, I have a loop, in which I do multiplications of the form ...
4
votes
2answers
141 views

Comparing various implementations/software packages for large-scale finite element simulations

I currently use FEniCS and Deal.II to solve various FEM problems. I am also writing my own implementation of these problems by directly implementing the data structures, routines, and solvers within ...
9
votes
4answers
2k views

Parallelizing a for-loop in Python

Are there any tools in Python that are like Matlab's parfor? I found this thread, but it's four years old. I thought maybe someone here might have more recent experience. Here's an example of the ...
2
votes
1answer
86 views

On an HPC machine, does one compute node equate to one motherboard?

I am using the Hydra Process Manager from MPICH and I plan on running jobs on multiple compute nodes. Currently working with an Intel Xeon E5-2680v2 with a total of 20 cores with 2 sockets (32 GB of ...
1
vote
0answers
31 views

Minimum effort merging of two sets

I have the following problem. I have two sequences of elements $A = [a_1,a_2,\cdots,a_n]$ and $B = [b_1,b_2,\cdots,b_m]$. I can build a matrix $D[n \times m]$ where $d_{ij} = d(a_i,b_j)$ My greedy ...
0
votes
1answer
181 views

Armadillo eig_sym() for extracting eigenvalues. Is it parallel at all? [closed]

After wasting 3 days with scalapack, I gave up and moved to Armadillo, considering it uses lapack underneath its beatiful and easy interface. I would like to calculate the eigen values and eigen ...
0
votes
0answers
130 views

parallel computatioan of a PDE in MATLAB

I want to solve a 1-D PDE $(\partial_{tt} + \alpha\partial_t)u(x,t)=\partial_{xx}u(x,t)-\sin(u(x,t))+f$, using method of lines and for this I defined a spatial grid of about n~1000 points. Since my ...
0
votes
1answer
78 views

Already exist? parallel computing [closed]

Parallel computing allow us to compute tasks in parallel. And if there was a system like torrent that unlike sharing files, it share computing power? For example I have a gpu and I use only 25% of its ...
2
votes
1answer
70 views

communication penalty when using wide stencils in parallel computations

When reading about discontinuous Galerkin methods one finds the argument that these methods allow higher-order accuracy while maintaining a compact stencil (a cell only communicates with its direct ...
0
votes
0answers
88 views

FFT parallel processing in MPI

I am working now in Beowulf Cluster and parallel processing, I want code for Fast Fourier transfer functions written in any language, e.g., C/C++. Without using FFTW library based on message passing ...
0
votes
0answers
170 views

Advantages and Disadvantages of PETSC vs HyPre?

What are some of the differences between using Petsc and using Hypre? Also what are some of the advantages and disadvantages of both? Does Petsc use more memory or run faster? Thank you very much for ...
0
votes
3answers
197 views

Determining efficiency in MFLOPS/s of a parallel program

I am running some scientific (parallel) code and would like to obtain some performance profiling measurements. I want to obtain the "efficiency" of the code in terms of flops/s over theoretical (peak) ...
1
vote
1answer
69 views

Shared memory parallel computing solutions compared

As of 2015, what shared memory parallel computing solutions are available? What are the advantages and disadvantages of each for various use cases arising in high performance scientific computing? I ...
1
vote
1answer
77 views

Parallel Monte Carlo simulation using PETSc

I am trying to do Monte Carlo simulation for a large problem which requires eigensolution of a matrix for each sample. The matrix itself is quite large so much so that I want the eigensolution itself ...
2
votes
1answer
171 views

Efficient Successive Over Relaxation using Cuda?

I recently implemented Successive Over Relaxation using Cuda as a part of my course project and was curious to know how I can make the code more efficient. I'm using Red/Black SOR scheme which is ...
4
votes
1answer
143 views

Is solving Verlet constraints parallelizable?

I am a computer science student and I came across the Verlet Cloth simulation. An implementation is here. I was interested in how such computation would lend itself to parallelization. The problem ...
1
vote
1answer
382 views

Parallel vs Serial Thomas Algorithm

I am currently writing a code that solves a large tridiagonal matrix every iteration and runs for 1,000's of iterations. I am currently using a Thomas algorithm to solve the matrix serially. I found a ...
3
votes
1answer
198 views

Designing a preconditioner for a very Ill-conditionned matrix

I am a physicist with limited numerical methods knowledge and I am trying to speed up the inversion of a very ill-conditioned problem ($rcond>10^{30}$). The same sparse square matrix is used ...
1
vote
2answers
192 views

Are there any popular (paralleled) implementations of Lanczos methods for SVD/eigendecomposion?

I want to use it in Matlab or Java. Will these two languages be much slower for computing the algorithm compared to C, C++, in case efficiency is an important factor? I'm aware of that there's a ...
1
vote
1answer
73 views

Efficiently swap vectors in PETSC

The PETSC function VecSwap after a few layers of abstraction calls the BLAS primitive swap, which essentially copies data around (unless I am missing something). ...
1
vote
0answers
61 views

Any book or teaching material on computational neuroscience using hardware acceleration?

I am looking for some material (be it in the form of a book, reader notes, teaching material freely available etc.) on the topic of computational neuroscience employing hardware acceleration? I guess ...
4
votes
2answers
615 views

Parallel optimization algorithms for a problem with very expensive objective function

I am optimizing a function of 10-20 variables. The bad news is that each function evaluation is expensive, approx 30 min of serial computation. The good news is that I have a cluster with a few dozen ...
3
votes
1answer
501 views

Solving Advection (Convection) - Diffusion - Reaction Partial Differential Equation in Python

I am looking for library written in Python which will enable me to solve the coupled nonlinear equations which looks like: I need the library which will enable me to couple this solver to other ...
1
vote
3answers
201 views

Parallel solver for sparse matrices on unstructured grids

I am trying to solve Euler equations on unstructured grids. Consequently, the problem reduces to solving Ax=b where A is a ...
2
votes
2answers
461 views

Using Trilinos and/or Petsc in Windows with Visual Studio?

I am currently working on a project of building a simulator. Up till now I have mostly played around with very simple examples, but with regards to the future I have been thinking that it might be a ...
4
votes
1answer
260 views

Parallel Gram-Schmidt algorithms

I've heard that classical Gram-Schmidt is more amenable to parallelization than modified Gram-Schmidt; apparently the reason has something to do with "level 2 BLAS", which I'm not familiar with. ...
2
votes
1answer
118 views

Optimal Discontinuous Galerkin (DG) solver on a parallel system

I am seeking the optimal method for implementing DG on a parallel system. For my research, I come across two types of problems. For the first problem, I am solving a time-independent (steady-state) ...
0
votes
1answer
329 views

Using OpenMP in Fortran for common array operations

In Fortran, given arrays A(1000000) and B(1000000), and to compute the sum you simply write: C = A + B However, when I want to implement OpenMP, I have to write ...
1
vote
1answer
177 views

Fast way to compute all eigenvalues of a dense Hermitian matrix

I am finding the eigenvalues of dense NxN Hermitian matrix which is calculated from a density operator in quantum physics. All the eigenvalues are needed as I need to calculate the sum of the absolute ...