def PrintDiffusiveFluxMatrix(G, params):
    """Auxiliary function to print the diffusive flux matrix (G)"""

    dim = params.dim
    print("The diffusive matrix is:\n")
    for ll in range(dim+2):
        for mm in range(dim):
            print("G[", ll, ",", mm, "]=", G[ll, mm], "\n")

    return 0