def METHOD_NAME( Et, Et_pml, Ez, Bt, Bt_pml, Bz,
                  damp_array, n_pml ) :
    """
    Damp the E and B fields in the PML cells (i.e. the last n_pml cells
    in r), in an anisotropic manner which is given by the PML principles
    Parameters :
    ------------
    Et, Et_pml, Ez, Bt, Bt_pml, Bz : 2darrays of complexs
        Contain the fields to be damped
        The first axis corresponds to z and the second to r
    damp_array: 1darray of floats
        An array of length n_guards, which contains the damping factors
    n_pml: int
        Number of PML cells
    """
    # Obtain Cuda grid
    iz, i_pml = cuda.grid(2)
    # Obtain the size of the array along z and r
    Nz, Nr = Et.shape
    # Modify the fields
    if i_pml < n_pml:
        # Apply the damping arrays
        if iz < Nz:
            # Get the damping factor
            damp_factor= damp_array[i_pml]
            # Get the index in the bigger field array
            ir = Nr - n_pml + i_pml
            # Substract the theta PML fields to the regular theta fields
            Et[iz,ir] -= Et_pml[iz,ir]
            Bt[iz,ir] -= Bt_pml[iz,ir]
            # Damp the theta PML fields
            Et_pml[iz,ir] *= damp_factor
            Bt_pml[iz,ir] *= damp_factor
            # Add the theta PML fields back to the regular theta fields
            Et[iz,ir] += Et_pml[iz,ir]
            Bt[iz,ir] += Bt_pml[iz,ir]
            # Damp the z fields
            Ez[iz,ir] *= damp_factor
            Bz[iz,ir] *= damp_factor