All Questions

Filter by
Sorted by
Tagged with
0
votes
0answers
47 views

Determining the indices of a VTK mesh

I am trying to assign fixed constraint at specific indices of VTK mesh, however, I only can view STL file in a blender as follows: Upon assigning the favored indices in my scene at the SOFA physics ...
0
votes
0answers
18 views

I have a trouble determining appropriate transition matrix

M.Sc students of the Department of statistics, FUTA are expected to do course work for a year and write their thesis the following year before graduating. A student has a probability of 0.25 of ...
0
votes
0answers
32 views

solving electron motion in undulator by Boris method

I am trying to use Boris scheme to solve the electron trajectory in undulator. The undulator field I used is: $$B_x = b_0\sin(2\pi \tfrac{z}{\lambda_u})$$ where $b_0 = \dfrac{2\pi c_{0}K}{q m_{e} \...
0
votes
0answers
33 views

Method for implementing QP solver with matrix terms?

I am trying to implement (in code) a QP solver for the following equation: $$\min_{u} u^{T} Wu$$ $$s.t. \; \beta u = \tau_{ref}$$ $$ Au \leq b $$ See this document, section 5.1 (Page 35) $u$ is a ...
0
votes
0answers
80 views

Smoothing FFT result

I am trying to calculate the spectrum of Bremmstrahlung, which involves calculating the Fourier transformed acceleration. I am solving a non-linear ODE to numerically calculate the acceleration in the ...
0
votes
0answers
38 views

Divergence issues when using intrinsic cohesive elements approach

When I model the strain localisation of a microscopic sample (or say RVE ) with cohesive elements approach, the convergence performance looks very terrible. I have to use extremely time increments (...
0
votes
0answers
172 views

What algorithm do BLAS and ATLAS use for matrix multiplication?

I have searched and what I understood was that they use the naive one with several memory and cache optimizations. However, I wanted to know whether they are using the Strassen or the Coppersmith-...
0
votes
0answers
23 views

Controllability - Maximum Matching

I found this image on wikipedia referring to Barabasi's work on Network Controllability. I tried to verify it. We have a A matrix of dimensions (20 by 20) made as the image suggests. According to the '...
0
votes
0answers
37 views

host vs nodes vs sockets vs thread vs processors vs cores in cluster computing [duplicate]

I am new to the computer cluster and I am so confused with these terms, could you please explain them? Regards, Vricc
0
votes
0answers
19 views

How to start coding for posterior inference

I am trying to implement the model given in http://proceedings.mlr.press/v84/andersen18a/andersen18a.pdf where they have used mean-field variational inference for posterior inference, but I want to ...
0
votes
0answers
67 views

FEniCS implementation of Maxwell equations for a dipole antenna

someone knows where I can find a FEniCS implementation of Maxwell equations for a dipole or other type of antenna? I mean a dipole antenna with an arbitrary geometry of every 'leg' in the dipole.
0
votes
1answer
83 views

The proper way to assess the error of Jacobi iteration (for 2D Poisson equation)?

Motivation: I'm using 2D regular grid (it's actually a quadtree but I can still treat it as a finite difference thing if I weight-average the solution over smaller scale cells for the purpose of ...
0
votes
0answers
28 views

On using Ritz Method to solve a Mindlin–Reissner plate

I am trying to replicate the method given in the this paper. I have written a Matlab program which determines the displacement field of Mindlin–Reissner plate theory using Ritz method. The limitation ...
0
votes
0answers
33 views

Grids for atmosphere simulation with finite volumes on the globe

I am currently in the early construction process of building a simple CFD model of a rotating planetary atmosphere. The planet should be allowed to tilt significantly, so that a time-dependent source ...
0
votes
0answers
70 views

How one can simulate a system given by differential equation?

I want to simulate a diffusion environment given by the differential equation $$\frac{\partial u(x,y,t)}{\partial t}=D\left(\frac{\partial^2 u(x,y,t)}{\partial x^2}+\frac{\partial^2 u(x,y,t)}{\...
0
votes
0answers
66 views

Reducing run time of a numerical calculation using a mex file in Matlab

I wrote a Matlab code that involves doing a numeric calculation (relaxation), but it is quite slow. I learned of the possibility of using a mex file to run a C code and integrate it into Matlab, so I ...
0
votes
0answers
33 views

Coupling 1D-3D problem in FENICS

I was wondering whether any of you knew if FENICS has the capability to solve coupled 1D-3D problems that are linear non-iteratively? As an example, pipes embedded in a porous domain will provide ...
0
votes
0answers
68 views

Can I use matplotlib to plot the surface of a 3D body?

If Matplotlib could volume render, I would not ask this question. But it can't. Can I however instead use Matplotlib to plot the surface of a 3D body? I.e. is there a way to (i) triangulate the ...
0
votes
1answer
65 views

Dealing neighbor list in NVT Monte Carlo (MC) simulation

I'm making a NVT Monte Carlo (MC) simulation code with only short range interaction. I found many MC tutorial codes (usually Lennard-Jones system) in online. However, most of them are doing energy ...
0
votes
0answers
123 views

How to do a Hoshen–Kopelman algorithm in Matlab with periodic boundary conditions?

Say I have an array s with elements of 0 and 1. From my understanding, a Hoshen–Kopelman algorithm for s will be (in 3 ...
0
votes
0answers
45 views

Dividing Point Cloud into voxels

I am reading a paper on implementing convolution neural network for a 3D point cloud. In this paper, they are dividing the point cloud into voxels. Is there any easy way to do it using point cloud ...
0
votes
0answers
37 views

Is semi-Lagrangian 1D advection identical to upwind Euler, when $|u|<\Delta x/\Delta t$?

This looked true to me, I worked through the algebra for 1D and confirmed it. But no one seems to mention it, so maybe I'm missing something... I'm using Bridson's SIGGRAPH 2007 course notes. [5MB ...
0
votes
1answer
35 views

Change of random variables and check by plot

Question As a test, I transform a uniform distribution over the unit square. But when I check the transformed distribution with Monte Carlo, it is wrong. What went wrong? Thanks. Problem Random ...
0
votes
0answers
48 views

Simulate friction arround steady-state point

When simulating a friction dynamic system, I have come to a problem when approaching steady-state. In the example below, a box moves along a surface with a certain friction coefficient: ...
0
votes
0answers
57 views

A non linear ode with boundary conditions at infinity

I want to solve the non-linear ODE $$\frac{d^2}{dx^2}y=a(y+y^3)$$ With the boundary conditions that $$\lim_{x\to \pm \infty} y(x) =0$$ I am not aware of any analytical method for solving this kind ...
0
votes
0answers
117 views

Shape functions in Euler Bernoulli Beam Equation

Does anyone have a intuitive explanation of why Hermite polynomials have to be utilized as the shape functions in the FEM solution of the Euler Bernoulli Beam 4th order ODE? I have been learning FEM ...
0
votes
0answers
67 views

Why does GMRES converge much slower for large Dirichlet boundary conditions?

I'm trying to numerically solve a simple Laplace equation in 2D, with a nonlinear source term: $\nabla^2 u = u^2$ with boundary conditions as $u=0$ everywhere except for $y=1$ where $u=u_0$. I'm ...
0
votes
0answers
79 views

How to read data from simply formatted text file (c++)?

Mesh information like points, faces and cells is to be stored into separate files: e.g. for points file: # points data: x y z N_Points 100 x1 y1 z1 x2 y2 z2 ... cell file: # cells data: ...
0
votes
0answers
57 views

What are the differences between these different forms of equation?

What are the differences between Conservative differential form, Non-conservative differential form, Conservative Integral form and Non-conservative integral form of differential equations? I know ...
0
votes
0answers
33 views

A preconditioner for self-consistent iteration

I tried to derive a preconditioner for self-consistent iteration similar to section IX in arXiv:0804.2583. For simplicity, consider here only one orbital (one or two electrons) systems. Suppose that ...
0
votes
0answers
28 views

how to estimate the number of people on a street within an hour?

I tried to use opencv to analyze a video filmed on the street. But the problem is the performance is not enough. I think the number of people must follow the poisson distribution. So I want to ...
0
votes
0answers
38 views

Techniques to optimise the integral of a function of known analytical form

I need to compute repeatedly a function that depends on an integral. The integral is not solvable analytically, but it depends on the argument of the function parametrically, like this: $$ f(x) = \...
0
votes
0answers
90 views

Super Computer cluster service for students [duplicate]

I would like to ask whether there is a free of charge supercomputer for students, in my current institution they don't have a supercomputer, I am running MD simulations, I would be grateful for ...
0
votes
0answers
86 views

Trouble with creating correct element matrices in Finite Element Analysis for a cantilever beam

I'm trying to solve for displacements of a cantilever beam numerically with FEA. The beam is modeled as a 3D-solid made up of a set of 8-noded hexahedral elements, which are in their undeformed state, ...
0
votes
0answers
42 views

Optimization (best input variables search) for a non-smooth non-linear unknown function

I am trying to optimize a system that monitors and advises a user multiple times over a certain period of time depending on changing outside factors. The systems behavior can be altered by 5 ...
0
votes
1answer
188 views

How to use LAPACK function (DGELSY) in Fortran

I am trying to use Least Squares Minimization to solve a the matrix problem: b = A*x for x. The system is overdetermined, and A is a dense matrix. In the LAPACK library, I believe the routine DGELSY ...
0
votes
0answers
41 views

Testing the SUPG method and other methods for hyperbolic equations

I am interesting in integrating the simple equation $$ \frac{\partial \phi}{\partial t} + \mathbf{u}\cdot\nabla \phi = 0 $$ with a Dirichlet boundary condition at the influx boundary ($\mathbf{u} \...
0
votes
0answers
28 views

Numerically computing deflection due to thermal expansion

Using linear elasticity formulation, I am attempting to numerically compute the displacement due to thermal expansion. This is done for a 3-D isotropic material. The governing equations are simply: ...
0
votes
0answers
14 views

Compute cell values and isosurface constant from density values of particles

I am trying to reconstruct the surface for a fluid simulation based on a list of particles using the Marching Cubes algorithm. From different resources, such as http://paulbourke.net/geometry/...
0
votes
0answers
30 views

Fit exponential convergence

I'm working with a numerical algorithm whose output $y$ asymptotically approaches a certain unknown value $a$. I expect an exponential convergence, i.e. the data $y$ given by my algorithm should be ...
0
votes
0answers
23 views

Simulating Anderson model, have problem with momentum representation (MATLAB)

I want to change from real-space representation to momentum-space representation I have a Hamilton-operator (Anderson-model), and I calculated some kind of entropy of its eigenstates (this is working, ...
0
votes
0answers
55 views

Computing Trajectory Equations of Kerr Geodesics

I want to numerically solve the trajectory equations of a Kerr geodesic given by wikipedia in Matlab. The trajectories look like: I implemented the equations and solved it with the standard Runge-...
0
votes
0answers
51 views

Stokes Equation fails to converge for an ellipse

This might be because of the mesh, but the following code blows up for all values of b not 1. Does anybody have any experience working with the ellipse mesh in Fenics? ...
0
votes
0answers
29 views

Finding the polynomial for the solution of an ODE

I’m stuck trying to solve part (b) and (c) of the below problem, but part (b) is the one of main concern here as I think (c) should follow easily once (b) is completed. I don’t know where to start ...
0
votes
0answers
52 views

Strange spectral (finite) element results of a solid plate

I tried to implement the spectral element method1 as proposed in [1]. I simulated an aluminum plate and a vertical concentrated force was applied at the middle of the top left quarter of the plate2. I ...
0
votes
0answers
44 views

Neumann boundary conditions on arbitrary surface for finite difference diffusion

I am facing the following problem, formulated in practical terms: I have a region $\Omega$ in two or three dimensions, represented as a binary mask, and an initial density $u_0$ within that region ...
0
votes
0answers
225 views

Solving a non-convex optimization problem using fmincon

I am trying to solve a non-convex optimization problem using fmincon(). At each iteration, I am iteratively looking for the optimum value and when the termination ...
0
votes
0answers
113 views

Double mach reflection at a inclined wedge

I am running into a strange problem when solving the 2D compressible Euler equations on a inclined wedge. To elaborate, my top boundary condition seems to emitting some type of instability. I have ...
0
votes
0answers
34 views

Extracting raw data from a graph

I am supposed to do black-box modeling, I need the input and output data for training, however, the authors didn't agree to share the raw data with me, I can only use their publication where these ...
0
votes
0answers
67 views

The final Boundary Condition is Unknown, Is Backward Euler is still valid to be implemented?

I am working on conductive polymer modeling and supposed to do one-dimensional diffusion model in the thickness of the polymer, however, due to the small value of thickness in micro, when I use the ...

15 30 50 per page
1
173 174
175
176 177
181