Source

PetIGA / src / petigainv.f90.in

Diff from to

src/petigainv.f90.in

             - A(2,1) * ( A(1,2)*A(3,3) - A(3,2)*A(1,3) ) &
             + A(3,1) * ( A(1,2)*A(2,3) - A(2,2)*A(1,3) )
   case default
-     detA = 0.0
+     detA = 0
   end select
 end function Determinant
 
   real   (kind=IGA_REAL_KIND   )  :: invA(dim,dim)
   select case (dim)
   case (1)
-     invA = 1.0/detA
+     invA = 1/detA
   case (2)
      invA(1,1) = + A(2,2)
      invA(2,1) = - A(2,1)
      invA(3,3) = + A(1,1)*A(2,2) - A(1,2)*A(2,1)
      invA = invA/detA
   case default
-     invA = 0.0
+     invA = 0
   end select
 end function Inverse