Questions tagged [projection]

The tag has no usage guidance.

Filter by
Sorted by
Tagged with
2
votes
1answer
54 views

How to robustly and numerically expand a $k$-order polynomial in two variables defined on a polygon domain?

Given a $k$-order polynomial in two variable $p(x, y)$ defined on a polygon domain $K$. And I want to numerically expand it to the following form $$ p(x, y) = c_0 + c_1 x + c_2 y + c_3 x^2 + c_4 xy + ...
0
votes
0answers
22 views

Reduce projection error while retaining similar amount of elements in CG-FEM

Based on the answers I got to my questions (Interpolation of function onto mesh gives different results, depending on mesh density and Solving a non-linear heat equation with the galerkin method gives ...
2
votes
2answers
69 views

projective reconstruction from orthogonal views

This is a problem from projective geometry. Suppose I have a vector $z \in R^k$ of unit length $\| z \| =1$ inside a $k$-dimensional hypercube. I don't know its value but do know its projection upto ...
2
votes
1answer
107 views

Projection onto the set of Orthogonal matrices

Let $M \in \mathbb{R}^{n \times n}$ and denote the set of Orthogonal matrices by \begin{equation} \mathcal{O}_{n} = \left\lbrace Q \in \mathbb{R}^{n \times n} \colon QQ^{T} = \mathbb{I}_{n} \right\...
2
votes
0answers
56 views

Geometric interpretation of lemma

I am currently studying eigenvalue problems. I already worked through the Minimax-principles, seen why $\lambda_{h, m} \geq \lambda_m$, when comparing the eigenvalues of a discretization and the ...
5
votes
0answers
65 views

Padua-type pointset for functions singular on line $x=y$

The Padua points $\mathrm{Pad}_{n} \subset [-1,1]^{2}$ are a unisolvent pointset with optimal growth of Lebesgue constant, described in detail here. With some work they can be used to generate a ...
0
votes
1answer
61 views

Wanting an explanation of the variables in Iterative PCA algorithm

I've been trying to implement the CPU GS-PCA algorithm in this paper . The code starts on page 28 I have a program written a script in python which gives the same output as the C++ program in the ...
1
vote
1answer
111 views

Questions about iterative projection methods in Saad book

I am reading Chapter 5 of Saad's iterative methods book, and I don't understand section 5.2.1 about the two propositions of optimality results. In the statements of the propositions, what does it mean ...
1
vote
0answers
61 views

Stable Method of orthogonal projection onto a subspace with the help of Moore-Penrose inverse,

Projection of a vector $v$ onto the column space of a matrix $A$ is given by $AA^\dagger v$. From the definition of Moore-Penrose Inverse we know that $AA^\dagger v = (A^T)^\dagger A^T v $. Below is ...
3
votes
0answers
73 views

Backward stable algorithm to get orthogonal projection onto the column space of a matrix

I have to find the orthogonal projection of a vector $b$ onto the matrix $A$ of size $m \times n$. In my application, I don't have the luxury of calculating the QR factorization. All I have are ...
4
votes
3answers
325 views

How can I make sure the flow is divergence-free when I use moving mesh?

I am using projection method and P2/P1 finite element method to solve the incompressible Navier-Stokes equations while the mesh is constantly adapted as the body moves (edge swapping, splitting and ...
5
votes
1answer
87 views

$L_2$ projection with integer constraints and prescribed sum

Suppose I am given a vector $v^0\in\mathbb{R}^n$ and integers $k,\ell\in\mathbb{Z}$. Assuming $k$ is close to zero (e.g. $0\leq k\leq5$), is there an algorithm for solving the following integer ...
1
vote
0answers
76 views

How does this Constrained Minimization algorithm work?

I don't fully understand the subsection 3.2 Constrained minimization of this paper. In particular, I don't understand the first step "Register active set" and the definition of projection $P(x)$. ...
5
votes
2answers
261 views

Efficiently removing projection to subspace without having an orthogonal basis

I have a number of vectors $v_1, …, v_n$ and another vector $w$, that all are linearly independent, but not orthogonal. Let $V := \mathop{\text{span}}(v_1, …, v_n)$. I need to remove $w$’s projection ...
2
votes
1answer
203 views

Minimizing linear objective on intersection of convex sets

Suppose I wish to solve the following optimization problem: $$ \begin{array}{rl} \min_{\mu\in\mathbb{R}^n} &\mu^\top c\\ \textrm{subject to} & \mu\in C_1\cap C_2\cap\cdots\cap C_k, \end{array} ...
0
votes
1answer
78 views

Regarding solution vector of the wave equation

I am simulating the wave equation using FEM. For a 2D wave equation, when I visualise my output in Paraview, I see a separate solution in 'x' and 'y' direction for each node on the mesh. Therefore, if ...
1
vote
0answers
171 views

How to obtain projections from sinogram in ART reconstruction technique?

I'm kind new in the Computed Tomography field and I'm trying to understand and implement ART technique. Said it, I started to read the book The Mathematics of Medical Imaging - A Beginners Guide by ...
1
vote
0answers
81 views

Oscillations in Chorin's method due to the BC

I am pretty new to the CFD and I wanted to start with Chorin's projection. The starting problem is just a free jet flowing in the investigated area. I got terrible oscillations almost immediately and ...
3
votes
2answers
171 views

Implicit projection method with inflow boundary conditions

I am trying to use a projection method that deals with the viscous effects implicitly to model flow around a cylinder. I'm having trouble figuring out what the boundary conditions should be, ...
1
vote
0answers
35 views

BC's for intermediate velocities in Implicit Fractional Step Methods

Lately, I was reading some seminal papers on Fractional Step Algorithms and I found this one: Kim, D., Choi, H. A Second-Order Time-Accurate Finite Volume Method for Unsteady Incompressible Flow on ...
2
votes
1answer
96 views

Projecting a vector field onto a H(div) space

I've got a uniform quadrangular mesh and for each node there's a vector quantity $u$ defined. I also have a non-aligned material interface across the mesh. Now I need that vector quantity to have a ...
1
vote
1answer
42 views

How to compute $\mathrm{proj}_{SDP}(C\odot X)./C$ without numerical problems?

I have a matrix, $X$, it is symmetric. I project $C \odot X$ and $D\odot X$ to semidefinite cone. $C$ is a Gramian matrix with some elements near zero and of course semidefinite, with one row and ...
7
votes
1answer
115 views

Approximation properties of FEM projections operators on a boundary

We have an elliptic projection $$P: V \rightarrow V_{h}$$ which satisfies $$\Vert u - Pu \Vert_{L^{2}(\Omega_{e})} \leq Ch^{k+1} \enspace .$$ Can we say anything about $\Vert u - Pu \Vert_{L^{2}( ...
2
votes
0answers
149 views

How to obtain the reduced model from a subspace projection method?

I have a system of ordinary differential equations (ODEs). It is a large system that has dozens of equations and hundreds of parameters. I wish to reduce its size so it becomes computationally more ...
2
votes
1answer
133 views

How to project a 0 genus mesh model on a sphere?

I have a mesh which represents a 0-genus model. My goal is to construct a homeomorphism from that mesh to its bounding sphere. I'm trying to understand a paragraph in http://citeseerx.ist.psu.edu/...
5
votes
2answers
655 views

Incompressible Navier-Stokes equations: Is projection method exact?

Is the projection method of integrating Navier-Stokes equations exact? Take the incompressible flow equations: $$ \frac{\partial\mathbf{u}}{\partial t} = -\mathbf{u}\cdot \nabla \mathbf{u} -\nabla ...
6
votes
1answer
496 views

How to project a vector into the H(div) space (in the context of finite elements)?

Say I have a simple elliptic PDE: $$ -\nabla\cdot(K\nabla p) = f \;\;\;\text{in}\;\Omega $$ with the appropriate boundary conditions. I solve for $p$ using a FEM (a discontinuous Galerkin method to ...
4
votes
1answer
277 views

Active Elements in Projected Newton's Method?

To those who are familiar with the projected Newton's method or projected gradient method... We consider a constrained optimization problem with simple bounds. Particularly, minimize f(x) subject to ...
8
votes
5answers
900 views

Motivation behind Galerkin method

I have a question about Galerkin method. I don't understand why the Galerkin method weights the residual by the shape functions and sets it equal to zero. I want to know what is reason of this. Why we ...
2
votes
0answers
233 views

Treatment of Neumann (Traction) boundary conditions using projection methods

I am looking to solve the incompressible Navier-Stokes equations in 3D, using an inflow boundary condition specifying a velocity: $\mathbf{u} = \mathbf{g}_0 \,\, \forall \,\, \mathbf{x} \in \Gamma_u$ ...
4
votes
2answers
1k views

Projection Method: Boundary condition on intermediate velocity field

I'm trying to solve variable density and viscosity Navier-Stokes equation using lagged pressure projection method. I'm solving for cavity problem as a test case now (once I get projection right, I ...
2
votes
2answers
173 views

Projection of vector field on to a gradient field

Say I have a vector field with non-zero curl, therefore the potential function depends on the path I choose to integrate. In this paper the authors proposed to project the vector field into a gradient ...
1
vote
0answers
60 views

Mass-conservative reprojection (on a sphere)

I have a 2D distribution of mass on a sphere given as a matrix of masses in latitude-longitude grid cells. I need these masses projected to another grid on the same sphere with different location of ...
1
vote
1answer
1k views

Lid driven cavity flow problem to solve with the help of projection method and use coding in fortran 90 [closed]

I want to solve Lid driven cavity flow with the help of projection method i want to make own fortran program to solve two dimensional unsteady flow. I propose to solve the equations in the following ...
4
votes
0answers
670 views

Generating pseudo-random orthonormal bases for random projection

I am performing series of random projections i.e. projecting the input matrix onto randomly generated orthonormal bases (of much lower dimensionality). The projection is just a matrix multiplication ...
3
votes
2answers
875 views

How to properly use polynomial projection to get values at visualization nodes?

I am trying to implementing a nodal discontinuous Galerkin spectral element method for linear and non-linear systems of equations. The solution at each time step is given at ...