Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Interpolation for h-refined meshes #154

Draft
wants to merge 3 commits into
base: master
Choose a base branch
from
Draft

Interpolation for h-refined meshes #154

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
68d465d
Initial submit for retaining dofs on h-refined meshes.
stefanhenneking Oct 8, 2023
ef9287d
Routine for initiating dofs on new nodes.
stefanhenneking Oct 8, 2023
9c6c539
updating routines to be compatible with most recent hp3D version.
stefanhenneking Oct 8, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
18 changes: 17 additions & 1 deletion trunk/src/meshmod/break.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,3 +1,4 @@
#include "typedefs.h"
!-------------------------------------------------------------------------
!> @brief routine breaks an element middle node according to a given
!! refinement flag. Compatibility is assumed across the faces
Expand All @@ -6,17 +7,24 @@
!> @param[in] Mdle - middle node
!> @param[in] Kref - refinement flag
!!
!> @date Feb 2023
!> @date Oct 2023
!-------------------------------------------------------------------------
subroutine break(Mdle,Kref)
!
use data_structure3D
use element_data
use refinements
use refinements_history
!
implicit none
!
integer, intent(in) :: Mdle, Kref
!
real(8) :: xnod (NDIMEN ,MAXbrickH)
VTYPE :: zdofH(MAXEQNH,MAXbrickH)
VTYPE :: zdofE(MAXEQNE,MAXbrickE)
VTYPE :: zdofV(MAXEQNV,MAXbrickV)
VTYPE :: zdofQ(MAXEQNQ,MAXbrickQ)
!
integer :: ntype
integer :: nodesl(27),norientl(27)
Expand All @@ -30,6 +38,11 @@ subroutine break(Mdle,Kref)
!
iprint = 0
!
!..save element dof
call nodcor(Mdle, xnod)
call solelm(Mdle, zdofH,zdofE,zdofV,zdofQ)
call save_element(Mdle,xnod,zdofH,zdofE,zdofV,zdofQ)
!
!..nodal connectivities
call elem_nodes(Mdle, nodesl,norientl)
!
Expand Down Expand Up @@ -153,6 +166,9 @@ subroutine break(Mdle,Kref)
call nr_mdle_sons(ntype,Kref, nrsons)
NRELES=NRELES+nrsons-1
!
!..generate new dof for the new and active nodes
!!!!! call initiate_dof(Mdle,Xnod,ZdofH,ZdofE,ZdofV,ZdofQ)
!
#if DEBUG_MODE
if (iprint.ge.1) then
write(*,7100) Mdle, S_Type(NODES(Mdle)%ntype), Kref
Expand Down
Loading