OneStopGate.Com
OnestopGate   OnestopGate
   Sunday, November 17, 2024 Login  
OnestopGate
Home | Overview | Syllabus | Tutorials | FAQs | Downloads | Recommended Websites | Advertise | Payments | Contact Us | Forum
OneStopGate

GATE Resources
Gate Articles
Gate Books
Gate Colleges 
Gate Downloads 
Gate Faqs
Gate Jobs
Gate News 
Gate Sample Papers
Training Institutes

GATE Overview
Overview
GATE Eligibility
Structure Of GATE
GATE Coaching Centers
Colleges Providing M.Tech/M.E.
GATE Score
GATE Results
PG with Scholarships
Article On GATE
Admission Process For M.Tech/ MCP-PhD
GATE Topper 2012-13
GATE Forum




GATE 2025 Exclusive
Organizing Institute
Important Dates
How to Apply
Discipline Codes
GATE 2025 Exam Structure

GATE 2025 Syllabus
Aerospace Engg..
Agricultural Engg..
Architecture and Planning
Chemical Engg..
Chemistry
Civil Engg..
Computer Science / IT
Electronics & Communication Engg..
Electrical Engg..
Engineering Sciences
Geology and Geophysics
Instrumentation Engineering
Life Sciences
Mathematics
Mechanical Engg..
Metallurgical Engg..
Mining Engg..
Physics
Production & Industrial Engg..
Pharmaceutical Sciences
Textile Engineering and Fibre Science

GATE Study Material
Aerospace Engg..
Agricultural Engg..
Chemical Engg..
Chemistry
Civil Engg..
Computer Science / IT
Electronics & Communication Engg..
Electrical Engg..
Engineering Sciences
Instrumentation Engg..
Life Sciences
Mathematics
Mechanical Engg..
Physics
Pharmaceutical Sciences
Textile Engineering  and Fibre Science

GATE Preparation
GATE Pattern
GATE Tips N Tricks
Compare Evaluation
Sample Papers 
Gate Downloads 
Experts View

CEED 2013
CEED Exams
Eligibility
Application Forms
Important Dates
Contact Address
Examination Centres
CEED Sample Papers

Discuss GATE
GATE Forum
Exam Cities
Contact Details
Bank Details

Miscellaneous
Advertisment
Contact Us


Home » GATE Study Material » Mathematics » Numerical Analysis » Linear Equations » PA = LU Factorization with Pivoting

PA = LU Factorization with Pivoting

Looking for GATE Preparation Material? Join & Get here now!

** Gate 2013 Question Papers.. ** CEED 2013 Results.. ** Gate 2013 Question Papers With Solutions.. ** GATE 2013 CUT-OFFs.. ** GATE 2013 Results.. **

PA = LU Factorization with Pivoting

PA = LU Factorization with Pivoting

     

Background

Definition ( LU-Factorization ).  The nonsingular matrix A has an LU-factorization if it can be expressed as the product of a lower-triangular matrix L and an upper triangular matrix U:  

          [Graphics:Images/LUFactorMod_gr_1.gif].  

When this is possible we say that A has an LU-decomposition.  It turns out that this factorization (when it exists) is not unique.  If L has 1's on it's diagonal, then it is called a Doolittle factorization.  If U has 1's on its diagonal, then it is called a Crout factorization.  When  [Graphics:Images/LUFactorMod_gr_2.gif], it is called a Cholesky decomposition.  In this module we will develop an algorithm that produces a Doolittle factorization.

 

Theorem ( LU Factorization with NO pivoting).  If row interchanges arenot needed to solve the linear system AX = B, then A has the LU factorization (illustrated with 4�4 matrices).  

        [Graphics:Images/LUFactorMod_gr_3.gif]= [Graphics:Images/LUFactorMod_gr_4.gif][Graphics:Images/LUFactorMod_gr_5.gif].
        
Remark 1.  This is not a linear system.
Remark 2.  The easy solution uses row vectors and is a modification of limited Gauss Jordan elimination.
Remark 3.  A sufficient condition for the factorization to exist is that all principal minors of A are nonsingular.

 

Mathematica Subroutine (Limited Gauss-Jordan Elimination).

[Graphics:Images/LUFactorMod_gr_6.gif]

Mathematica Subroutine (LandU).

[Graphics:Images/LUFactorMod_gr_7.gif]

 

Theorem (A = LU;  Factorization with NO Pivoting).  Assume that A has a Doolittle, Crout or Cholesky factorization.  The solution X to the linear system  [Graphics:Images/LUFactorMod_gr_64.gif], is found in three steps:  

    
1.  Construct the matrices  [Graphics:Images/LUFactorMod_gr_65.gif], if possible.  
    
2.  Solve  [Graphics:Images/LUFactorMod_gr_66.gif]  for  [Graphics:Images/LUFactorMod_gr_67.gif]  using forward substitution.
    
3.  Solve  [Graphics:Images/LUFactorMod_gr_68.gif]  for  [Graphics:Images/LUFactorMod_gr_69.gif]  using back substitution.

 

    The above theorem assumes that there are no row interchanges.  We have seen in Example 3 an example of a nonsingular matrix A could not be directly factored as A = LU.  If row interchanges are permitted then a factorization of a "permuted matrix" will be obtained.  A permutation of the first n positive integers  [Graphics:Images/LUFactorMod_gr_70.gif].  is an arrangement  [Graphics:Images/LUFactorMod_gr_71.gif]  of these integers in a definite order.  For example,  [Graphics:Images/LUFactorMod_gr_72.gif]  is a permutation of the five integers  [Graphics:Images/LUFactorMod_gr_73.gif].  The standard base vectors  [Graphics:Images/LUFactorMod_gr_74.gif]  for  [Graphics:Images/LUFactorMod_gr_75.gif]  are used in the next definition.  

 

Definition ( Permutation Matrix).  An n�n permutation matrix P is a matrix with precisely one entry whose value is "1" in each column and row, and all of whose other entries are "0."  The rows of  P  are a permutation of the rows of the identity matrix and  P  can be written as  

        
[Graphics:Images/LUFactorMod_gr_76.gif].

The elements of  [Graphics:Images/LUFactorMod_gr_77.gif]  have the form  

        [Graphics:Images/LUFactorMod_gr_78.gif]  

 

Theorem.  Suppose that  [Graphics:Images/LUFactorMod_gr_79.gif]  is a permutation matrix.  The product  PA  is a new matrix whose rows consists of the rows of  A  rearranged in the new order  [Graphics:Images/LUFactorMod_gr_80.gif].  

 

For example, the permutation matrix  [Graphics:Images/LUFactorMod_gr_81.gif]  will interchange rows 1 and 2 and also interchange rows 3 and 4.  

[Graphics:Images/LUFactorMod_gr_82.gif]
[Graphics:Images/LUFactorMod_gr_83.gif]

Theorem.  If  A  is a nonsingular matrix, then there exists a permutation matrix  P  so that  PA  has an LU-factorization

        PA = LU.  

 

Theorem (PA = LU;  Factorization with Pivoting).  Given that A is nonsingular.  The solution X to the linear system  [Graphics:Images/LUFactorMod_gr_84.gif], is found in four steps:  

    
1.  Construct the matrices  [Graphics:Images/LUFactorMod_gr_85.gif].  
    
2.  Compute the column vector   [Graphics:Images/LUFactorMod_gr_86.gif].
    
3.  Solve  [Graphics:Images/LUFactorMod_gr_87.gif]  for  [Graphics:Images/LUFactorMod_gr_88.gif]  using forward substitution.
    
4.  Solve  [Graphics:Images/LUFactorMod_gr_89.gif]  for  [Graphics:Images/LUFactorMod_gr_90.gif]  using back substitution.

 

Mathematica Subroutine (LUfactor).
Remark.  The matrix  A  is a global variable and elements are changed when the  LUfactor  subroutine is executed.  Hence, it is important to save a copy of  A  in the variable  A0.
Use one of the following two versions of the  LUfactor  subroutine and the  SolveLU  subroutine.  The first version of  LUfactor  uses parallel programming and "row operations" instead of loops.

[Graphics:Images/LUFactorMod_gr_91.gif]

This is the second version of  LUfactor  and it uses more loops and traditional programming.

Mathematica Subroutine (LUfactor).

[Graphics:Images/LUFactorMod_gr_92.gif]

Use the subroutine SolveLU which is similar to the forward substitution and back substitution subroutines.

Mathematica Subroutine (SolveLU).

[Graphics:Images/LUFactorMod_gr_93.gif]

Remark.  Everything has been carefully set up so that  L,  U, and  P  can all be studied.

 

Application to Polynomial Curve Fitting

Theorem ( Least-Squares Polynomial Curve Fitting ).
Given the  [Graphics:Images/LUFactorMod_gr_188.gif]  data points  [Graphics:Images/LUFactorMod_gr_189.gif],  the least squares polynomial of degree  m  of the form  

        [Graphics:Images/LUFactorMod_gr_190.gif]

that fits the n data points is obtained by solving the following linear system


        [Graphics:Images/LUFactorMod_gr_191.gif][Graphics:Images/LUFactorMod_gr_192.gif]  

for the m+1 coefficients [Graphics:Images/LUFactorMod_gr_193.gif].  These equations are referred to as the "normal equations".

 



Discussion Center

Discuss/
Query

Papers/
Syllabus

Feedback/
Suggestion

Yahoo
Groups

Sirfdosti
Groups

Contact
Us

MEMBERS LOGIN
  
Email ID:
Password:

  Forgot Password?
 New User? Register!

INTERVIEW EBOOK
Get 9,000+ Interview Questions & Answers in an eBook. Interview Question & Answer Guide
  • 9,000+ Interview Questions
  • All Questions Answered
  • 5 FREE Bonuses
  • Free Upgrades
GATE RESOURCES
 
  • Gate Books
  • Training Institutes
  • Gate FAQs
  • GATE BOOKS
     
  • Mechanical Engineeering Books
  • Robotics Automations Engineering Books
  • Civil Engineering Books
  • Chemical Engineering Books
  • Environmental Engineering Books
  • Electrical Engineering Books
  • Electronics Engineering Books
  • Information Technology Books
  • Software Engineering Books
  • GATE Preparation Books
  • Exciting Offers



    GATE Exam, Gate 2009, Gate Papers, Gate Preparation & Related Pages


    GATE Overview | GATE Eligibility | Structure Of GATE | GATE Training Institutes | Colleges Providing M.Tech/M.E. | GATE Score | GATE Results | PG with Scholarships | Article On GATE | GATE Forum | GATE 2009 Exclusive | GATE 2009 Syllabus | GATE Organizing Institute | Important Dates for GATE Exam | How to Apply for GATE | Discipline / Branch Codes | GATE Syllabus for Aerospace Engineering | GATE Syllabus for Agricultural Engineering | GATE Syllabus for Architecture and Planning | GATE Syllabus for Chemical Engineering | GATE Syllabus for Chemistry | GATE Syllabus for Civil Engineering | GATE Syllabus for Computer Science / IT | GATE Syllabus for Electronics and Communication Engineering | GATE Syllabus for Engineering Sciences | GATE Syllabus for Geology and Geophysics | GATE Syllabus for Instrumentation Engineering | GATE Syllabus for Life Sciences | GATE Syllabus for Mathematics | GATE Syllabus for Mechanical Engineering | GATE Syllabus for Metallurgical Engineering | GATE Syllabus for Mining Engineering | GATE Syllabus for Physics | GATE Syllabus for Production and Industrial Engineering | GATE Syllabus for Pharmaceutical Sciences | GATE Syllabus for Textile Engineering and Fibre Science | GATE Preparation | GATE Pattern | GATE Tips & Tricks | GATE Compare Evaluation | GATE Sample Papers | GATE Downloads | Experts View on GATE | CEED 2009 | CEED 2009 Exam | Eligibility for CEED Exam | Application forms of CEED Exam | Important Dates of CEED Exam | Contact Address for CEED Exam | CEED Examination Centres | CEED Sample Papers | Discuss GATE | GATE Forum of OneStopGATE.com | GATE Exam Cities | Contact Details for GATE | Bank Details for GATE | GATE Miscellaneous Info | GATE FAQs | Advertisement on GATE | Contact Us on OneStopGATE |
    Copyright © 2024. One Stop Gate.com. All rights reserved Testimonials |Link To Us |Sitemap |Privacy Policy | Terms and Conditions|About Us
    Our Portals : Academic Tutorials | Best eBooksworld | Beyond Stats | City Details | Interview Questions | India Job Forum | Excellent Mobiles | Free Bangalore | Give Me The Code | Gog Logo | Free Classifieds | Jobs Assist | Interview Questions | One Stop FAQs | One Stop GATE | One Stop GRE | One Stop IAS | One Stop MBA | One Stop SAP | One Stop Testing | Web Hosting | Quick Site Kit | Sirf Dosti | Source Codes World | Tasty Food | Tech Archive | Software Testing Interview Questions | Free Online Exams | The Galz | Top Masala | Vyom | Vyom eBooks | Vyom International | Vyom Links | Vyoms | Vyom World
    C Interview Questions | C++ Interview Questions | Send Free SMS | Placement Papers | SMS Jokes | Cool Forwards | Romantic Shayari