aselimov/CACmb
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Add check for orientation to make sure it's right handed and orthogonal

Browse Source
This commit is contained in:
Alex Selimov 2020-03-04 21:06:12 -05:00
parent 9038242b98
commit d2f3c9e3ec
1 changed files with 11 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
src/opt_orient.f90
View File

@ -20,6 +20,7 @@ module opt_orient
integer, intent(inout) :: arg_pos integer, intent(inout) :: arg_pos
integer :: i, ibasis, inod integer :: i, ibasis, inod
logical :: isortho, isrighthanded
real(kind=dp) :: inv_sub_box_ori(3,3,sub_box_num) real(kind=dp) :: inv_sub_box_ori(3,3,sub_box_num)
!First parse the orient command !First parse the orient command
@ -70,6 +71,7 @@ module opt_orient
integer, intent(inout) :: arg_pos integer, intent(inout) :: arg_pos
integer :: i, arg_len integer :: i, arg_len
logical :: isortho, isrighthanded
character(len=8) :: ori_string character(len=8) :: ori_string
!Pull out the new user orientation !Pull out the new user orientation
@ -82,6 +84,15 @@ module opt_orient
!Normalize the orientation matrix !Normalize the orientation matrix
new_orient = matrix_normal(new_orient,3) new_orient = matrix_normal(new_orient,3)
!Check right hand rule and orthogonality
call check_right_ortho(new_orient, isortho, isrighthanded)
if (.not.isortho) then
stop "Directions in orient are not orthogonal"
else if (.not.isrighthanded) then
stop "Directions in orient are not righthanded"
end if
arg_pos = arg_pos + 1 arg_pos = arg_pos + 1
end subroutine parse_orient end subroutine parse_orient