initialize_type_crs Subroutine

    subroutine initialize_type_crs(self, num_nodes, row, col)
        implicit none
        class(type_crs), intent(inout) :: self
        integer(int32), intent(in) :: num_nodes
        integer(int32), intent(in), optional :: row(:)
        integer(int32), intent(in), optional :: col(:)

        integer(int32) :: i

        if (.not. present(row) .or. .not. present(col)) then
            print *, "Error: row and col must be provided for CRS initialization."
            stop
        end if

        self%nnz = size(col)
        self%num_row = num_nodes
        self%num_ptr = self%num_row + 1

        call allocate_array(self%ptr, self%num_ptr)
        do i = 1, self%num_ptr
            self%ptr(i) = row(i)
        end do

        call allocate_array(self%ind, self%nnz)
        call allocate_array(self%val, self%nnz)
        do i = 1, self%nnz
            self%ind(i) = col(i)
            self%val(i) = 0.0d0
        end do

    end subroutine initialize_type_crs