10
1
mirror of https://gitlab.com/scemama/QCaml.git synced 2024-07-25 04:07:24 +02:00

AngularMomentum has + and -

This commit is contained in:
Anthony Scemama 2018-03-14 21:58:55 +01:00
parent e5da95100a
commit 99742b85ac
4 changed files with 15 additions and 3 deletions

View File

@ -87,8 +87,8 @@ let create ?cutoff p_a p_b =
Cs.(center p_a = center p_b)
in
let totAngMomInt =
Am.to_int (Cs.totAngMom p_a) +
Am.to_int (Cs.totAngMom p_b)
Am.(Cs.totAngMom p_a + Cs.totAngMom p_b)
|> Am.to_int
in
Some {
Sp.i ; j ;

View File

@ -101,7 +101,8 @@ let contracted_class_shell_pair ~zero_m shell_p geometry : float Zmap.t =
and shell_b = Csp.shell_b shell_p
in
let maxm =
Am.to_int (Cs.totAngMom shell_a) + Am.to_int (Cs.totAngMom shell_b)
Am.(Cs.totAngMom shell_a + Cs.totAngMom shell_b)
|> Am.to_int
in
(* Pre-computation of integral class indices *)

View File

@ -144,3 +144,10 @@ let zkey_array a =
Hashtbl.add zkey_array_memo a result;
result
let ( + ) a b =
of_int ( (to_int a) + (to_int b) )
let ( - ) a b =
of_int ( (to_int a) - (to_int b) )

View File

@ -85,3 +85,7 @@ val zkey_array : kind -> Zkey.t array
]}
*)
val ( + ) : t -> t -> t
val ( - ) : t -> t -> t