Commit c912910c authored by Florian Goth's avatar Florian Goth
Browse files

fix some debugging stuff that made it somehow out...

parent d590116f
Pipeline #12970 passed with stages
in 15 minutes and 15 seconds
...@@ -128,25 +128,19 @@ contains ...@@ -128,25 +128,19 @@ contains
use Operator_mod use Operator_mod
use graphdata_mod use graphdata_mod
implicit none implicit none
class(CmplxmscbOpT) :: this class(CmplxmscbOpT) :: this
Type(Operator), intent(in) :: Op_T type(Operator), intent(in) :: Op_T
integer, intent(in) :: method integer, intent(in) :: method
Integer :: i, k integer :: i
type(GraphData) :: gd type(GraphData) :: gd
Complex(kind=kind(0.D0)), allocatable, dimension(:,:) :: tmp Complex(kind=kind(0.D0)), allocatable, dimension(:,:) :: tmp
Real(kind=kind(0.D0)), allocatable, dimension(:) :: diags Real(kind=kind(0.D0)), allocatable, dimension(:) :: diags
write (*,*) "[CmplxMSCB]: init"
this%Zero = 1.E-12 this%Zero = 1.E-12
this%Ndim_hop = Op_T%N this%Ndim_hop = Op_T%N
this%g = Op_T%g this%g = Op_T%g
do i = 1, 20
Op_T%O(i, i+1) = 0
Op_T%O(i+1, i) = 0
enddo
call Op_T_to_graphdata(Op_T, gd, diags) call Op_T_to_graphdata(Op_T, gd, diags)
! some sanity checks and status informations ! some sanity checks and status informations
call determine_used_colors_of_graph(gd) call determine_used_colors_of_graph(gd)
...@@ -161,9 +155,9 @@ contains ...@@ -161,9 +155,9 @@ contains
! check wether it is supported behaviour ! check wether it is supported behaviour
do i = 1, size(Op_T%P) do i = 1, size(Op_T%P)
if (Op_T%P(i) /= i) then if (Op_T%P(i) /= i) then
write (*,*) "P unsupported" write (*,*) "P unsupported"
endif endif
enddo enddo
this%P = Op_T%P this%P = Op_T%P
call dealloc_graphdata(gd) call dealloc_graphdata(gd)
...@@ -178,9 +172,8 @@ contains ...@@ -178,9 +172,8 @@ contains
If ( dble(this%g*conjg(this%g)) > this%Zero ) then If ( dble(this%g*conjg(this%g)) > this%Zero ) then
call this%fe%adjoint_over_two(arg) call this%fe%adjoint_over_two(arg)
Endif Endif
end subroutine end subroutine
subroutine CmplxmscbOpT_rmult(this, arg) subroutine CmplxmscbOpT_rmult(this, arg)
class(CmplxmscbOpT), intent(in) :: this class(CmplxmscbOpT), intent(in) :: this
Complex(kind=kind(0.D0)), intent(inout), dimension(:,:) :: arg Complex(kind=kind(0.D0)), intent(inout), dimension(:,:) :: arg
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment