All Questions

Filter by
Sorted by
Tagged with
0
votes
1answer
66 views

How can PDE matrices be identified?

I need to include experimental results for lots of PDE (partial differential equation) matrices in my research work. How can I identify PDE matrices? For example, matrices in the UFL Sparse Matrix ...
0
votes
1answer
299 views

Is there general algorithms to solve such 3D cutting problems?

Suppose a cuboid $\mathbb{A}$ has $L$,$M$ and $N$ as its length, width and height respectively, where $L\ge{M}\ge{N}>0$; Now we want to cut $\mathbb{A}$ into smaller cuboids with length $x$, width $...
0
votes
1answer
162 views

fast gradient method for convex piecewise linear function

What is the state of art gradient based algorithms in convex optimization solving non-smooth piece-wise linear functions? Thank you. EDIT: It is different from one of my previous post in the sense ...
0
votes
1answer
80 views

Mathematical error when attempting to represent step function using fourier series

I am attempting to work through a very simple problem. Determine the Fourier series expansion for the following heat PDE problem with ICS and BCS: $$ u_{t} = \alpha^2u_{xx}$$ $$ u(0, t) = u(L, t) =...
0
votes
1answer
2k views

Gonzalez algorithm

I'm getting confused about implementing the Gonzalez algorithm for k-center clustering. The algorithm is : let $S$ be a data set Choose $z_1 \in S$ arbitrarily. Number $z$ is center. (In general ...
0
votes
1answer
212 views

Strict Feasibility in Interior Point Methods

As we know, in the interior point methods, all the iterates have to be strictly feasible. I implemented an affine scaling interior point for nonlinear objective functions. For small examples (2D), it ...
0
votes
1answer
624 views

Interpolating 3D Array non-monotonic data in MatLab

I am working on creating a program for simulations where three variables are parametrized, and we modify one parameter while keeping the other two constant. An example array looks like this when ...
0
votes
1answer
1k views

Temperature dependent 1-d conduction in Python?

I'm trying to write a Python code that is a numerical solver for 1-d heat conduction (using FVM) with a temperature dependent thermal conductivity. The solver has three functions I need to iterate ...
0
votes
3answers
1k views

Solving a double integral in Matlab

I need to solve this integral: $M = \rho h \int_0^a \int_0^b N^2 \ dx dy$ where $N$ is defined as: $N = sin\left(\dfrac{\pi}{4}\eta +\dfrac{3\pi}{4}\right)sin\left(\dfrac{\pi}{4}\eta +\dfrac{3\pi}{...
0
votes
2answers
394 views

Faster methods for projecting a mesh onto a hierachally unrelated mesh?

I have a set of independent meshes whose results I would like to project onto another non-hierachally related mesh. Until now, I've been accomplishing this by finding the nearest-distance node in the ...
0
votes
1answer
694 views

How to solve Energy Balance equation by numerical method

Good Day I am new to heat transfer technique please give me some suggestion on solving energy balance equation $$a \frac{\partial T_p}{\partial t}=\frac{\partial}{\partial x}\left(b\frac{\partial ...
0
votes
1answer
77 views

How to present a polynomial interpolation if the first order information of function f(x) is given?

Suppose $f'(x_1),\ f'(x_2),\ f'(x_3)$ are given, how to give a polynomial interpolation $p(x)$ such that $p'(x_1)= f'(x_1),\ p'(x_2)=f'(x_2),\ p'(x_3)=f'(x_3)$? And how to give an error analysis?
0
votes
2answers
157 views

Computing eigenpairs of singular matrix with ZGEEV?

I've never run into a singular matrix before, so bear with me. I have a complex non-symmetric matrix (about 1000 x 1000) that I know has a couple zero eigenvalues. It isn't guaranteed to be ...
0
votes
1answer
980 views

What is Mesh Independence Report?

I am performing analysis on chassis (Static Structural) and for optimization purpose i am asked to generate MESH-INDEPENDENCE REPORT,of which i have no idea. I have tried going through research papers ...
0
votes
1answer
70 views

Unexpected behavior, values tend to converge instead of fluctuate. (MD)

I am writing a molecular dynamics program to create an lattice and populate it with atoms/molecules. These then are given random velocities and the system is initialized. Then throughout time the ...
0
votes
1answer
84 views

Largest Cylinder inside Polyhedron

Imagine you have a piece of wood and from that piece you want to get the largest cylinder possible. How do you determine the position and orientation of the cylinders axis, to maximize its radius? ...
0
votes
1answer
116 views

CVX : Obtaining the minimizing parameter at the optimum

In CVX, how do we return the value of the parameter over which the problem is minimized at the optimal value? By this, I mean, how do we obtain $$x^* = \arg\min_x f(x)$$ when solving the problem ...
0
votes
1answer
106 views

Actually calculating the rate of convergence of iteratvie methods when exact solution is unknown

When solving a system of nonlinear equations using iterative methods, the rate of convergence usually is defined by the following formula: (1) where x* is the exact solution. However usually we ...
0
votes
2answers
1k views

Simulating a traveling sine wave

I'm trying to make an animation of a travelling sine wave (amplitude vs. position) would anyon here happen to know how to do so?
0
votes
1answer
155 views

MAX-SAT and MAX-cut

I have been using MAX-SAT solver to obtain the exact ground state of ising spin glass model: For 1D periodic model, for systems with 50 binary variables and interaction range of 15th nearest ...
0
votes
1answer
255 views

Why doesn't the integral of a smoothing kernel equal 1?

I was checking the unity condition on a smoothing kernel for SPH, however I don't understand why the integral is not giving 1. The kernel function is $\displaystyle W(||\vec r||,h) = \frac{315}{64\pi ...
0
votes
2answers
385 views

How to represent molecules and compare equality

I originally asked this question at StackOverflow, and was suggested to bring it here. I've seen this question about the representation of molecules in memory, and it makes sense to me (tl;dr ...
0
votes
2answers
242 views

Statistical analysis of optimization algorithms

If we optimize some parameter using 4 optimization algorithms, 2 of which are population based (say A and B) and 2 trajectory methods (single point search)(say C and D); what statistical test can be ...
0
votes
1answer
38 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: ...
0
votes
2answers
112 views

How to determine the number of c points in algebraic multi grid

I am trying to write an algebraic multi-grid solver (in c++). At a given level I determine which nodes are c-points and which nodes are f-points (where the total number of c and f points equals the ...
0
votes
1answer
199 views

jump conditions for Poisson/Darcy equation in primal form versus mixed form

Consider the Darcy equation, $$\mathbf{v} + \dfrac{k}{\mu_0}\nabla p = \mathbf{f} \\ \mathrm{div}\; \mathbf{v} = 0$$ If the coefficient $k$ is piecewise constant across an interface $\Gamma$ in the ...
0
votes
1answer
475 views

Python: Computation time Issue with mpi4py

I am using in Python mpi4py to process in parallel 20 minimization functions. Each of the 20 worker processes the same algorithm but with different random initial ...
0
votes
1answer
138 views

Writing a 3D array from Petsc

I am trying to do something fairly simple somehow I have made it hard. Is there an example of how to send a 3D array to a binary file?
0
votes
1answer
65 views

FEM with soil slope

I what to calculate displacements, stress and strain in a soil slope with a FEM script. The slope moves like a laminar flow. Can you suggest me some bibliography on this problem? I've already look on ...
0
votes
1answer
59 views

Should I use bestman2 or dcache CLIs (srmping, srmls, etc..)

I am fairly new to using SRM, and I am wondering if I should be using srm CLI tools (srmping, srmls, srmcp, etc..) provided by dcache, or bestman2 counter parts (srm-ping, srm-ls). I used to think ...
0
votes
1answer
73 views

How many years it would take our laptops to be as fast as the fastest super computer in 2000 [closed]

How many years it would take our laptops to be as fast as the fastest super computer in 2000 according to Moore's law?
0
votes
1answer
70 views

Finding roots without knowing much about the function

Consider solving numerically for roots: $( x_0, y_0): f(x_0, y_0) = 0, g(x_0, y_0) = 0$ where you only know that f, g continuously differentiable but the theoretical differentiation is not a ...
0
votes
1answer
69 views

Particle Collision to Static paticles

I have a system of particles with equal distance with each other and another at random positions which is moving with time. I want to know: a) The method by which I can reduce the number of particles ...
0
votes
1answer
536 views

Solver for large non-linear system of equations

I am curently using R package nleqslv for solving a non-linear system of equations with 300 variables. I need to scale this to the system with ~50k variables and naturally this does not scale very ...
0
votes
1answer
472 views

MATLAB Newton non-linear equation

I have the following non-linear equation: where $w0=0.25,w0=0.5,w0=0.75$. I have to prove that if $k$ is a root, then also $−k$ is a root and that there exists only one $k∈(0,1)$ root, but my MATLAB ...
0
votes
1answer
420 views

Elastic LP Programming

Say I have an LP that is unfeasible and that I want to find the solution that makes it feasible without strongly violating the current constraints. What is a principled way of solving this problem, ...
0
votes
1answer
80 views

Rank of image intensity matrix

I've been reading a paper about using Matrix Completion for Photometric Stereo but I am having some troubles in section 2.2 trying to understand why irrespective of the number of pixels and the number ...
0
votes
2answers
283 views

Systems of nonlinear equations

Consider the nonlinear system of equations $$ (1) \quad qk^2a_1^2E^2+wna_0a_1AE+pnka_0^2a_1E+rn^2a_0^2A^2-rn^2a_0^3A^2+qk^2a_0a_1ABE-qk^2na_0^2E^2=0, $$ $$ (2) \quad nwa_1^2AE+nwa_0a_1AB+2pnka_0a_1^...
0
votes
2answers
289 views

WENO reconstruction of flux involving derivative terms

I have a set of modified compressible Euler equations that I would like to solve using a WENO method. The issue is that the modified flux function involves derivative and filtering terms and I'm not ...
0
votes
1answer
505 views

Truncated SVD implementation in Java

I need the Truncated SVD implementation in java. I need to pass a matrix of doubles and an integer value representing the rank where to filter out noise. In output i need a filtered matrix of doubles. ...
0
votes
1answer
2k views

How to detect smiles using OpenCV and IntraFace?

I have installed OpenCV and IntraFace, now the code can detect landmarks on the face, like this : Now how can I convert this to a smile detector? I have the landmark points as an array, so I should ...
0
votes
1answer
400 views

Custom CVX Functions - Overriding Errors [closed]

I have constructed a function in Matlab that is convex and increasing (qualitatively similar to an exponential function but I am hoping to avoid the successive approximation requirement of exp). In my ...
0
votes
1answer
77 views

Which version of FFTW

At the moment, I'm installing the "GADGET-2" application. In its documentation it says "Note that the MPI-capable version 2.x of FFTW is required, the new version 3 lacks MPI capability at this point. ...
0
votes
1answer
52 views

Best way of getting updates on intermediate simulation results

I will be interfacing our fluid dynamics code with an optimizer. I expect the optimization process to take 3-4 days for any decent decrease. But I will want to have a look at the intermediate results ...
0
votes
1answer
68 views

SVM Math Question

I'm studying support vector machines and came across this paper. The following equation doesn't make sense to me, especially the part with the 0 ∀i. Any help understanding the basics of SVMs? yi * (...
0
votes
1answer
305 views

How could I implement this neural network in R? [closed]

In Establishing, versus maintaining, brain function: a neurocomputational model of cortical reorganisation after injury to the immature brain, Varier et al develop a neural network model of motor ...
0
votes
1answer
761 views

Degradation Model in MATLAB

I am trying (using MATLAB) to generate the following image (from the section 4.1 Numerical Examples) of Wu Tian Chen research article 'Condition-based Maintenance Optimization Using Neural Network-...
0
votes
1answer
66 views

Can I use Engineering codes like ASHRAE , NFPA and UPC in my software without permission? [closed]

Say, I want to use the UPC code Water Supply Fixture Unit tables in my software that I'm writing ... is it okay or there will be copyright issues?
0
votes
1answer
112 views

max speed <--> time discretization

I'm working on a heat diffusion problem, $$ \frac{\partial T}{\partial t}=\vec{\nabla}\cdot\left(\kappa T^{5/2}\,\vec{\nabla}T\right) $$ I know that, after discretization, the time step for the 1D ...
0
votes
1answer
235 views

Implementation of the MARTINI force field

Given an input parameters for the MARTINI force field such as: ...

15 30 50 per page
1
160 161
162
163 164
181