Skip to content

Commit 11860a6

Browse files
committed
issues in vonMises addressed
1 parent f2a9e0f commit 11860a6

File tree

3 files changed

+10
-10
lines changed

3 files changed

+10
-10
lines changed

src/common/m_variables_conversion.fpp

Lines changed: 5 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -1032,12 +1032,12 @@ contains
10321032
qK_prim_vf(i)%sf(j, k, l) = qK_cons_vf(i)%sf(j, k, l)/rho_K
10331033
end do
10341034
end if
1035-
#ifndef POST_PROCESS
1035+
#ifndef MFC_POST_PROCESS
10361036
! to save von Mises stress instead of elastic internal energy
1037-
qK_prim_vf(xiend+1)%sf(j, k, l) = sqrt((3/2)*qK_prim_vf(1)%sf(j, k, l)**2 + &
1038-
2*qK_prim_vf(2)%sf(j, k, l)**2 + qK_prim_vf(3)%sf(j, k, l)**2 + &
1039-
2*qK_prim_vf(4)%sf(j, k, l)**2 + 2*qK_prim_vf(5)%sf(j, k, l)**2 + &
1040-
qK_prim_vf(6)%sf(j, k, l)**2)
1037+
qK_prim_vf(xiend+1)%sf(j, k, l) = sqrt((3d0/2d0)*(qK_prim_vf(strxb)%sf(j, k, l)**2d0 + &
1038+
2d0*qK_prim_vf(strxb+1)%sf(j, k, l)**2d0 + qK_prim_vf(strxb+2)%sf(j, k, l)**2d0 + &
1039+
2d0*qK_prim_vf(strxb+3)%sf(j, k, l)**2d0 + 2d0*qK_prim_vf(strxb+4)%sf(j, k, l)**2d0 + &
1040+
qK_prim_vf(strxe)%sf(j, k, l)**2d0))
10411041
#endif
10421042

10431043
!$acc loop seq

src/post_process/m_start_up.f90

Lines changed: 4 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -394,7 +394,10 @@ subroutine s_save_data(t_step, varname, pres, c, H)
394394
call s_write_variable_to_formatted_database_file(varname, t_step)
395395
end if
396396
if (prim_vars_wrt) then
397-
q_sf = q_prim_vf(xiend+1)%sf(j, k, l)
397+
q_sf = q_prim_vf(xiend+1)%sf( &
398+
-offset_x%beg:m + offset_x%end, &
399+
-offset_y%beg:n + offset_y%end, &
400+
-offset_z%beg:p + offset_z%end)
398401
write (varname, '(A,I0)') 'vonMises'
399402
call s_write_variable_to_formatted_database_file(varname, t_step)
400403
end if

submit_hyper.sh

Lines changed: 1 addition & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -1,7 +1,4 @@
11
#!/bin/bash
22

3-
./mfc.sh run ./examples/3D_hyperelasticity/hyper_gel.py -p batch -N 1 -n 4 -g 0 -w 01:00:00 -# test1 -t pre_process -c oscar
4-
./mfc.sh run ./examples/3D_hyperelasticity/hyper_gel.py -p batch -N 1 -n 4 -g 0 -w 01:00:00 -# test1 -t simulation -c oscar
53
#./mfc.sh run /scratch/bciv/rodrigu1/cav2024/hyper/hyper_gel.py -p gpuA40x4 -N 1 -n 1 -g 1 -w 01:00:00 -# test1 -t post_process -a bciv-delta-gpu -c delta
6-
7-
4+
./mfc.sh run /projects/bciv/mcarcanabarbosa/kymo/geldef.py -p gpuA100x4 -N 1 -n 4 -g 1 -w 01:00:00 -# test1 -t post_process -a bciv-delta-gpu -c delta

0 commit comments

Comments
 (0)