14 |
|
#ifdef IS_MPI |
15 |
|
real( kind = dp ), allocatable, dimension(:,:), public :: q_Row |
16 |
|
real( kind = dp ), allocatable, dimension(:,:), public :: q_Col |
17 |
+ |
real( kind = dp ), allocatable, dimension(:,:), public :: qcom_Row |
18 |
+ |
real( kind = dp ), allocatable, dimension(:,:), public :: qcom_Col |
19 |
|
real( kind = dp ), allocatable, dimension(:,:), public :: u_l_Row |
20 |
|
real( kind = dp ), allocatable, dimension(:,:), public :: u_l_Col |
21 |
|
real( kind = dp ), allocatable, dimension(:,:), public :: A_Row |
77 |
|
thisStat = -1 |
78 |
|
return |
79 |
|
endif |
80 |
+ |
|
81 |
+ |
allocate(qcom_Row(ndim,nrow),stat=alloc_stat) |
82 |
+ |
if (alloc_stat /= 0 ) then |
83 |
+ |
thisStat = -1 |
84 |
+ |
return |
85 |
+ |
endif |
86 |
+ |
|
87 |
+ |
allocate(qcom_Col(ndim,ncol),stat=alloc_stat) |
88 |
+ |
if (alloc_stat /= 0 ) then |
89 |
+ |
thisStat = -1 |
90 |
+ |
return |
91 |
+ |
endif |
92 |
|
|
93 |
|
allocate(u_l_Row(ndim,nrow),stat=alloc_stat) |
94 |
|
if (alloc_stat /= 0 ) then |
250 |
|
if (allocated(A_Row)) deallocate(A_Row) |
251 |
|
if (allocated(u_l_Col)) deallocate(u_l_Col) |
252 |
|
if (allocated(u_l_Row)) deallocate(u_l_Row) |
253 |
+ |
if (allocated(qcom_Col)) deallocate(qcom_Col) |
254 |
+ |
if (allocated(qcom_Row)) deallocate(qcom_Row) |
255 |
|
if (allocated(q_Col)) deallocate(q_Col) |
256 |
|
if (allocated(q_Row)) deallocate(q_Row) |
257 |
|
#else |