Questions tagged [projection]

For questions about implementing mappings from a set of elements (points, vectors, etc) to some subset of these elements.

Filter by
Sorted by
Tagged with
1 vote
2 answers
88 views

Efficient projection onto the kernel of a matrix

Suppose I have a positive semidefinite matrix $M = \sum_i^N A_i^T A_i$ where each $A_i$ is a fat matrix of shape (m,n) and $m << n$, we can also assume that $A_i$ is full rank(but the stacked ...
HRI's user avatar
  • 113
1 vote
0 answers
35 views

How to find the formula of a projected circle in a pencil of conics structure?

Hi this is my first question on the platform so feel free to comment if I have a mistake regarding the question. I'm working on an ellipse detection scheme in which I have markers consisted of 3 ...
kemal alperen cetiner's user avatar
1 vote
0 answers
122 views

From 3D to 2D with a STL file

I would like to do a 2D projection from a 3D geometry saved in a stl file and know the distance between the two projected planes. In order to explain better the concept I will start with an almost ...
Daniel's user avatar
  • 99
0 votes
1 answer
109 views

How is the integral of a projection over an element $T$ computed in practice? (deal.II related)

I'm studying an error estimator for the equation $\nabla\cdot(\beta u) + cu = f$ and it contains the following term $$||f - cU_h - \Pi(f-c U_h) ||_T$$ where : $\Pi$ is the local orthogonal $L^2$ ...
FEGirl's user avatar
  • 281
0 votes
0 answers
168 views

Do the projection (with Jacobian) and marginalisation (inversion of matrix and remove a row/column and reinversion) commute?

I try to check the equality or the inequality between 2 Fisher matrices. The goal is too see if the projection (with Jacobian) and marginalisation (inversion of matrix and remove a row/column and ...
user avatar
1 vote
0 answers
60 views

Why does including the pressure in this FVM for Stokes 2nd Problem lead to wrong solutions?

I'm trying to learn how to use finite volume methods and I want to solve a more general case of Stokes' second problem i.e. an infinite half-plane oscillating harmonically with no-slip boundary ...
Wihtedeka's user avatar
  • 136
6 votes
1 answer
248 views

Projection method FVM poisson part, adding source term

The idea of the method is to decompose the Navier-Stokes equation into the solenoidal and irrotational parts. $$\frac{\partial u}{\partial t}+u(\nabla \cdot u)=-\frac{1}{\rho}\nabla p+\nabla ^2 u$$ ...
2Napasa's user avatar
  • 362
0 votes
1 answer
57 views

Projecting nodal solution to gauss points with certain accuracy

I am having a problem that was also mentioned at the accepted answer to this question by Wolfgang Bangerth. I need to calculate, as it was specified at the question at the link, F1 integral and for ...
noname 's user avatar
1 vote
0 answers
226 views

Lumped mass matrices for higher-order finite elements for CFD

Given that some of the mass lumping techniques, for example, row-sum lumping does not produce practically viable lumped mass matrices for all the element shapes, what are the techniques used for mass ...
Chenna K's user avatar
  • 914
2 votes
1 answer
225 views

Efficiently compute a projection matrix from Householders reflectors

Let $A \in \mathbb{R}^{m \times n}$ where $m \geq n$. Let $B$ and $\tau$ be the result of applying LAPACK's dgeqrfp method (R on the upper right triangle, and the ...
Matthias Beaupère's user avatar
3 votes
1 answer
199 views

Project to nearest point on convex polytope

I have a point $y \in \mathbb{R}^d$ and a convex polytope $\mathcal{P}$ given as the intersection of half-spaces: $$\mathcal{P} = \{x \in \mathbb{R}^d \mid a_1 \cdot x \le b_1, \dots, a_n \cdot x \le ...
D.W.'s user avatar
  • 400
1 vote
0 answers
100 views

Stokes problem with imposed acceleration on boundaries (projection scheme)

I am trying to solve FSI problems with finite elements and using a projection scheme (I am taking as reference the review of Guermond: Guermond, J. L.; Minev, P.; Shen, Jie, An overview of ...
gc11's user avatar
  • 11
2 votes
1 answer
86 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 + ...
Huayi Wei's user avatar
  • 143
2 votes
2 answers
101 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 ...
drw's user avatar
  • 203
2 votes
1 answer
550 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\...
JKay's user avatar
  • 121
2 votes
0 answers
78 views

Geometric interpretation of lemma

I am currently studying eigenvalue problems. I already worked through the minimax-principles and seen why $\lambda_{h, m} \geq \lambda_m$, when comparing the eigenvalues of a discretization and the ...
mdot's user avatar
  • 155
5 votes
0 answers
80 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 ...
user14717's user avatar
  • 2,155
0 votes
1 answer
73 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 ...
sav's user avatar
  • 257
1 vote
1 answer
135 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 ...
user123's user avatar
  • 679
1 vote
0 answers
87 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 ...
Abhay Gupta's user avatar
3 votes
0 answers
143 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 ...
Abhay Gupta's user avatar
4 votes
3 answers
927 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 ...
shidi.yan1992's user avatar
5 votes
1 answer
119 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 ...
Justin Solomon's user avatar
1 vote
0 answers
100 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)$. ...
sarasvati's user avatar
  • 185
5 votes
2 answers
744 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 ...
Wrzlprmft's user avatar
  • 2,032
2 votes
1 answer
270 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} ...
Justin Solomon's user avatar
0 votes
1 answer
102 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 ...
CRG's user avatar
  • 347
1 vote
0 answers
345 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 ...
Flávio Schuindt's user avatar
1 vote
0 answers
93 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 ...
Victor Pira's user avatar
3 votes
2 answers
266 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, ...
Lukas Bystricky's user avatar
1 vote
0 answers
38 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 ...
Kbzon's user avatar
  • 251
2 votes
1 answer
147 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 ...
balborian's user avatar
  • 601
1 vote
1 answer
46 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 ...
user85361's user avatar
  • 221
7 votes
1 answer
159 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}( ...
Mike Harmon's user avatar
2 votes
0 answers
165 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 ...
milez's user avatar
  • 121
2 votes
1 answer
181 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/...
Brainless's user avatar
  • 203
5 votes
2 answers
1k 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 ...
Abhilash Reddy M's user avatar
6 votes
1 answer
789 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 ...
Justin Dong's user avatar
3 votes
1 answer
398 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 ...
Linh Huynh's user avatar
11 votes
5 answers
3k 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 ...
mohammad's user avatar
  • 111
2 votes
0 answers
347 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$ ...
computanjohn's user avatar
4 votes
2 answers
2k 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 ...
Pranav's user avatar
  • 309
2 votes
2 answers
407 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 ...
BRabbit27's user avatar
  • 1,029
1 vote
0 answers
64 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 ...
Roux's user avatar
  • 111
1 vote
1 answer
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 ...
sajjad ur's user avatar
5 votes
0 answers
708 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 ...
Adam's user avatar
  • 151
3 votes
2 answers
1k 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 ...
Michael Schlottke-Lakemper's user avatar