scale_element_modalfilter Module Reference

module FElib / Element/ ModalFilter More...

Data Types
type	modalfilter
	Derived type representing a modal filter. More...

Functions/Subroutines
subroutine	modalfilter_init_line (this, elem, etac, alpha, ord, tend_flag)

Detailed Description

module FElib / Element/ ModalFilter

Description

A module for modal filtering

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

modalfilter_init_line()

subroutine scale_element_modalfilter::modalfilter_init_line	(	class(modalfilter), intent(inout)	this,
		class(lineelement), intent(in)	elem,
		real(rp), intent(in)	etac,
		real(rp), intent(in)	alpha,
		integer, intent(in)	ord,
		logical, intent(in), optional	tend_flag )

Definition at line 61 of file scale_element_modalfilter.F90.

 
    implicit none
    class(ModalFilter), intent(inout) :: this
    class(LineElement), intent(in) :: elem
    real(RP), intent(in) :: etac
    real(RP), intent(in) :: alpha
    integer, intent(in) :: ord
    logical, intent(in), optional :: tend_flag
 
    real(RP) :: filter1D(elem%Np)
    integer :: p
    logical :: tend_flag_
    !----------------------------------------------------
    
    tend_flag_ = .false.
    if ( present(tend_flag) ) tend_flag_ =  tend_flag
 
    call get_exp_filter( filter1d,               & ! (out)
      etac, alpha, ord, elem%Np, elem%PolyOrder, & ! (in)
      tend_flag_ )                                 ! (in)
 
    allocate( this%FilterMat(elem%Np,elem%Np) )
    this%FilterMat(:,:) = 0.0_rp
    do p=1, elem%Np
      this%FilterMat(p,p) = filter1d(p)
    end do
    this%FilterMat(:,:) = matmul(this%FilterMat, elem%invV)
    this%FilterMat(:,:) = matmul(elem%V, this%FilterMat)
    !$acc enter data copyin(this%FilterMat)
    
    return