SlideShare a Scribd company logo
International Journal of Engineering Science Invention
ISSN (Online): 2319 – 6734, ISSN (Print): 2319 – 6726
www.ijesi.org ||Volume 6 Issue 1|| January 2017 || PP. 38-39
www.ijesi.org 38 | Page
“ Implimentation of SD Processor Based On CRDC Algorithm ”
Raj Kumar Verma#1
, Mrs.Akanksha Awasthi#2
#1
Research Scholar, C.V.R.U, Bilaspurc.g. India.
#2
Asst.Prof.C.V.R.U., Bilaspurc.g. India.
Abstract: In Digital Signal Processing (DSP) there are many complex algorithms for which an efficient
hardware implementation is required in real time applications. One such complex algorithm is Singular-value
Decomposition (SD) which is an important algorithm with applications in varied domains of signal processing
such as direction estimation, spectrum analysis and systems identification. It is a generalized extension to the
eigen-decomposition for non-square matrices and is hence of great importance, particularly for subspace based
algorithms in signal processing. But SD is known to be a very complicated algorithm with computational
complexity ~O(N3
) (for a NxN square matrix). For real-time computation of such a complex algorithm the use of
a parallel and direct mapped hardware solution is indeed desired.
Hardware and software resources:
1. OS: Windows 9x or upper
2. RAM : Minimum 512 MB
3. Programming Language: XILINX 8.6 or upper ( VERILOG)
I. Introduction
In this chapter basic idea of the project undertaken is explained. Starting from the aim of the project,
Importance and the problem associated with the design is also discussed. Finally line of approach to solve the
problem is presented and organization of rest of the report is briefly shown.
Importance
The Singular Value Decomposition (SD) is an important matrix factorization with applications in signal
processing, image processing and robotics. It is generally acknowledged that the SD is the only generally
reliable method for determining the rank of a matrix numerically. The SD is a very useful tool, for example, in
analyzing data matrices from sensor arrays for adaptive beam forming , and low rank approximations to
matrices in image enhancement . The wide variety of applications for the SD coupled with its computational
complexity justifies dedicating hardware to this computation. Designed 2X2 CRDC based SD processor can be
used as basic building block for an array of processors of N X N matrix.
Problem statement
The importance of hardware implementation is now known but there are certain complexities that come on the
way of implementation. In a typical computer algorithm for the SD, the sines and cosines of the rotation angles
are computed through formulas that require division and square root operations. Also the explicit angles are not
required and only thesines and cosines are computed. The rotations are then applied to the 2x 2 sub-matrix using
standard matrix multiplication techniques. But in this method many time-consuming operations such as
multiplication, division, and square rootare needed.
The SD is more closely mapped ontohardware through the useof the CRDC algorithms. CRDC algorithm
provides an iterative scheme, consisting of simple addition and binary shift operations, to compute trigonometric
values to any desired precision. The CRDC algorithms can provide the calculation of vector rotation and inverse
tangent without costly hardware that proves very beneficial for implementation of SD processor.
Line of approach
The complexity of SD computation and its real-time applications were the inspiration behind efficient hardware
implementation. In SD processor time consuming operations like multiplication, square root and division are
eliminated by the use of CRDC algorithm. CRDC algorithm only requires simple shift and additions which
efficiently maps SD processor.
“ Implimentation of SD Processor Based On CRDC Algorithm ”
www.ijesi.org 39 | Page
Basic CRDC architecture
An iterative CRDC architecture can be obtained simply by duplicating each of the three difference equations in
hardware as shown in Figure.
Figure: Basic CRDC hardware
The decision function, i
d is driven by the sign of the y or z register depending on whether it is operated in
rotation ( i
z ) or vectoring mode ( y i
). In operation, the initial values are loaded via multiplexers into the x, y
and z registers. Then on each of the next n clock cycles, the values from the registers are passed through the
shifters and adders-subtractors and the results placed back in the registers. The shifters are modified on each
iteration to cause the desired shift for the iteration. Likewise, the ROM (Look-up table) address is incremented
on each iteration so that the appropriate elementary angle value is presented to the z adder-subtractor.
On the last iteration, the results are read directly from the adder-subtractors. A simple state machine is required
to keep track of the current iteration, and to select the degree of shift and ROM address for each iteration.As
shown above, CRDC iteration requires 2 shifters, 1 table lookup and 3 adders. For n bits of precision, n
iterations are needed.
II. Conclusion
In linear algebra, the singular value decomposition (SD) is an important factorization of a rectangular real or
complex matrix, with several applications in signal processing and statistics.
A 2X2 CRDC based SD processor was designed and successfully simulated. CRDC algorithm was used to
simplify the computational complexities of SD calculation. Apart from calculating the SD of a 2x2 matrix the
processor is able to calculate basic CRDC functions also that are sine/cosine, inverse Tan.
References
[1]. Golub, G. H. and Van Loan, C. F. (1965), “Calculating the singular values and pseudo-inverse of a matrix”, SIAM J. Numer. Anal.
2, Ser. B (1965) pp 205–224
[2]. Sibul, L. H. (1984), “Application of Singular-value Decomposition to Adaptive beamforming”, IEEE Int. Conf. Acoustics, Speech,
and Signal Processing, 2:33.11.1-33.11.4
[3]. Andrews, H. C. and Patterson, C. L. (1976), “Singular Value Decompositions and Digital Image Processing”, IEEE Trans.
Acoustics, Speech, and Signal Processing, ASSP- 24(1), pp 26-53
[4]. Volder, J. (1959), “The CORDIC Trigonometric Computing Technique”, IRE Trans. Electronic Computers, EC-8(3), pp 330-334
[5]. Cavallaro, J.R. and Luk, F.T. (1988), “CORDIC Arithmetic for an SVD Processor”, Journal for Parallel and Distributed Computing,
vol. 5, pp. 271-290
[6]. Kota, K. (1991), “Architectural, Numerical and Implementation Issues in the VLSI Design of an Integrated CORDIC SVD
Processor”, Master's thesis, Rice University, Department of Electrical and Computer Engineering.
[7]. Luk, F. T. (), “Architectures for Computing Eigenvalues and SVDs”, Proc. SPIE Highly Parallel Signal Processing Architectures,
614.
[8]. Luk, F. T. (1985), “A Parallel Method for Computing the Generalized Singular Value Decomposition”, Journal of Parallel and
Distributed Computing, 2: pp 250-260
[9]. Luk, F. T. (1980), “Computing the Singular Value Decomposition on the ILLIAC IV”, ACM Transactions on Mathematical
Software, 6(4): pp 524-539
[10]. Meggitt, J. E. (1962), “Pseudo Division and Pseudo Multiplication Processes”, IBM Journal, pages 210-226
[11]. Speiser J. M. and Whitehouse, H. J. (1983), “A Review of Signal Processing with Systolic Arrays”, Proc. SPIE Real-Time Signal
Processing, 431(VI) pp 2-6
[12]. Ercegovac, M. D. and Lang, T. (1990), "Redundant and on-line CORDIC: application to matrix triangularization and SVD", IEEE
Trans. Computers, vol. 39, pp. 725-740

More Related Content

PDF
Hv2514131415
PDF
Bivariatealgebraic integerencoded arai algorithm for
PDF
Area Time Efficient Scaling Free Rotation Mode Cordic Using Circular Trajectory
PDF
Paper id 25201467
PDF
Reconfigurable CORDIC Low-Power Implementation of Complex Signal Processing f...
PDF
Bf36342346
PDF
Fpga sotcore architecture for lifting scheme revised
PDF
IRJET-ASIC Implementation for SOBEL Accelerator
Hv2514131415
Bivariatealgebraic integerencoded arai algorithm for
Area Time Efficient Scaling Free Rotation Mode Cordic Using Circular Trajectory
Paper id 25201467
Reconfigurable CORDIC Low-Power Implementation of Complex Signal Processing f...
Bf36342346
Fpga sotcore architecture for lifting scheme revised
IRJET-ASIC Implementation for SOBEL Accelerator

What's hot (18)

PDF
IIIRJET-Implementation of Image Compression Algorithm on FPGA
PDF
Non standard size image compression with reversible embedded wavelets
PDF
PDF
International Journal of Engineering Research and Development
PPTX
QUALITY ASSESSMENT FOR LIDAR POINT CLOUD REGISTRATION USING IN-SITU CONJUGATE...
PDF
IRJET- Implementation of Radix-16 and Binary 64 Division VLSI Realization...
PDF
“FIELD PROGRAMMABLE DSP ARRAYS” - A NOVEL RECONFIGURABLE ARCHITECTURE FOR EFF...
PDF
Accelerating Real Time Applications on Heterogeneous Platforms
PPTX
Vedic multiplier
PDF
Dynamic Texture Coding using Modified Haar Wavelet with CUDA
PDF
CORDIC Algorithm for WLAN
PDF
International Journal of Engineering Research and Development
PDF
Implementation of an arithmetic logic using area efficient carry lookahead adder
PDF
Design of High Speed and Low Power Veterbi Decoder for Trellis Coded Modulati...
PDF
D. Vulcanov: Symbolic Computation Methods in Cosmology and General Relativity...
DOCX
IEEE 2014 MATLAB IMAGE PROCESSING PROJECTS An efficient-parallel-approach-fo...
PDF
HIGH SPEED REVERSE CONVERTER FOR HIGH DYNAMIC RANGE MODULI SET
PDF
Fpga based efficient multiplier for image processing applications using recur...
IIIRJET-Implementation of Image Compression Algorithm on FPGA
Non standard size image compression with reversible embedded wavelets
International Journal of Engineering Research and Development
QUALITY ASSESSMENT FOR LIDAR POINT CLOUD REGISTRATION USING IN-SITU CONJUGATE...
IRJET- Implementation of Radix-16 and Binary 64 Division VLSI Realization...
“FIELD PROGRAMMABLE DSP ARRAYS” - A NOVEL RECONFIGURABLE ARCHITECTURE FOR EFF...
Accelerating Real Time Applications on Heterogeneous Platforms
Vedic multiplier
Dynamic Texture Coding using Modified Haar Wavelet with CUDA
CORDIC Algorithm for WLAN
International Journal of Engineering Research and Development
Implementation of an arithmetic logic using area efficient carry lookahead adder
Design of High Speed and Low Power Veterbi Decoder for Trellis Coded Modulati...
D. Vulcanov: Symbolic Computation Methods in Cosmology and General Relativity...
IEEE 2014 MATLAB IMAGE PROCESSING PROJECTS An efficient-parallel-approach-fo...
HIGH SPEED REVERSE CONVERTER FOR HIGH DYNAMIC RANGE MODULI SET
Fpga based efficient multiplier for image processing applications using recur...
Ad

Viewers also liked (20)

PDF
Use of Texting in Teaching and Its Effect on On-Time Homework Submission Rate...
PDF
Optimization of ZLD in Distillery Industry by Reverse Osmosis Process for Pre...
PDF
Experimental Investigations to Study the Air Flow Patterns on the Headlight D...
PDF
The Addition of Aluminum Nanoparticles to Polypropylene Increases Its Thermal...
PDF
Digital Heritage Documentation Via TLS And Photogrammetry Case Study
PDF
Enhancing the Mean Ratio Estimator for Estimating Population Mean Using Conve...
PDF
Distance Estimation to Image Objects Using Adapted Scale
PDF
Revision of the Rainfall Intensity Duration Frequency Curves for the City of ...
PDF
Sulphonamides: A Pharmaceutical Review
PDF
New Analytical Technique For The Determination Of Mercury (II) By Synergistic...
PDF
Software Architecture Design Recovery through Run-Time Source Code Collaborat...
PDF
A Study of Neural Network Learning-Based Recommender System
PDF
Method Data Envelopment Analysis for Estimating of Efficiency in Indication o...
PDF
Influence of Environmental Factors on the Production of Violacein Synthesized...
PDF
A Novel Polymeric Prodrugs Synthesized by Mechanochemical Solid-State Copolym...
PDF
A New Approach on the Log - Convex Orderings and Integral inequalities of the...
PDF
Oral Contraceptive Pills Use and Hypertension
PDF
Design and Prototype Development of Hybrid Ploughing, Seeding and Fertilizing...
PDF
Advanced Approach to Play a Atari Game
PDF
Analysis of G+4 Multistoried Building for Various Locations of Shear Wall Con...
Use of Texting in Teaching and Its Effect on On-Time Homework Submission Rate...
Optimization of ZLD in Distillery Industry by Reverse Osmosis Process for Pre...
Experimental Investigations to Study the Air Flow Patterns on the Headlight D...
The Addition of Aluminum Nanoparticles to Polypropylene Increases Its Thermal...
Digital Heritage Documentation Via TLS And Photogrammetry Case Study
Enhancing the Mean Ratio Estimator for Estimating Population Mean Using Conve...
Distance Estimation to Image Objects Using Adapted Scale
Revision of the Rainfall Intensity Duration Frequency Curves for the City of ...
Sulphonamides: A Pharmaceutical Review
New Analytical Technique For The Determination Of Mercury (II) By Synergistic...
Software Architecture Design Recovery through Run-Time Source Code Collaborat...
A Study of Neural Network Learning-Based Recommender System
Method Data Envelopment Analysis for Estimating of Efficiency in Indication o...
Influence of Environmental Factors on the Production of Violacein Synthesized...
A Novel Polymeric Prodrugs Synthesized by Mechanochemical Solid-State Copolym...
A New Approach on the Log - Convex Orderings and Integral inequalities of the...
Oral Contraceptive Pills Use and Hypertension
Design and Prototype Development of Hybrid Ploughing, Seeding and Fertilizing...
Advanced Approach to Play a Atari Game
Analysis of G+4 Multistoried Building for Various Locations of Shear Wall Con...
Ad

Similar to “ Implimentation of SD Processor Based On CRDC Algorithm ” (20)

PDF
Designing of cordic processor in verilog using xilinx ise simulator
PDF
Comparative analysis of multi stage cordic using micro rotation techniq
PDF
Comparative analysis of multi stage cordic using micro rotation techniq
PDF
Comparative analysis of multi stage cordic using micro rotation techniq
PDF
2016.04.21 - State of the Art
PDF
Comparative analysis of multi stage cordic using micro rotation technique
PDF
DSM Based low oversampling using SDR transmitter
PDF
Implementation of Rotation and Vectoring-Mode Reconfigurable CORDIC
PDF
Realization of Direct Digital Synthesis in Cordic Algorithm
PPTX
Presentation 7.pptx successive approximate register
DOCX
Novel design algorithm for low complexity programmable fir filters based on e...
PDF
FPGA Implementation of Optimized CIC Filter for Sample Rate Conversion in Sof...
PDF
Implementation of first order statistical processor on FPGA for feature extra...
PDF
Ax03303120316
PDF
Implementation of Efficiency CORDIC Algorithmfor Sine & Cosine Generation
PDF
FPGA Implementation of CORDIC Algorithm Architecture
PDF
Design and Implementation of Variable Radius Sphere Decoding Algorithm
PDF
J010234960
PDF
Titan X Research Paper
PDF
High speed cordic design for fixed angle of rotation
Designing of cordic processor in verilog using xilinx ise simulator
Comparative analysis of multi stage cordic using micro rotation techniq
Comparative analysis of multi stage cordic using micro rotation techniq
Comparative analysis of multi stage cordic using micro rotation techniq
2016.04.21 - State of the Art
Comparative analysis of multi stage cordic using micro rotation technique
DSM Based low oversampling using SDR transmitter
Implementation of Rotation and Vectoring-Mode Reconfigurable CORDIC
Realization of Direct Digital Synthesis in Cordic Algorithm
Presentation 7.pptx successive approximate register
Novel design algorithm for low complexity programmable fir filters based on e...
FPGA Implementation of Optimized CIC Filter for Sample Rate Conversion in Sof...
Implementation of first order statistical processor on FPGA for feature extra...
Ax03303120316
Implementation of Efficiency CORDIC Algorithmfor Sine & Cosine Generation
FPGA Implementation of CORDIC Algorithm Architecture
Design and Implementation of Variable Radius Sphere Decoding Algorithm
J010234960
Titan X Research Paper
High speed cordic design for fixed angle of rotation

Recently uploaded (20)

PDF
composite construction of structures.pdf
PPTX
Lecture Notes Electrical Wiring System Components
PPTX
M Tech Sem 1 Civil Engineering Environmental Sciences.pptx
DOCX
ASol_English-Language-Literature-Set-1-27-02-2023-converted.docx
PDF
Digital Logic Computer Design lecture notes
PPTX
web development for engineering and engineering
DOCX
573137875-Attendance-Management-System-original
PDF
keyrequirementskkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk
PPTX
UNIT 4 Total Quality Management .pptx
PPTX
OOP with Java - Java Introduction (Basics)
PPTX
MET 305 2019 SCHEME MODULE 2 COMPLETE.pptx
PDF
Well-logging-methods_new................
PPTX
Current and future trends in Computer Vision.pptx
PPT
Introduction, IoT Design Methodology, Case Study on IoT System for Weather Mo...
PDF
The CXO Playbook 2025 – Future-Ready Strategies for C-Suite Leaders Cerebrai...
PPTX
Sustainable Sites - Green Building Construction
PPT
introduction to datamining and warehousing
PPTX
Infosys Presentation by1.Riyan Bagwan 2.Samadhan Naiknavare 3.Gaurav Shinde 4...
PPTX
CYBER-CRIMES AND SECURITY A guide to understanding
PPTX
UNIT-1 - COAL BASED THERMAL POWER PLANTS
composite construction of structures.pdf
Lecture Notes Electrical Wiring System Components
M Tech Sem 1 Civil Engineering Environmental Sciences.pptx
ASol_English-Language-Literature-Set-1-27-02-2023-converted.docx
Digital Logic Computer Design lecture notes
web development for engineering and engineering
573137875-Attendance-Management-System-original
keyrequirementskkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk
UNIT 4 Total Quality Management .pptx
OOP with Java - Java Introduction (Basics)
MET 305 2019 SCHEME MODULE 2 COMPLETE.pptx
Well-logging-methods_new................
Current and future trends in Computer Vision.pptx
Introduction, IoT Design Methodology, Case Study on IoT System for Weather Mo...
The CXO Playbook 2025 – Future-Ready Strategies for C-Suite Leaders Cerebrai...
Sustainable Sites - Green Building Construction
introduction to datamining and warehousing
Infosys Presentation by1.Riyan Bagwan 2.Samadhan Naiknavare 3.Gaurav Shinde 4...
CYBER-CRIMES AND SECURITY A guide to understanding
UNIT-1 - COAL BASED THERMAL POWER PLANTS

“ Implimentation of SD Processor Based On CRDC Algorithm ”

  • 1. International Journal of Engineering Science Invention ISSN (Online): 2319 – 6734, ISSN (Print): 2319 – 6726 www.ijesi.org ||Volume 6 Issue 1|| January 2017 || PP. 38-39 www.ijesi.org 38 | Page “ Implimentation of SD Processor Based On CRDC Algorithm ” Raj Kumar Verma#1 , Mrs.Akanksha Awasthi#2 #1 Research Scholar, C.V.R.U, Bilaspurc.g. India. #2 Asst.Prof.C.V.R.U., Bilaspurc.g. India. Abstract: In Digital Signal Processing (DSP) there are many complex algorithms for which an efficient hardware implementation is required in real time applications. One such complex algorithm is Singular-value Decomposition (SD) which is an important algorithm with applications in varied domains of signal processing such as direction estimation, spectrum analysis and systems identification. It is a generalized extension to the eigen-decomposition for non-square matrices and is hence of great importance, particularly for subspace based algorithms in signal processing. But SD is known to be a very complicated algorithm with computational complexity ~O(N3 ) (for a NxN square matrix). For real-time computation of such a complex algorithm the use of a parallel and direct mapped hardware solution is indeed desired. Hardware and software resources: 1. OS: Windows 9x or upper 2. RAM : Minimum 512 MB 3. Programming Language: XILINX 8.6 or upper ( VERILOG) I. Introduction In this chapter basic idea of the project undertaken is explained. Starting from the aim of the project, Importance and the problem associated with the design is also discussed. Finally line of approach to solve the problem is presented and organization of rest of the report is briefly shown. Importance The Singular Value Decomposition (SD) is an important matrix factorization with applications in signal processing, image processing and robotics. It is generally acknowledged that the SD is the only generally reliable method for determining the rank of a matrix numerically. The SD is a very useful tool, for example, in analyzing data matrices from sensor arrays for adaptive beam forming , and low rank approximations to matrices in image enhancement . The wide variety of applications for the SD coupled with its computational complexity justifies dedicating hardware to this computation. Designed 2X2 CRDC based SD processor can be used as basic building block for an array of processors of N X N matrix. Problem statement The importance of hardware implementation is now known but there are certain complexities that come on the way of implementation. In a typical computer algorithm for the SD, the sines and cosines of the rotation angles are computed through formulas that require division and square root operations. Also the explicit angles are not required and only thesines and cosines are computed. The rotations are then applied to the 2x 2 sub-matrix using standard matrix multiplication techniques. But in this method many time-consuming operations such as multiplication, division, and square rootare needed. The SD is more closely mapped ontohardware through the useof the CRDC algorithms. CRDC algorithm provides an iterative scheme, consisting of simple addition and binary shift operations, to compute trigonometric values to any desired precision. The CRDC algorithms can provide the calculation of vector rotation and inverse tangent without costly hardware that proves very beneficial for implementation of SD processor. Line of approach The complexity of SD computation and its real-time applications were the inspiration behind efficient hardware implementation. In SD processor time consuming operations like multiplication, square root and division are eliminated by the use of CRDC algorithm. CRDC algorithm only requires simple shift and additions which efficiently maps SD processor.
  • 2. “ Implimentation of SD Processor Based On CRDC Algorithm ” www.ijesi.org 39 | Page Basic CRDC architecture An iterative CRDC architecture can be obtained simply by duplicating each of the three difference equations in hardware as shown in Figure. Figure: Basic CRDC hardware The decision function, i d is driven by the sign of the y or z register depending on whether it is operated in rotation ( i z ) or vectoring mode ( y i ). In operation, the initial values are loaded via multiplexers into the x, y and z registers. Then on each of the next n clock cycles, the values from the registers are passed through the shifters and adders-subtractors and the results placed back in the registers. The shifters are modified on each iteration to cause the desired shift for the iteration. Likewise, the ROM (Look-up table) address is incremented on each iteration so that the appropriate elementary angle value is presented to the z adder-subtractor. On the last iteration, the results are read directly from the adder-subtractors. A simple state machine is required to keep track of the current iteration, and to select the degree of shift and ROM address for each iteration.As shown above, CRDC iteration requires 2 shifters, 1 table lookup and 3 adders. For n bits of precision, n iterations are needed. II. Conclusion In linear algebra, the singular value decomposition (SD) is an important factorization of a rectangular real or complex matrix, with several applications in signal processing and statistics. A 2X2 CRDC based SD processor was designed and successfully simulated. CRDC algorithm was used to simplify the computational complexities of SD calculation. Apart from calculating the SD of a 2x2 matrix the processor is able to calculate basic CRDC functions also that are sine/cosine, inverse Tan. References [1]. Golub, G. H. and Van Loan, C. F. (1965), “Calculating the singular values and pseudo-inverse of a matrix”, SIAM J. Numer. Anal. 2, Ser. B (1965) pp 205–224 [2]. Sibul, L. H. (1984), “Application of Singular-value Decomposition to Adaptive beamforming”, IEEE Int. Conf. Acoustics, Speech, and Signal Processing, 2:33.11.1-33.11.4 [3]. Andrews, H. C. and Patterson, C. L. (1976), “Singular Value Decompositions and Digital Image Processing”, IEEE Trans. Acoustics, Speech, and Signal Processing, ASSP- 24(1), pp 26-53 [4]. Volder, J. (1959), “The CORDIC Trigonometric Computing Technique”, IRE Trans. Electronic Computers, EC-8(3), pp 330-334 [5]. Cavallaro, J.R. and Luk, F.T. (1988), “CORDIC Arithmetic for an SVD Processor”, Journal for Parallel and Distributed Computing, vol. 5, pp. 271-290 [6]. Kota, K. (1991), “Architectural, Numerical and Implementation Issues in the VLSI Design of an Integrated CORDIC SVD Processor”, Master's thesis, Rice University, Department of Electrical and Computer Engineering. [7]. Luk, F. T. (), “Architectures for Computing Eigenvalues and SVDs”, Proc. SPIE Highly Parallel Signal Processing Architectures, 614. [8]. Luk, F. T. (1985), “A Parallel Method for Computing the Generalized Singular Value Decomposition”, Journal of Parallel and Distributed Computing, 2: pp 250-260 [9]. Luk, F. T. (1980), “Computing the Singular Value Decomposition on the ILLIAC IV”, ACM Transactions on Mathematical Software, 6(4): pp 524-539 [10]. Meggitt, J. E. (1962), “Pseudo Division and Pseudo Multiplication Processes”, IBM Journal, pages 210-226 [11]. Speiser J. M. and Whitehouse, H. J. (1983), “A Review of Signal Processing with Systolic Arrays”, Proc. SPIE Real-Time Signal Processing, 431(VI) pp 2-6 [12]. Ercegovac, M. D. and Lang, T. (1990), "Redundant and on-line CORDIC: application to matrix triangularization and SVD", IEEE Trans. Computers, vol. 39, pp. 725-740