C++ Program to Perform Matrix Division

Posted in C++ Tutorial

Updated on Nov 27, 2024

By Mari Selvan

👁️ 108 - Views

⏳ 4 mins

💬 1 Comment

Photo Credit to CodeToFun

🙋 Introduction

Matrix operations are fundamental in linear algebra and various scientific and engineering applications. One essential matrix operation is division.

In this tutorial, we'll explore a C++ program that performs matrix division, providing a foundational understanding of how to handle matrices in a programming context.

📄 Example

Let's delve into the C++ code that performs matrix division.

matrixDivision.cpp

Copied

#include <iostream>
#include <iomanip>
#include <vector>

// Function to print a matrix
void printMatrix(const std::vector < std::vector < double >> & matrix) {
  for (const auto & row: matrix) {
    for (double value: row) {
      std::cout << std::fixed << std::setprecision(2) << value << "\t";
    }
    std::cout << std::endl;
  }
}

// Function to perform matrix division
std::vector < std::vector < double >> matrixDivision(const std::vector < std::vector < double >> & matrixA,
  const std::vector < std::vector < double >> & matrixB) {
  int rows = matrixA.size();
  int cols = matrixA[0].size();

  std::vector < std::vector < double >> result(rows, std::vector < double > (cols, 0.0));

  // Check if matrixB is invertible
  // Additional logic for inverse calculation is required here

  // Perform matrix division
  for (int i = 0; i < rows; ++i) {
    for (int j = 0; j < cols; ++j) {
      result[i][j] = matrixA[i][j] / matrixB[i][j];
    }
  }

  return result;
}

// Driver program
int main() {
  // Define matrices A and B
  std::vector < std::vector < double >> matrixA = {{4.0, 8.0},{2.0, 6.0}};
  std::vector < std::vector < double >> matrixB = {{2.0, 4.0},{1.0, 3.0}};

  // Call the function to perform matrix division
  std::vector < std::vector < double >> result = matrixDivision(matrixA, matrixB);

  // Print the result
  std::cout << "Result of matrix division:" << std::endl;
  printMatrix(result);

  return 0;
}

💻 Testing the Program

To test the program with different matrices, replace the values of matrixA and matrixB in the main function.

Output

Result of matrix division:
2.00   2.00
2.00   2.00

Compile and run the program to see the result of matrix division.

🧠 How the Program Works

The program defines a function printMatrix to print a matrix for better visualization.
The matrixDivision function checks if matrix B is invertible (additional logic for inverse calculation is required).
It performs matrix division element-wise and prints the result.

🧐 Understanding the Concept of Matrix Division

Matrix division is not a straightforward operation like addition or multiplication.

Matrix division involves multiplying one matrix by the inverse of another. Ensure that the second matrix is invertible, and additional logic for inverse calculation may be required.

For two matrices A and B, the division A/B is equivalent to A * B^(-1), where B^(-1) is the inverse of matrix B.