vector_routines.f90

Go to the documentation of this file.

!/*************************************************************************
! *
! *  Project
! *                         _____ _____  __  __ _____
! *                        / ____|  __ \|  \/  |  __ \
! *  ___  _ __   ___ _ __ | |  __| |__) | \  / | |__) |
! * / _ \| '_ \ / _ \ '_ \| | |_ |  ___/| |\/| |  ___/
! *| (_) | |_) |  __/ | | | |__| | |    | |  | | |
! * \___/| .__/ \___|_| |_|\_____|_|    |_|  |_|_|
! *      | |
! *      |_|
! *
! * Copyright (C) Akiel Aries, <akiel@akiel.org>, et al.
! *
! * This software is licensed as described in the file LICENSE, which
! * you should have received as part of this distribution. The terms
! * among other details are referenced in the official documentation
! * seen here : https://akielaries.github.io/openGPMP/ along with
! * important files seen in this project.
! *
! * You may opt to use, copy, modify, merge, publish, distribute
! * and/or sell copies of the Software, and permit persons to whom
! * the Software is furnished to do so, under the terms of the
! * LICENSE file. As this is an Open Source effort, all implementations
! * must be of the same methodology.
! *
! *
! *
! * This software is distributed on an AS IS basis, WITHOUT
! * WARRANTY OF ANY KIND, either express or implied.
! *
! ************************************************************************/
! vector_routines.f90
 
SUBROUTINE mtx_add_routine_float_(A, B, C, mtx_size) bind(C)
    USE :: iso_fortran_env
    USE :: iso_c_binding
 
    INTEGER, INTENT(IN) :: mtx_size
    REAL(KIND=c_float), DIMENSION(mtx_size, mtx_size), INTENT(IN) :: a, b
    REAL(KIND=c_float), DIMENSION(mtx_size, mtx_size), INTENT(OUT) :: c
 
    c = a + b
END SUBROUTINE mtx_add_routine_float_
 
SUBROUTINE mtx_add_routine_int_(A, B, C, mtx_size) bind(C)
    USE :: iso_fortran_env
    USE :: iso_c_binding
 
    INTEGER, INTENT(IN) :: mtx_size
    INTEGER(C_INT), DIMENSION(mtx_size, mtx_size), INTENT(IN) :: A, B
    INTEGER(C_INT), DIMENSION(mtx_size, mtx_size), INTENT(OUT) :: C
 
    c = a + b
END SUBROUTINE mtx_add_routine_int_
 
SUBROUTINE mtx_mult(matrix1, matrix2, result, nrows1, ncols1, ncols2)
    implicit none
    INTEGER, INTENT(IN) :: nrows1, ncols1, ncols2
    REAL, INTENT(IN) :: matrix1(nrows1, ncols1), matrix2(ncols1, ncols2)
    REAL, INTENT(OUT) :: result(nrows1, ncols2)
    INTEGER :: i, j, k
 
    ! Perform matrix multiplication
    do i = 1, nrows1
        do j = 1, ncols2
            result(i, j) = 0.0
            do k = 1, ncols1
                result(i, j) = result(i, j) + matrix1(i, k)*matrix2(k, j)
            end do
        end do
    end do
END SUBROUTINE mtx_mult