Skip to content

Commit

Permalink
Remove double operators (+ +) and (+ -)
Browse files Browse the repository at this point in the history 
These cause warnings with the gnu compiler, and are unnecessary (tested
with gfortran 5.3.0). For example:

inv.f90:292:14:
              - ( amat(3,1) * bmat(1,3) &
              1
Warning: Extension: Unary operator following arithmetic operator (use parentheses) at (1)
  • Loading branch information
@mark-petersen
mark-petersen committed Jul 2, 2019
1 parent 3565130 commit 4e62647
Showing 1 changed file with 18 additions and 18 deletions.
36 changes: 18 additions & 18 deletions src/inv.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -212,83 +212,83 @@ pure subroutine mul_3x3(amat,adim,bmat,bdim, &

if (scal .eq. +1.d0) then

cmat(1,1) = cmat(1,1) + &
cmat(1,1) = cmat(1,1) &
+ ( amat(1,1) * bmat(1,1) &
+ amat(1,2) * bmat(2,1) &
+ amat(1,3) * bmat(3,1) )
cmat(2,1) = cmat(2,1) + &
cmat(2,1) = cmat(2,1) &
+ ( amat(2,1) * bmat(1,1) &
+ amat(2,2) * bmat(2,1) &
+ amat(2,3) * bmat(3,1) )
cmat(3,1) = cmat(3,1) + &
cmat(3,1) = cmat(3,1) &
+ ( amat(3,1) * bmat(1,1) &
+ amat(3,2) * bmat(2,1) &
+ amat(3,3) * bmat(3,1) )

cmat(1,2) = cmat(1,2) + &
cmat(1,2) = cmat(1,2) &
+ ( amat(1,1) * bmat(1,2) &
+ amat(1,2) * bmat(2,2) &
+ amat(1,3) * bmat(3,2) )
cmat(2,2) = cmat(2,2) + &
cmat(2,2) = cmat(2,2) &
+ ( amat(2,1) * bmat(1,2) &
+ amat(2,2) * bmat(2,2) &
+ amat(2,3) * bmat(3,2) )
cmat(3,2) = cmat(3,2) + &
cmat(3,2) = cmat(3,2) &
+ ( amat(3,1) * bmat(1,2) &
+ amat(3,2) * bmat(2,2) &
+ amat(3,3) * bmat(3,2) )

cmat(1,3) = cmat(1,3) + &
cmat(1,3) = cmat(1,3) &
+ ( amat(1,1) * bmat(1,3) &
+ amat(1,2) * bmat(2,3) &
+ amat(1,3) * bmat(3,3) )
cmat(2,3) = cmat(2,3) + &
cmat(2,3) = cmat(2,3) &
+ ( amat(2,1) * bmat(1,3) &
+ amat(2,2) * bmat(2,3) &
+ amat(2,3) * bmat(3,3) )
cmat(3,3) = cmat(3,3) + &
cmat(3,3) = cmat(3,3) &
+ ( amat(3,1) * bmat(1,3) &
+ amat(3,2) * bmat(2,3) &
+ amat(3,3) * bmat(3,3) )

else &
if (scal .eq. -1.d0) then

cmat(1,1) = cmat(1,1) + &
cmat(1,1) = cmat(1,1) &
- ( amat(1,1) * bmat(1,1) &
+ amat(1,2) * bmat(2,1) &
+ amat(1,3) * bmat(3,1) )
cmat(2,1) = cmat(2,1) + &
cmat(2,1) = cmat(2,1) &
- ( amat(2,1) * bmat(1,1) &
+ amat(2,2) * bmat(2,1) &
+ amat(2,3) * bmat(3,1) )
cmat(3,1) = cmat(3,1) + &
cmat(3,1) = cmat(3,1) &
- ( amat(3,1) * bmat(1,1) &
+ amat(3,2) * bmat(2,1) &
+ amat(3,3) * bmat(3,1) )

cmat(1,2) = cmat(1,2) + &
cmat(1,2) = cmat(1,2) &
- ( amat(1,1) * bmat(1,2) &
+ amat(1,2) * bmat(2,2) &
+ amat(1,3) * bmat(3,2) )
cmat(2,2) = cmat(2,2) + &
cmat(2,2) = cmat(2,2) &
- ( amat(2,1) * bmat(1,2) &
+ amat(2,2) * bmat(2,2) &
+ amat(2,3) * bmat(3,2) )
cmat(3,2) = cmat(3,2) + &
cmat(3,2) = cmat(3,2) &
- ( amat(3,1) * bmat(1,2) &
+ amat(3,2) * bmat(2,2) &
+ amat(3,3) * bmat(3,2) )

cmat(1,3) = cmat(1,3) + &
cmat(1,3) = cmat(1,3) &
- ( amat(1,1) * bmat(1,3) &
+ amat(1,2) * bmat(2,3) &
+ amat(1,3) * bmat(3,3) )
cmat(2,3) = cmat(2,3) + &
cmat(2,3) = cmat(2,3) &
- ( amat(2,1) * bmat(1,3) &
+ amat(2,2) * bmat(2,3) &
+ amat(2,3) * bmat(3,3) )
cmat(3,3) = cmat(3,3) + &
cmat(3,3) = cmat(3,3) &
- ( amat(3,1) * bmat(1,3) &
+ amat(3,2) * bmat(2,3) &
+ amat(3,3) * bmat(3,3) )
Expand Down

0 comments on commit 4e62647

Please sign in to comment.