Skip to main content
Stack Overflow
  1. About
  3. For Teams
Filter by
Sorted by
Tagged with
1 vote
1 answer
88 views

Summary Using the Eigen library, I see significantly different numbers of BiCGStab iterations when I solve Ax = b with a right-preconditioner, i.e., [A M^(-1)] Mx = b, (a) by precomputing A M^(-1) ...
-1 votes
2 answers
190 views

I'm attempting to implement an algorithm to compute an array of values. Unfortunately, I've found the literal implementation of that algorithm to be too slow for my use case. I suspect there's an ...
1 vote
0 answers
64 views

I have some troubles in numerically calculating PDF below: SNR of particular system model This PDF is the expression of SNR of practical RIS-aided system: https://ieeexplore.ieee.org/document/9387559 (...
1 vote
2 answers
136 views

Context I am trying to improve my ability to recognize the number of correct decimal digits after performing short calculations involving floating-point numbers. In particular, I find it challenging ...
1 vote
1 answer
286 views

From the seminal paper by Hughes I learnt that any B-Spline curve (or surface) can be described by a linear combination of B-Spline basis and control points, i.e. C(x) = \sum_{i=1}^n N_{i,p} B_i ...
0 votes
0 answers
47 views

I am moving a code base to run using PyTorch tensor on GPU. For loops are terrible on GPU especially for small-sized data. I am trying to vectorize the function below, i.e. to have no for loops. I ...
0 votes
3 answers
130 views

I'm trying to numerically solve a 2-D Poisson equation with around a 200x200 grid. I'm trying to implement the diagonal method to enable parallelism: #include <math.h> #include <omp.h> #...
0 votes
0 answers
161 views

I am studying computer science at the moment. We have the following exercise: Execute LU decomposition given a matrix A and vector b using NumPy. We had to derive L and U by hand and came up the ...
0 votes
0 answers
166 views

I'm using the mpmath library to compute a self defined function f(x) and needed to compute its higher order derivatives. I found the Automatic differentiation on Wikipedia and found it useful. # ...
0 votes
0 answers
26 views

Computing higher order derivatives is a nested function and can be quite consuming. I tried to use the parallel programming to speed up the nested function. from concurrent.futures import ...
1 vote
1 answer
406 views

I'm doing a project required to compute the precise(up to a certain precession N) numerical derivatives of a function. The usual approach was to use the finite difference types of algorithms, i.e. (f(...
0 votes
1 answer
54 views

I have C code for generic-sub-matrix-vector (gesubmv) multiplication that computes y(rinds) = A(rinds, cinds)*x(cinds) + y(rinds), where rinds and cinds are vectors of row and column (not necessarily ...
0 votes
1 answer
248 views

I want to obtain the bases for kernel (nullspace) of a matrix using the Row Reduced Echelon Form (RREF). Although scipy has the function to calculate nullspace, it does not give me sparse vectors ...
-1 votes
2 answers
716 views

I have an double integral: Because exp(−x^2) is a non-integrable function, I have tried to solve this using the quadgk function in MATLAB, but I don't get a good result. Changing the integral's ...
0 votes
1 answer
161 views

I am trying to put together a numerical simulation (specifically, Beta cell dynamics) based on Betram et al. 2007 (https://www.sciencedirect.com/science/article/pii/S0006349507709621). The model ...

15 30 50 per page
1
2 3 4 5
...
10

AltStyle によって変換されたページ (->オリジナル) /