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

learn more… | top users | synonyms

0
votes
1answer
41 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
64 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
71 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
67 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
44 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
51 views

Computational problems needed

I'm developing a new distributed computing solution. I'm looking for problems that need a lot of computational power in order to validate the technology and get the first use case. We're currently ...
0
votes
0answers
74 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
93 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
56 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
64 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
92 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 ...
3
votes
1answer
106 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
179 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
158 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
126 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
64 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
46 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 ...
3
votes
2answers
265 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 ...
2
votes
1answer
209 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
144 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
185 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 ...
3
votes
1answer
181 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
89 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
140 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
91 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 ...
1
vote
1answer
91 views

Solving a problem using multiple threads using CLP

Is it possible to run Clp on multiple threads using Julia and JuMP?
3
votes
1answer
101 views

How does matrix-matrix product scale with multiple CPUs?

These days, one can have 64 cores in a single node. I wonder how well the dense matrix-matrix product (SGEMM and DGEMM) scales ...
4
votes
1answer
343 views

Python: Multistart optimization using parallel programming

[I previously asked my question on StackOverflow but this site may be more appropriate] In Matlab, I am currently using the MultiStart as an optimization algo in a parallel setup for a computer ...
2
votes
1answer
188 views

Maxwell vs Kepler for GPU computation

I am looking at the Nvidia GT-860M which comes in both the old Kepler and new Maxwell architecture. The old one (1152) seems to have almost twice the cores as the new one (640). The new one has a ...
3
votes
2answers
183 views

C - OpenMP, MPI, Serial Program

I'm part of a Computational Science course and come from a non-programming background, so please forgive me my ignorance. I'm working on a set of code in C to numerically solve the Navier Stokes ...
3
votes
2answers
77 views

Local versus remote development

I know there are similar questions asked on other StackExchange sites, but I think development in Computational Science may work differently because of the computing environments that are common. I ...
3
votes
1answer
75 views

Threaded QR with column pivoting

My program needs to perform pivoted QR decomposition on tall (e.g. 1e9 by 100) matrices. I run into the bottleneck that the major computational time of my program is spent on doing serial pivoted-QR ...
1
vote
0answers
80 views

Why would this a pair of E5-2670 cpus be faster than 4 E5-2640 cpus for large matrix diagonalisation problems?

Not really much more to say; Of several computers available for use, these are two of the larger ones; one has 2x E5-2670, and the other has 4x E5-4640. The problems we're looking at essentially boil ...
4
votes
4answers
121 views

Removing non-determinism from molecular dynamics code

I've been looking through the other answers, and I haven't found a good answer yet. I'll describe the simulation I am running, and then the problems. The program simulates particles undergoing ...
1
vote
0answers
118 views

How to solve sparse linear systems using CPUs and GPUs?

There are several libraries to (iteratively) solve large sparse linear equation systems in parallel on a number of CPUs. Our parallel cluster also has attached powerful GPUs, but so far, I did not ...
1
vote
3answers
227 views

Scripting for High Performance Computing

I would like links to resources which provide codes for scripting (preferably in python) for high performance computing. I'm looking for sample codes and examples which can, (1) automatically create ...
1
vote
1answer
105 views

Is it possible to simulate fluid dynamics in a time-based and deterministic manner?

The Problem Domain I have a number of network-connected PCs. I want to be able to simulate and replicate the same simple fluid dynamics simulation (E.g. Navier-Stokes), in real-time, between them. ...
4
votes
2answers
1k views

GPU Computing: Which GPU Brand and Type gives best value for money?

Recently I started to get interested in GPU-computing. I am considering buying a GPU and am inclined to buy an NVIDIA GPU since this allows for coding in both OpenCL and CUDA. But is NVIDIA really the ...
2
votes
3answers
202 views

Dense distributed matrix

A dense matrix is distributed for parallel computation column-wise, then multiplied from left & right by sparse matrices. What would be appropriate c++ libraries for these tasks?
3
votes
2answers
271 views

Best incremental multidimensional Delaunay tessellation algorithm

I'm looking for a specific type of Delaunay tessellation algorithm. The algorithm should be: incremental so that I can add new sites inside known simplexes (i.e. no searching for the right simplex ...
2
votes
1answer
146 views

Installing GPAW (Python DFT package) on a computing cluster

I'm in the middle of installing GPAW on a cluster but I'm running into problems building the "gpaw-python" interpreter, its job is to allow users to execute parallelised GPAW calculations. Following ...
3
votes
2answers
213 views

Why is computational cost measured in Floating Pt. Ops. in times of parallel computing?

In times of parallel computing, it seems to me that algorithms (also basic ones, like matrix-vector multiplication) should be measured by their dependent steps (that use results from steps before) ...
3
votes
0answers
58 views

BLACS context value and multiple MPI communicators

I am reposting here a question previously asked in stackoverflow: as suggested in a reply to my question, scicomp could be a better place to obtain some useful comments/suggestions for my problem. I ...
7
votes
5answers
276 views

Is discontinuous Galerkin really any more parallelizable than continuous Galerkin?

I've always heard that easy parallelization was one of the advantages of DG methods, but I don't really see why any of those reasons don't also apply to continuous Galerkin.
3
votes
3answers
248 views

Leave openMP theads ready

I have a C++ program that has a double for loop that I'd like to parallelize using opeMP. It is simple (a polynomial product) of varying length (between 100 and ...
3
votes
1answer
127 views

Increasing thread number degrades performance

I have a problem with parallelization and the brownian dynamics (molecular dynamics) code that I am using. We have our own home-grown framework at the university, and recently we've made the change ...
2
votes
1answer
85 views

Optimizing MPICH2 on FX 9370 processor

Short version of question: We have a Opteron 6272 cluster, an we have recently started bilding FX 9370 based small cluster. Problem is that we are not able to get desired speed, i.e scaling with ...
5
votes
3answers
278 views

Does IBM LSF facilitate parallel computing?

I have been browsing through the documentation of the IBM Platform Load Sharing Facility, to find out if it supports parallel task computing. Below is what I think is relevant for this question, ...
4
votes
3answers
192 views

Is there a general framework for solving PDEs on uniform grid in parallel

Hej, I want to simulate a partial differential equation (a modified Cahn-Hilliard equation, but the details do not matter much. The questions also applies to the diffusion equation). I'm looking for ...
2
votes
1answer
104 views

Multiple independent random number streams

Having multiple streams of pseudo-random numbers known to be independent and with a uniform distribution I want to do Monte Carlo simulations in parallel. In other words, one thread will have a ...