All Questions

Filter by
Sorted by
Tagged with
-1
votes
0answers
20 views

What is the meaning of triangles color in the result of Tipping Problem in scikit-fuzzy (fuzzy logic)?

I am following this example https://scikit-fuzzy.github.io/scikit-fuzzy/auto_examples/plot_tipping_problem_newapi.html from documentation of scikit-fuzzy library,but I have a question in the figure ...
1
vote
0answers
17 views

Is it possible to use a fixed point iteration for solving this nonlinear system?

Consider the following differential equation \begin{align} \frac{\partial f(u)}{\partial x} &= g(x), \ \ x\in [x_{L},x_{R}] \label{Eq2.2} \\ u(x_{L}) &= g_{1} \\ u(x_{R}) &= g_{2} \end{...
3
votes
2answers
403 views

Is there any open-source code for a hybrid 2D mesh (triangles and quadrilaterals)?

The question is pretty much the title. Note that I have lots of experience using open-source meshing tool, e.g. Gmsh and OpenFoam blockMesh & snappyHexMesh. Nevertheless, I have no idea on how to ...
2
votes
0answers
11 views

Is there any reliable free/open source tool for structured mesh smoothing?

I have been using Pointwise for grid generation and found the quality of smoothed grids to be stunning. I am not aware of any free/open source alteranative that offers the same capabilities for ...
2
votes
2answers
78 views

Is there a Python version of the ODE tool pplane?

This is the same question as this one, except for Python instead of Mathematica. Basically, the MATLAB software pplane is a staple in ODE courses. Is there a Python equivalent? Sample outputs from ...
1
vote
0answers
15 views

How do you correctly implement Scipy's FFT procedures to produce a low-pass filter - image processing

I'm following this low-pass filter example in the text "Image Operators: Image Processing in Python 1st Edition" by Jason M. Kinser, but can't seem to duplicate their results. The text's ...
11
votes
2answers
885 views

Higher precision floating-point arithmetic in numerical PDE

I have the impression, from very different resources and talks with researches, that there is a growing demand for high precision computations in numerical partial differential equations. Here, high ...
2
votes
0answers
27 views

3D Cooley-Tukey FFT

To compute the $N$-point DFT $$ X[k] = \sum_{n=0}^{N-1} x[n] W_N^{kn} $$ where $N = N_1 N_2$, we can write the indices as $n = N_2 n_1 + n_2$ and $k = k_1 + N_1 k_2$, (effectively packing the data ...
3
votes
1answer
336 views

When do not use preconditioners for sparse linear system of equations?

I'm implementing a solver of Finite Element Method, and to solve the linear system of equations I'm using gmres from MKL of Intel. Exists the option with and without a preconditioning. In what case it ...
4
votes
1answer
74 views

Roundoff errors in FEM computations - generalized eigenvalues

This is a continuation of my previous question. I am trying to effectively compute a bound for the roundoff errors in some FEM computation (2d polygons, triangular meshes). Below I will write some of ...
5
votes
1answer
101 views

Algorithms for computing winding numbers of 2-sphere maps

I have a question concerning computational geometry which arises in the simulation of fields with topological defects, and I'd like to know whether there's an efficient algorithm (or any algorithm) to ...
6
votes
0answers
72 views

FEM : energy minimization VS PDE solving

Engineering FEM When I studied engineering, I learned the traditional approach for finite elements for elasticity. The point was to solve the PDE $-div(\sigma)=f$ as: Multiply your PDE with a test ...
1
vote
0answers
21 views

2D DFT for lower frequencies only; is there something significantly faster than numpy.fft.fft2 (throwing away high frequencies)?

I do a lot of 2D discrete FFT in python using np.fft.fftshift(np.fft.fft2(y)), then throw away 90% or more of the array, keeping only the central low-frequency area....
0
votes
1answer
172 views

Solve non-linear equation in R

I need to solve the following equation for $x$ in [0, 1]. Assume $0<\alpha<1$ and $0<\lambda$. $$(1 - x)^{\alpha+1} - \lambda (x+1)^{\alpha+1} = -2\lambda (\alpha + 1) x^\alpha$$ Would very ...
2
votes
1answer
168 views

How to avoid gsl root finder evaluate function outside its domain

When I use the newton's method or hybrid solver in the GSL package to deal with 1-D or multidimensional root solving problems, the code frequently crashes when the solver requests function value ...
0
votes
1answer
283 views

Global to local transformation matrix in 2D frame structures

In section 3.2 of this paper [1], where 2D planar frame structures are being analyzed, the authors mentioned a transformation matrix to be used in extracting the element displacement vector from the ...
2
votes
0answers
70 views

Regularisation of ill-conditioned matrix-vector problem

I have a linear* problem which arises from an integro-differential system, and writes: $$ (\mathbf{I}+\lambda \mathbf{A})x = b $$ where $\mathbf{A}$ is a real full matrix, size $n\times n$, but is not ...
8
votes
4answers
3k views

What can a computational scientist do in the fourth industrial revolution?

This question is neither scientific nor technical but more career related. I am at a junction in my professional life where I need to make a decision with regard to the future of my career. At the ...
0
votes
1answer
71 views

Computing eigenvalues of Schrodinger equation with spin

I want to solve a 2-dimensional particle in box problem with two electrons in the quantum well.I would like to take into account spin of electrons and Coulomb interactions to compute singlet and ...
2
votes
1answer
255 views

Is there a simple way to avoid carbuncles for FD WENO methods?

I have implemented finite-difference WENO scheme for Euler equations (with some variants - WENO-JS, WENO-Z, WENO-M, different flux splitting). It works well, but have problem with so-called carbuncles ...
1
vote
0answers
112 views

How to compute the Eigenvalue and Eigenstates of Quantum well with Effective mass using finite difference method in Python?

I want to compute the eigenvalues and eigenstates of a quantum well with different effective masses of electron in the barrier and in the quantum well. As can be seen [1]: https://github.com/mholtrop/...
-1
votes
0answers
34 views

Null Christoffel symbols associated to the FLRW metric obtained via Mathematica [closed]

I'm trying to make a mathematica notebook that computes the Christoffel symbols associated to the Friedmann-Lemaître-Robertson-Walker metric (noted FLRW in the notebook) describing an homogenous and ...
1
vote
1answer
1k views

Split-step Fourier method applied on Schrodinger equation

I'm trying to solve a Schrodinger equation of the form $i\frac{\partial}{\partial t}\psi=-\frac{\partial^2}{\partial x^2}\psi + (V(x)+\alpha|\psi|^2)\psi$ using the split-step Fourier method ...
1
vote
2answers
55 views

Two-dimensional ordering issue – alternate sort order ascending/descending to reduce fluctuations - trivial?

I have a solution in search of a problem that some of you could perhaps help me with. Let $L$ be a list of elements. Each element has two inherent properties/attributes ($a$, $b$) that can each be ...
3
votes
2answers
545 views

Time Reversibility of Velocity Verlet Algorithm

I'm very new to computational Physics and am finding conflicting statements on whether the velocity Verlet algorithm, defined as: $\begin{align} x_{n+1} &= x_n + v_n \Delta t + \frac{1}{2} a_n \...
3
votes
0answers
45 views

What determines the order of a finite volume scheme?

I often hear that cell centred finite volume is second order accurate but at the same time I come across notions of high order FVM flux schemes. Is there a distinction between the two? If I were to ...
2
votes
1answer
89 views

Determining the voxels between two boundary surfaces

Issue description I am working on human brain tACS simulations where I have the models of the skin, skull, csf, brain and ventricles in STL format. The shape does not matter and there are no ...
1
vote
2answers
110 views

Why aren't face integrals for an element calculated in FEM but they show up in FVM?

Consider the Laplace problem: \begin{align} -\nabla^2 u = f \qquad \text{in } \Omega \\ u = 0 \qquad \text{on } \Gamma \end{align} The weak problem is find $u_h \in V \subset H^1$ such that $\...
0
votes
0answers
31 views

Help with debugging block GMRES

I have written block version of GMRES by referring [1] and MATLAB implementation of gmres. I need to write it for complex matrices. My block implementation when run on single RHS is giving correct ...
2
votes
0answers
185 views

Implementation of Z^2 error estimator in Abaqus for adaptive mesh refinement

Currently, I am working on a remeshing routine for my simulations (Abaqus 6.14-1) using python scripts. The simulation deals with the Brinell indentation test and as the remeshing software I use Gmsh ...
0
votes
0answers
42 views

State change of input-output system

Edited Given a computer model $F:\mathbb{R}^3 \to \mathbb{R}$, with inputs $x, w$ and $z$, and output $y=F(x,w,z)$, where for any input, we are able to evaluate the output, my goal is to tune the ...
1
vote
3answers
106 views

Fitting line to a staircase function

I have a staircase/step function $n(E)$. I know the points $\{E_i\}$ at which each "step" occurs and all steps are of constant height 1. I need to fit a line $a + bE$ to this function and ...
0
votes
2answers
262 views

Is it possible to partition 2D data into bins such that each bin contains the same number of samples?

I am trying to sort data following a bivariate distribution into a numpy histogramdd, where each bin should contain the same number of data points (to the nearest whole sample). I expect that some ...
4
votes
2answers
72 views

Backward Euler + Quasi Newton(Broyden) method fails to solve Van der Pol's equation(Stiff ODE)

The first guess is using the forward Euler approach. The first jacobian is using finite differences. Then NR method is used to solve for the next iteration and Broyden's method is used to update the ...
-1
votes
0answers
17 views

If L is regular so is collapsing doubles [closed]

Suppose L is a regular language. I need to prove that $$ L'=\{x_0\cdots x_n:x_0x_0x_1x_1\cdots x_nx_n\in L\}$$ I thought I could take a DFA which computes L, and take each accepting state, together ...
7
votes
3answers
1k views

Tanh-sinh quadrature numerical integration method converging to wrong value

I'm trying to write a Python program to use Tanh-sinh quadrature to compute the value of \begin{equation} \int_{-1}^1 \frac{dx}{\sqrt{1-x^2}} \end{equation} but although the program converges to a ...
2
votes
2answers
138 views

Different sources of error in Finite Element computations

Consider the problem $-\Delta u = f$ in $\Omega$, with $u=0$ on $\partial \Omega$. Suppose that $\Omega$ is a polygon and that we approximate the solutions to the previous problem using Lagrange ...
1
vote
1answer
45 views

Importance Sampling for multidimensional integrals and random numbers from multivariable pdf's

I am aiming to get a numerical value for a five-dimensional integral using Monte Carlo Integration. I am getting good results using the Mean Value Method, but I would like to try to use Importance ...
6
votes
1answer
132 views

General approach to infinite sums

My question is specific to algorithms and models of computation. I would like to write code to evaluate the following expression quickly and accurately: $$\log \left( \sum_{i=1}^{\infty}{I_{\nu+i}(2\...
3
votes
0answers
64 views

The implicit form of a NURBS curve

I am trying to evaluate and analyse a NURBS curve to generate a mechanism. I understand that the general form of a NURBS curve is commonly written as a parametric equation in the form of $f_{par}(t)$. ...
-1
votes
0answers
20 views

how to unscramble a .wav file to find the actual voice?

this question is a bonus question for my HW, and it's supposed to be something I am suppose to look up to solve. However, I have no idea what to look up on here to help me write this code. Can anyone ...
1
vote
0answers
30 views

Digital beamforming: how payload manipulation can change beam direction without manipulating the carrier?

I'm interested in how digital beamforming works and I can not find an answer for a lot of time. I googled, asked teammates, and couldn't get it. Let me describe my question. From my understanding, the ...
0
votes
0answers
36 views

Normalizing the right-hand side in Jacobi-preconditioned conjugate gradients

I have been reading the following paper: CG versus MINRES: An empirical comparison. In it a conjugate gradient solver is applied to a system matrix $A$ Jacobi-preconditioned on both sides. ...
-2
votes
0answers
47 views

How to run PGI cuda on an M1 Macbook? [closed]

I am trying to run a piece of code that was written in PGI cuda FORTRAN. I usually run it on a cluster that has all the packages and software already installed. I would like to run this code on my M1 ...
3
votes
1answer
101 views

Multi threaded finite element assembly implementation

What is typically the best way to multi thread the assembly loop in a finite element code? Does anyone have experience with implementing this, that they can share? I can think of a couple of ways of ...
1
vote
0answers
33 views

Is there a way to generate a sample $(X_i, Y_i, Z_i)$ from custom distribution?

I'm newbie here. I'm wondering if it's possible to generate $(X_i, Y_i, Z_i)$ from my own distribution function? I know that there is a way to make own class for 1D variable. But what about 3D case?
10
votes
3answers
631 views

Are there any libraries out there that implement block Krylov subspace methods?

Question Are there libraries out there that implement block Krylov subspace methods? (I was not able to find any from a simple Google search.) Background Right now, I am working with a code that ...
0
votes
1answer
172 views

surface: rows (Z) must be the same as length (Y) and columns (Z) must be the same as length (X) in octave

I am having a dataset and trying to plot a 3-d plot between the independent and dependent variables.but, I am getting this error whenever I am trying to plot here is my code: ...
0
votes
1answer
54 views

Manufactured solution to 2d convection-diffusion with homogeneous Robin boundary conditions

I am looking for a manufactured (or analytical if it exists) solution to the 2d boundary-value problem $$\frac{\partial u}{\partial t} = \mathbf{a} \cdot \nabla u + D \nabla^2 u \quad \quad \mbox{in } ...

15 30 50 per page
1
2 3 4 5
198