All Questions

5 views

Max weighted subset (max sum diversification)

Given a set of elements $V$, with known cost $\pi_S$ for each subset $S \subset V$ and a monotone increasing function on the subsets $f(S)$ . I'm wondering if there is a pseudo-polynomial algorithm ...
21 views

Domain Decomposition with PETSc

Does anyone have any experience on Domain Decomposition using PETSc library? I have used PETSc for creating my vectors and matrix within my C++ code. I also used KSP to solve the linear system. I ...
21 views

Implicitly defined univariate function

So my fellow numerical computational peeps it may be that I am suffering from sleep deprivation but I'm struggling to numerically compute a function $u \rightarrow h(u)$ defined implicitly as follows: ...
19 views

20 views

Cardiac Excitation Threshold in C++ modelling

So I am trying to write a code in C++ about the cardiac excitation threshold. I know that this excitation threshold is the shortest stimulus2 value at which it can conduct an action potential (known ...
40 views

How to calculate collision force with no future knowledge

For a personal project, I am attempting to write a fairly realistic collision simulator (for relatively large objects, not quantum stuff). As I was consulting my physics textbook and various online ...
30 views

total memory usage of MPI shared memory

I am trying to use the MPI share memory feature. I have several SMP nodes, and each of them has four cores. I need an array of size N for each node that should be accessed by all four cores in each ...
21 views

In COMSOL, what is the difference between all the “Studies”? [closed]

When I right-click on (root)->Add Study, I can choose between: Eigenfrequency Frequency Domain Frequency-Domain Modal Mode Analysis What is the difference between all of them? Isn't the ...
78 views

How does the number of iteration until optimization begins depends on the dimension of the problem?

I am optimizing a function of 10-20 variables by running algorithm such as BOBYQA and a few other derivative-free algorithms. The bad news is that each function evaluation is very expensive, approx 30 ...
88 views

Condition number from incomplete Cholesky factorization

I'm having difficulties patching together from what I read about obtaining the condition number of a real, symmetric, positive definite sparse matrix. In my code, I found that there is incomplete ...
135 views

matlab: []*[] == [] makes vectorized code difficult

In my (admittedly brief) experience, there have been instances in which the code got uglier because []*[] == [], and it would be possible to write more elegant code if it was defined to be 0 instead. ...
11 views

Double integrals and parameters in GSL

I'm trying to calculate a numerical integral in C. I have a function that depends on several variables, two of which I want to integrate over while leaving rest of them as parameters. I was going to ...
44 views

Is it possible to output the matrix condition number from pardiso (MKL)? [closed]

I am assuming the pardiso solver calculates (or estimates) the condition number before proceeding to the solution phase. Is there a way to make pardiso output the condition number? Alternatively, ...
123 views

Puzzling remark about stability region of fifth-order Runge-Kutta method

I came across a puzzling remark in the paper P. J. van der Houwen, The development of Runge-Kutta methods for partial differential equations, Appl. Num. Math. 20:261, 1996 On lines 8ff on page 264, ...
29 views

Produce one smooth curve on one triangle mesh

I hope to get one smooth curve on one triangle mesh. I get one path on the mesh at first. The path consists of vertices of the mesh. I can see the path from the image below. Each one green dot ...
37 views

Derivative-free nonlinear optimization of discrete objective function with linear constraints (simplex)

I am trying to optimize a constrained-problem with a discrete, non-linear objective function. Evaluating this function is also fairly expensive. Nevertheless, despite the above two factors, I hope, ...
20 views

computational science fortran 90 solar cycle fortran model [migrated]

I'm a new fortran user. I've been given a task and I'm not quite sure how to tackle it. The task is: **Solar radiation is an energy source which plays a central role for atmospheric processes. Half ...
107 views

Test Case with Known Solution for 3D Navier Stokes Equations

Is there a test case for 3D incompressible Navier Stokes Equations like the Taylor vortex in two dimensions? I know, I can easily construct 3D manufactured solutions but I would like to have ...
52 views

Error implementing Robin boundary conditions in toy ODE problem

I am attempting to solve the following ODE problem: $$-u''+ u = x$$ $$u(0) = 0$$ $$u'(1) = -u(1)$$ The exact solution is: $u(x) = e^{-x-1} - e^{x-1} + x$ I have a Dirichlet at $x = 0$ and a Robin ...
79 views

Solving Generalization of Saddle point problem

I am interested in knowing if there is a generalization of the Uzawa iteration for the linear problems of the form \left[ \begin{array}{ccc}A& B^T&0\\ B&0&C^T\\ 0&C&0 ...
70 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 ...
78 views

Minimize quadratic form with equality constraints

I want to minimize function: $f(x) = x^T \cdot A \cdot x + b \cdot x$ given constraints: $B \cdot x = 0$. Here: $x$ is a vector ($x \in \mathbb{R}^n$), $A$ is a matrix of size $n \times n$, $b$ ...