cleaning

src/crystallite.f90

@@ -22,7 +22,6 @@ module crystallite
   use constitutive
   use discretization
   use lattice
-  use future
   use plastic_nonlocal
   use geometry_plastic_nonlocal, only: &
     nIPneighbors    => geometry_plastic_nonlocal_nIPneighbors, &

src/lattice.f90

@@ -11,7 +11,6 @@ module lattice
   use IO
   use config
   use math
-  use future
  
   implicit none
   private
@@ -493,11 +492,6 @@ module lattice
   integer(kind(LATTICE_undefined_ID)),        dimension(:),       allocatable, public, protected :: &
     lattice_structure
  
-  
-  interface lattice_forestProjection ! DEPRECATED, use lattice_forestProjection_edge
-    module procedure slipProjection_transverse
-  end interface lattice_forestProjection
-  
   interface lattice_forestProjection_edge
     module procedure slipProjection_transverse
   end interface lattice_forestProjection_edge
@@ -529,7 +523,6 @@ module lattice
    lattice_characteristicShear_Twin, &
    lattice_C66_twin, &
    lattice_C66_trans, &
-   lattice_forestProjection, &
    lattice_forestProjection_edge, &
    lattice_forestProjection_screw, &
    lattice_slip_normal, &

src/material.f90

@@ -474,9 +474,9 @@ subroutine material_parseHomogenization
   enddo
  
   do h=1, size(config_homogenization)
-    homogenization_typeInstance(h)  = count(homogenization_type(1:h)  == homogenization_type(h))
+    homogenization_typeInstance(h)  = count(homogenization_type(1:h) == homogenization_type(h))
-    thermal_typeInstance(h)         = count(thermal_type       (1:h)  == thermal_type       (h))
+    thermal_typeInstance(h)         = count(thermal_type       (1:h) == thermal_type       (h))
-    damage_typeInstance(h)          = count(damage_type        (1:h)  == damage_type        (h))
+    damage_typeInstance(h)          = count(damage_type        (1:h) == damage_type        (h))
   enddo
  
   homogenization_maxNgrains = maxval(homogenization_Ngrains,homogenization_active)
@@ -496,9 +496,9 @@ subroutine material_parseMicrostructure
   character(len=65536) :: &
     tag
  
-  allocate(microstructure_crystallite(size(config_microstructure)),          source=0)
+  allocate(microstructure_crystallite(size(config_microstructure)),   source=0)
-  allocate(microstructure_Nconstituents(size(config_microstructure)),        source=0)
+  allocate(microstructure_Nconstituents(size(config_microstructure)), source=0)
-  allocate(microstructure_active(size(config_microstructure)),               source=.false.)
+  allocate(microstructure_active(size(config_microstructure)),        source=.false.)
  
   if(any(discretization_microstructureAt > size(config_microstructure))) &
    call IO_error(155,ext_msg='More microstructures in geometry than sections in material.config')
@@ -531,7 +531,7 @@ subroutine material_parseMicrostructure
            case('texture')
              microstructure_texture(c,m) =  IO_intValue(strings(c),chunkPos,i+1)
            case('fraction')
-             microstructure_fraction(c,m) =  IO_floatValue(strings(c),chunkPos,i+1)
+             microstructure_fraction(c,m) = IO_floatValue(strings(c),chunkPos,i+1)
          end select
       
       enddo

src/math.f90

@@ -7,7 +7,6 @@
 !--------------------------------------------------------------------------------------------------
 module math
   use prec
-  use future
   use IO
   use debug
   use numerics
@@ -16,8 +15,7 @@ module math
 #if __INTEL_COMPILER >= 1900
   ! do not make use associated entities available to other modules
   private :: &
-    prec, &
+    prec
-    future
 #endif
 
   real(pReal),    parameter :: PI = acos(-1.0_pReal)                                                !< ratio of a circle's circumference to its diameter
@@ -75,15 +73,7 @@ module math
       3,3  &
       ],[2,9])                                                                                      !< arrangement in Plain notation
      
-!--------------------------------------------------------------------------------------------------
-! Provide deprecated name for compatibility
-  interface math_mul3x3
-    module procedure math_inner
-  end interface math_mul3x3
-  public :: &
-    math_mul3x3
 !---------------------------------------------------------------------------------------------------
-
  private :: &
    unitTest
 

src/mesh_abaqus.f90

@@ -10,6 +10,7 @@ module mesh
  use mesh_base
  use geometry_plastic_nonlocal
  use discretization
+ use math
 
  implicit none
  private
@@ -1758,9 +1759,6 @@ end subroutine mesh_build_sharedElems
 !> @brief build up of IP neighborhood, allocate globals '_ipNeighborhood'
 !--------------------------------------------------------------------------------------------------
 subroutine mesh_build_ipNeighborhood
- use math, only: &
-   math_mul3x3
-
  
  integer      ::           myElem, &                                                           ! my CP element index
                                  myIP, &
@@ -1908,7 +1906,7 @@ subroutine mesh_build_ipNeighborhood
          do pointingToMe = 1,FE_NipNeighbors(FE_celltype(neighboringType))                      ! find neighboring index that points from my neighbor to myself
            if (    myElem == mesh_ipNeighborhood(1,pointingToMe,neighboringIP,neighboringElem) &
                .and. myIP == mesh_ipNeighborhood(2,pointingToMe,neighboringIP,neighboringElem)) then ! possible candidate
-             if (math_mul3x3(mesh_ipAreaNormal(1:3,neighbor,myIP,myElem),&
+             if (math_inner(mesh_ipAreaNormal(1:3,neighbor,myIP,myElem),&
                              mesh_ipAreaNormal(1:3,pointingToMe,neighboringIP,neighboringElem)) < 0.0_pReal) then ! area normals have opposite orientation (we have to check that because of special case for single element with two ips and periodicity. In this case the neighbor is identical in two different directions.)
                mesh_ipNeighborhood(3,neighbor,myIP,myElem) = pointingToMe                            ! found match
                exit                                                                                  ! so no need to search further

src/mesh_base.f90

@@ -11,7 +11,6 @@ module mesh_base
   use, intrinsic :: iso_c_binding
   use prec
   use element
-  use future 
 
   implicit none
 

src/quaternions.f90

@@ -35,7 +35,6 @@
 !---------------------------------------------------------------------------------------------------
 module quaternions
   use prec
-  use future
 
   implicit none
   public