NDLI: Pseudospectral picture of the seismic algorithm for Surface Multiple Attenuation

Content Provider	IEEE Xplore Digital Library
Author	Drmac, Z. Grubisic, L. Jericevic, Z.
Copyright Year	2011
Description	Author affiliation: Department of Mathematics, University of Zagreb, Bijenička 30, 10000 Zagreb, Croatia (Drmac, Z.; Grubisic, L.) \|\| Department of Computer Engineering, Faculty of Engineering, Vukovarska 58, 51000 Rijeka, Croatia (Jericevic, Z.)
Abstract	Given a two dimensional signal, e.g. surface recorded 2-D wavefield representing marine seismic data, we consider the task of removing interference due to multiple reflection (the so-called multiple events). This problem has been successfully modeled in the previous work by the third author as a matrix optimization problem. Seen from an abstract perspective, this optimization problem is related to the mathematical notion of the pseudo-spectrum of a matrix. Pseudospectrum is sometimes also called the spectral portrait of a matrix and it presents a powerful visualization tool for analyzing the spectral properties of a nonhermitian matrix. Furthermore, there are several freely available visualization tools to plot the pseudospectrum, most notably Tom Wright's eigtool. The original mathematical problem can be seen as a quest for a minimum of the pseudo-spectrum for a certain nonhermitian matrix pair in the Frobenius matrix norm. The main contribution in the original algorithm was the use of matrix Eigenvalue Decomposition in the preprocessing stage of the algorithm to reduce the computational cost of the inner optimization loop to $O(n^{2}).$ Here n is the dimension of the matrix (e.g. 2-D signal). Motivated by the work of Trefethen and Embree we present several alternatives to this algorithm which also result with inner loop procedures having $O(n^{2})$ cost. Such algorithms are based on a Lanczos type singular value algorithms and solve the optimization problem with regard to the spectral norm. We further discuss issues related to the numerical stability of the complete algorithm with regard to the choice of the used matrix factorization in the preprocessing module. Notably, we consider the cost and the numerical stability for the use of Schur, Hessenberg and Eigenvalue matrix decompositions in this context. Let us close by noting that the structure of the computational problem is such that it falls into the class of “embarrassingly parallel” algorithms. This will also be exploited in the presented work.
Starting Page	320
Ending Page	325
File Size	6079218
Page Count	6
File Format	PDF
ISBN	9781457709968
e-ISBN	9789532330595
Language	English
Publisher	Institute of Electrical and Electronics Engineers, Inc. (IEEE)
Publisher Date	2011-05-23
Publisher Place	Croatia
Access Restriction	Subscribed
Rights Holder	MIPRO
Subject Keyword	Context Algorithm design and analysis Symmetric matrices Approximation algorithms Eigenvalues and eigenfunctions Matrix decomposition Optimization
Content Type	Text
Resource Type	Article

Sl.	Authority	Responsibilities	Communication Details
1	Ministry of Education (GoI), Department of Higher Education	Sanctioning Authority	https://www.education.gov.in/ict-initiatives
2	Indian Institute of Technology Kharagpur	Host Institute of the Project: The host institute of the project is responsible for providing infrastructure support and hosting the project	https://www.iitkgp.ac.in
3	National Digital Library of India Office, Indian Institute of Technology Kharagpur	The administrative and infrastructural headquarters of the project	Dr. B. Sutradhar bsutra@ndl.gov.in
4	Project PI / Joint PI	Principal Investigator and Joint Principal Investigators of the project	Dr. B. Sutradhar bsutra@ndl.gov.in Prof. Saswat Chakrabarti will be added soon
5	Website/Portal (Helpdesk)	Queries regarding NDLI and its services	support@ndl.gov.in
6	Contents and Copyright Issues	Queries related to content curation and copyright issues	content@ndl.gov.in
7	National Digital Library of India Club (NDLI Club)	Queries related to NDLI Club formation, support, user awareness program, seminar/symposium, collaboration, social media, promotion, and outreach	clubsupport@ndl.gov.in
8	Digital Preservation Centre (DPC)	Assistance with digitizing and archiving copyright-free printed books	dpc@ndl.gov.in
9	IDR Setup or Support	Queries related to establishment and support of Institutional Digital Repository (IDR) and IDR workshops	idr@ndl.gov.in

An Implementation of Matrix Eigenvalue Decomposition with Improved Jacobi Algorithm

Bit-level systolic algorithm for the symmetric eigenvalue problem

An approximate polynomial matrix eigenvalue decomposition algorithm for para-Hermitian matrices

New fast optimal window design algorithm based on the eigen-decomposition of the symmetric Toeplitz matrix

Parallel QR algorithm for the complete eigensystem of symmetric matrices

Incremental Nyström Low-Rank Decomposition for Dynamic Learning

Causality-constrained multiple shift sequential matrix diagonalisation for parahermitian matrices

An algorithm for computing singular values of large matrices for use in the analysis of large systems

Closed-form optimization of covariance intersection for low-dimensional matrices

Pseudospectral picture of the seismic algorithm for Surface Multiple Attenuation

Similar Documents

An Implementation of Matrix Eigenvalue Decomposition with Improved Jacobi Algorithm

Bit-level systolic algorithm for the symmetric eigenvalue problem

An approximate polynomial matrix eigenvalue decomposition algorithm for para-Hermitian matrices

New fast optimal window design algorithm based on the eigen-decomposition of the symmetric Toeplitz matrix

Parallel QR algorithm for the complete eigensystem of symmetric matrices

Incremental Nyström Low-Rank Decomposition for Dynamic Learning

Causality-constrained multiple shift sequential matrix diagonalisation for parahermitian matrices

An algorithm for computing singular values of large matrices for use in the analysis of large systems

Closed-form optimization of covariance intersection for low-dimensional matrices

Pseudospectral picture of the seismic algorithm for Surface Multiple Attenuation