Fix poly%clone allocation of poly_beta

6 days ago · b7ac77da69
parent 8a66a7b214
commit b7ac77da69
2 changed files with 12 additions and 2 deletions
--- a/amgprec/impl/smoother/amg_d_poly_smoother_clone_settings.f90
+++ b/amgprec/impl/smoother/amg_d_poly_smoother_clone_settings.f90
@ -61,7 +61,12 @@ subroutine amg_d_poly_smoother_clone_settings(sm,smout,info)
    smout%rho_ba    = sm%rho_ba
    smout%rho_estimate  = sm%rho_estimate
    smout%rho_estimate_iterations  = sm%rho_estimate_iterations
-    smout%poly_beta = sm%poly_beta
+    if (allocated(sm%poly_beta)) then 
+      smout%poly_beta = sm%poly_beta
+    else
+      if (allocated(smout%poly_beta)) deallocate(smout%poly_beta)
+    end if
+
    
    if (allocated(smout%sv)) then
      if (.not.same_type_as(sm%sv,smout%sv)) then
--- a/amgprec/impl/smoother/amg_s_poly_smoother_clone_settings.f90
+++ b/amgprec/impl/smoother/amg_s_poly_smoother_clone_settings.f90
@ -61,7 +61,12 @@ subroutine amg_s_poly_smoother_clone_settings(sm,smout,info)
    smout%rho_ba    = sm%rho_ba
    smout%rho_estimate  = sm%rho_estimate
    smout%rho_estimate_iterations  = sm%rho_estimate_iterations
-    smout%poly_beta = sm%poly_beta
+    if (allocated(sm%poly_beta)) then 
+      smout%poly_beta = sm%poly_beta
+    else
+      if (allocated(smout%poly_beta)) deallocate(smout%poly_beta)
+    end if
+
    
    if (allocated(smout%sv)) then
      if (.not.same_type_as(sm%sv,smout%sv)) then