1. Lisandro Dalcin
  2. PetIGA

Source

PetIGA / src / petigaval.F90

Diff from to

File src/petigaval.F90

  • Ignore whitespace
      end select
   case (1)
      select case (axis)
-     case (0); N = (/ 1.0 /)
+     case (0); N = 1
      end select
   end select
   dS = sqrt(sum(N*N))
   real(kind=IGA_REAL_KIND)             :: n(3)
   real(kind=IGA_REAL_KIND), intent(in) :: s(3)
   real(kind=IGA_REAL_KIND), intent(in) :: t(3)
-  ! c_i = eps_ijk s_j tj
+  ! n_i = eps_ijk s_j tj
   n(1) = s(2) * t(3) - s(3) * t(2)
   n(2) = s(3) * t(1) - s(1) * t(3)
   n(3) = s(1) * t(2) - s(2) * t(1)
   scalar (kind=IGA_SCALAR_KIND ), intent(out)      :: V(dof)
   integer(kind=IGA_INTEGER_KIND)  :: a, i
   ! V = matmul(N,transpose(U))
-  V = 0.0
+  V = 0
   do a = 1, nen
      V = V + N(a) * U(:,a)
   end do
   scalar (kind=IGA_SCALAR_KIND ), intent(out)      :: V(dim,dof)
   integer(kind=IGA_INTEGER_KIND)  :: a, c
   ! V = matmul(N,transpose(U))
-  V = 0.0
+  V = 0
   do a = 1, nen
      do c = 1, dof
         V(:,c) = V(:,c) + N(:,a) * U(c,a)
   scalar (kind=IGA_SCALAR_KIND ), intent(out)      :: V(dim*dim,dof)
   integer(kind=IGA_INTEGER_KIND)  :: a, i
   ! V = matmul(N,transpose(U))
-  V = 0.0
+  V = 0
   do a = 1, nen
      do i = 1, dof
         V(:,i) = V(:,i) + N(:,a) * U(i,a)
   scalar (kind=IGA_SCALAR_KIND ), intent(in)       :: U(dof,nen)
   scalar (kind=IGA_SCALAR_KIND ), intent(out)      :: V(dof)
   integer(kind=IGA_INTEGER_KIND)  :: a, c, i
-  V = 0.0
+  V = 0
   do a = 1, nen
      do c = 1, dof
         do i = 1, dim
   scalar (kind=IGA_SCALAR_KIND ), intent(out)      :: V(dim*dim*dim,dof)
   integer(kind=IGA_INTEGER_KIND)  :: a, i
   ! V = matmul(N,transpose(U))
-  V = 0.0
+  V = 0
   do a = 1, nen
      do i = 1, dof
         V(:,i) = V(:,i) + N(:,a) * U(i,a)
 !  scalar (kind=IGA_SCALAR_KIND ), intent(out)      :: V(dim**der,dof)
 !  integer(kind=IGA_INTEGER_KIND)  :: a, i
 !  ! V = matmul(N,transpose(U))
-!  V = 0.0
+!  V = 0
 !  do a = 1, nen
 !     do i = 1, dof
 !        V(:,i) = V(:,i) + N(:,a) * U(i,a)
   scalar (kind=IGA_SCALAR_KIND ), intent(out)      :: V(dim**der,dof)
   integer(kind=IGA_INTEGER_KIND)  :: a, i
   ! V = matmul(N,transpose(U))
-  V = 0.0
+  V = 0
   do a = 1, nen
      do i = 1, dof
         V(:,i) = V(:,i) + N(:,a) * U(i,a)