FICO
FICO Xpress Optimization Examples Repository
FICO Optimization Community FICO Xpress Optimization Home
Back to examples browserPrevious exampleNext example

Invert a Mosel matrix with NumPy

Description
Invert a Mosel matrix in Python using NumPy.

Further explanation of this example: see Chapter 'Python' in the 'Mosel Language Reference Manual'

python3_invert_matrix.zip[download all files]

Source Files

Data Files





invert_matrix.mos

(!*********************************************************
   Mosel Python Example Problems
   =============================

   file invert_matrix.mos
   ``````````````````````
   Invert a Mosel matrix in Python using pandas and NumPy.

   !!! This example requires an installation of Python 3, see
   !!! chapter 'python3' of the 'Mosel Language Reference' for
   !!! compatible versions and setup instructions.

  (c) 2018 Fair Isaac Corporation
      author: J. Müller
*********************************************************!)
model "invert_matrix"
  options noimplicit

  uses "python3"
  uses "mmjobs"

  ! Input data
  declarations
    I, J: range
    A, A_inverse: array(I, J) of real
  end-declarations

  procedure show_matrix(B: array(I: range, J: range) of real)
    forall (i in I) do
      forall (j in J) do
        write(B(i, j), " ")
      end-do
      writeln
    end-do
  end-procedure

  writeln("Run Python script that defines invert_matrix function.")
  writeln("Python DSO version:     ", getdsoprop("python3", PROP_VERSION))
  pyinitpandas
  pyrun("invert_matrix.py")

  I := 0..2
  J := 0..2
  A :: [1,0,3,
        0,1,2,
        0,0,1]

  writeln("Mosel matrix A:")
  show_matrix(A)

  ! Optional: Delete old values from result matrix before calling function.
  delcell(A_inverse)

  writeln("Invert matrix with NumPy.")
  pycall("invert_matrix", A_inverse, A)

  writeln("Matrix A_inverse:")
  show_matrix(A_inverse)
end-model

Back to examples browserPrevious exampleNext example