extend case of exponentiation, but make failure be a pattern match is…

…sues for now
camfort · Sep 16, 2024 · f0c07b3 · f0c07b3
1 parent d146080
commit f0c07b3
Showing 1 changed file with 6 additions and 4 deletions.
diff --git a/src/Language/Fortran/Repr/Eval/Value.hs b/src/Language/Fortran/Repr/Eval/Value.hs
@@ -273,10 +273,12 @@ evalBOp bop l r = do
           (FSVInt li, FSVInt ri) ->
             pure $ MkFScalarValue $ FSVInt $ fIntBOpInplace (^) li ri
           (FSVReal lr, FSVReal ri) ->
-            pure $ MkFScalarValue $ FSVReal $ fRealBOpInplace' exp exp lr ri
-              where
-                exp x y = x ** y
-          _ -> err $ ELazy "exponentiation: unsupported types"
+            pure $ MkFScalarValue $ FSVReal $ fRealBOpInplace' (**) (**) lr ri
+          (FSVReal lr, FSVInt ri) ->
+            -- Handle case of a real raised to an integer power.
+            pure $ MkFScalarValue $ FSVReal $ fRealBOpInplace' (**) (**) lr (FReal8 $ withFInt ri)
+
+--          _ -> err $ ELazy "exponentiation: unsupported types"
 
       F.Concatenation  ->
         case (l', r') of