DEV Community

Cover image for Caching an Inverse of a Matrix Using R
Masteramuk
Masteramuk

Posted on

Caching an Inverse of a Matrix Using R

Matrix inversion is usually a costly computation and there may be some benefit to caching the inverse of a matrix rather than computing it repeatedly.

This is an example code to cache an inverse of a matrix which then called if:

  1. The new matrix computation similar to the previous matrix
  2. Matrix is invertible — Read Invertible Matrix
  3. The structure and content of the new matrix is equal to previous The program consists of two functions

A function to cache the matrix and the inverse. It is similar to getter and setter in many programming languages like Java and C.
A function which will be call for the computational process of the matrix

Here are parts of the code:

Setting the matrix

invMatrix <- NULL #define the set function for the matrix set <- function(y) { x <<- y invMatrix <<- NULL } 

Setting the inverse matrix

setInvMatrix <- function(invM){ invMatrix <<- invM ##assign the inverse matrix to the environment var } ##define the function for returning the inverse matrix getInvMatrix <- function(){ invMatrix ##the environment var inverse matrix is return } 

Verification on the matrix

if (!is.null(invM)){ ## check if the return invMatrix is identical if ( identical( x$get() %*% invM, invM %*% x$get() ) ){ ## get it from the cache and skips the computation. print("getting cached data") return(invM) } } 

Setting the new inverse

data <- x$get() invM <- solve(data, ...) # sets the value of the inverse in the cache via the setinv function. x$setInvMatrix(invM) 

Complete source code is reachable at https://github.com/masteramuk/ProgrammingAssignment2

If you like the post, https://www.buymeacoffee.com/masteramuk :)

Top comments (0)