scale_meshfield_analysis_numerror_base Module Reference

module FElib / Data / Statistics / numerical error More...

Data Types
type	meshfieldanalysisnumerrorbase

Functions/Subroutines
subroutine	meshfield_analysis_numerror_base_init (this, porder_error_check, ndim, np, intrp_np, log_fname_base, log_step_interval)
	Setup.

Detailed Description

module FElib / Data / Statistics / numerical error

Description

This module provides a base class useful for evaluating numerical errrors

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

meshfield_analysis_numerror_base_init()

subroutine scale_meshfield_analysis_numerror_base::meshfield_analysis_numerror_base_init	(	class(meshfieldanalysisnumerrorbase), intent(inout)	this,
		integer, intent(in)	porder_error_check,
		integer, intent(in)	ndim,
		integer, intent(in)	np,
		integer, intent(in)	intrp_np,
		character(len=*), intent(in)	log_fname_base,
		integer, intent(in)	log_step_interval )

Setup.

Definition at line 86 of file scale_meshfield_analysis_numerror_base.F90.

    implicit none
    class(MeshFieldAnalysisNumerrorBase), intent(inout) :: this
    integer, intent(in) :: porder_error_check
    integer, intent(in) :: ndim
    integer, intent(in) :: np
    integer, intent(in) :: intrp_np
    character(len=*), intent(in) :: log_fname_base
    integer, intent(in) :: log_step_interval
 
    character(len=H_MID) :: fname
    integer :: ierr
    !---------------------------------------------------------------------------
 
    this%var_num = 0
 
    this%log_step_interval = log_step_interval
    this%log_rstep         = log_step_interval
 
    this%output_error_first = .true. 
 
    !--
    if ( prc_ismaster ) then
      this%log_fid = io_get_available_fid()
      fname = trim(log_fname_base)//".peall"
      open( unit   = this%log_fid,                 &
        file   = fname,                            &
        form   = 'formatted',                      &
        iostat = ierr                              )
      if ( ierr /= 0 ) then
        log_error('MeshField_NumErrorAnalysis_Init',*) 'File open error! :', trim(fname)
        call prc_abort
      endif 
    end if
 
    !--
    this%PolyOrderErrorCheck = porder_error_check
    this%intrp_np = intrp_np
    this%ndim = ndim
    allocate( this%IntrpMat(this%intrp_np,np) )
    allocate( this%intw_intrp(this%intrp_np) )
    allocate( this%epos_intrp(this%intrp_np,this%ndim) )
 
    return