Skip to main content
added Mike's formulation as an example
Source Link
kglr
  • 403.4k
  • 18
  • 501
  • 959

For a general square matrix m and arbitrary partition of it into conformable parts m={{a,b},{c,d}} (i.e., a and d are square matrices, and b and c have appropriate dimensions),

partitioned matrix

the formula for the inverse (which can be found in, for example, Review of Matrix Algebra) is

 m={{a,b},{c,d}}; e = d - c.Inverse[a].b; minv={{Inverse[a] + Inverse[a].b.Inverse[e].c.Inverse[a], -Inverse[a].b.Inverse[e]}, {-Inverse[e].c.Inverse[a], Inverse[e]}}; 

enter image description here

This can be used to define a function that takes a matrix and a partition specification and does the needed type and conformability checks before applying the formula.

The formula You can be found inuse @Mike's formulation (and, for exampleoff course, Review of Matrix Algebra.heed his warnings) as a first step:

Clear[sInverse2]; sInverse2[{{a_, b_}, {c_, d_}}] := With[{e = d - c.Inverse[a].b}, {{Inverse[a] + Inverse[a].b.Inverse[e].c.Inverse[a], -Inverse[a].b.Inverse[e]}, {-Inverse[e].c.Inverse[a], Inverse[e]}}] 

For a general square matrix m and arbitrary partition of it into conformable parts m={{a,b},{c,d}} (i.e., a and d are square matrices, and b and c have appropriate dimensions),

partitioned matrix

the formula for the inverse is

 m={{a,b},{c,d}}; e = d - c.Inverse[a].b; minv={{Inverse[a] + Inverse[a].b.Inverse[e].c.Inverse[a], -Inverse[a].b.Inverse[e]}, {-Inverse[e].c.Inverse[a], Inverse[e]}}; 

enter image description here

This can be used to define a function that takes a matrix and a partition specification and does the needed type and conformability checks before applying the formula.

The formula can be found in, for example, Review of Matrix Algebra.

For a general square matrix m and arbitrary partition of it into conformable parts m={{a,b},{c,d}} (i.e., a and d are square matrices, and b and c have appropriate dimensions),

partitioned matrix

the formula for the inverse (which can be found in, for example, Review of Matrix Algebra) is

 m={{a,b},{c,d}}; e = d - c.Inverse[a].b; minv={{Inverse[a] + Inverse[a].b.Inverse[e].c.Inverse[a], -Inverse[a].b.Inverse[e]}, {-Inverse[e].c.Inverse[a], Inverse[e]}}; 

enter image description here

This can be used to define a function that takes a matrix and a partition specification and does the needed type and conformability checks before applying the formula. You can use @Mike's formulation (and, off course, heed his warnings) as a first step:

Clear[sInverse2]; sInverse2[{{a_, b_}, {c_, d_}}] := With[{e = d - c.Inverse[a].b}, {{Inverse[a] + Inverse[a].b.Inverse[e].c.Inverse[a], -Inverse[a].b.Inverse[e]}, {-Inverse[e].c.Inverse[a], Inverse[e]}}] 
edited body
Source Link
kglr
  • 403.4k
  • 18
  • 501
  • 959

For a general square matrix m and arbitrary partition of it into conformable parts m={{a,b},{c,d}} (i.e., a and d are square matrices, and cb and dc have appropriate dimensions),

partitioned matrix

the formula for the inverse is

 m={{a,b},{c,d}}; e = d - c.Inverse[a].b; minv={{Inverse[a] + Inverse[a].b.Inverse[e].c.Inverse[a], -Inverse[a].b.Inverse[e]}, {-Inverse[e].c.Inverse[a], Inverse[e]}}; 

enter image description here

This can be used to define a function that takes a matrix and a partition specification and does the needed type and conformability checks before applying the formula.

The formula can be found in, for example, Review of Matrix Algebra.

For a general square matrix m and arbitrary partition of it into conformable parts m={{a,b},{c,d}} (i.e., a and d are square matrices, and c and d have appropriate dimensions),

partitioned matrix

the formula for the inverse is

 m={{a,b},{c,d}}; e = d - c.Inverse[a].b; minv={{Inverse[a] + Inverse[a].b.Inverse[e].c.Inverse[a], -Inverse[a].b.Inverse[e]}, {-Inverse[e].c.Inverse[a], Inverse[e]}}; 

enter image description here

This can be used to define a function that takes a matrix and a partition specification and does the needed type and conformability checks before applying the formula.

The formula can be found in, for example, Review of Matrix Algebra.

For a general square matrix m and arbitrary partition of it into conformable parts m={{a,b},{c,d}} (i.e., a and d are square matrices, and b and c have appropriate dimensions),

partitioned matrix

the formula for the inverse is

 m={{a,b},{c,d}}; e = d - c.Inverse[a].b; minv={{Inverse[a] + Inverse[a].b.Inverse[e].c.Inverse[a], -Inverse[a].b.Inverse[e]}, {-Inverse[e].c.Inverse[a], Inverse[e]}}; 

enter image description here

This can be used to define a function that takes a matrix and a partition specification and does the needed type and conformability checks before applying the formula.

The formula can be found in, for example, Review of Matrix Algebra.

added hyperlink for ref
Source Link
kglr
  • 403.4k
  • 18
  • 501
  • 959

For a general square matrix m and arbitrary partition of it into conformable parts m={{a,b},{c,d}} (i.e., a and d are square matrices, and c and d have appropriate dimensions),

partitioned matrix

the formula for the inverse is

 m={{a,b},{c,d}}; e = d - c.Inverse[a].b; minv={{Inverse[a] + Inverse[a].b.Inverse[e].c.Inverse[a], -Inverse[a].b.Inverse[e]}, {-Inverse[e].c.Inverse[a], Inverse[e]}}; 

partitioned matrix inverseenter image description here

This can be used to define a function that takes a matrix and a partition specification and does the needed type and conformability checks before applying the formula.

The formula can be found in, for example, Review of Matrix Algebra.

For a general square matrix m and arbitrary partition of it into conformable parts m={{a,b},{c,d}} (i.e., a and d are square matrices, and c and d have appropriate dimensions),

partitioned matrix

the formula for the inverse is

 m={{a,b},{c,d}}; e = d - c.Inverse[a].b; minv={{Inverse[a] + Inverse[a].b.Inverse[e].c.Inverse[a], -Inverse[a].b.Inverse[e]}, {-Inverse[e].c.Inverse[a], Inverse[e]}}; 

partitioned matrix inverse

For a general square matrix m and arbitrary partition of it into conformable parts m={{a,b},{c,d}} (i.e., a and d are square matrices, and c and d have appropriate dimensions),

partitioned matrix

the formula for the inverse is

 m={{a,b},{c,d}}; e = d - c.Inverse[a].b; minv={{Inverse[a] + Inverse[a].b.Inverse[e].c.Inverse[a], -Inverse[a].b.Inverse[e]}, {-Inverse[e].c.Inverse[a], Inverse[e]}}; 

enter image description here

This can be used to define a function that takes a matrix and a partition specification and does the needed type and conformability checks before applying the formula.

The formula can be found in, for example, Review of Matrix Algebra.

Source Link
kglr
  • 403.4k
  • 18
  • 501
  • 959
Loading