fft-split

run.sh
View Options

			#
			# @file run.sh
			#
			# @brief
			#
			# @copyright
			# Copyright (©) 2021 EPFL (Ecole Polytechnique Fédérale de Lausanne)
			# SPC (Swiss Plasma Center)
			#
			# spclibs is free software: you can redistribute it and/or modify it under
			# the terms of the GNU Lesser General Public License as published by the Free
			# Software Foundation, either version 3 of the License, or (at your option)
			# any later version.
			#
			# spclibs is distributed in the hope that it will be useful, but WITHOUT ANY
			# WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
			# FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
			#
			# You should have received a copy of the GNU Lesser General Public License
			# along with this program. If not, see <https://www.gnu.org/licenses/>.
			#
			# @authors
			# (in alphabetical order)
			# @author Trach-Minh Tran <trach-minh.tran@epfl.ch>
			#
			#!/bin/bash

			EXE=/home/ttran/bsplines/multigrid/src/poisson_mg
			TMP=/misc/multigrid
			[ -e $TMP ] \|\| mkdir -p $TMP

			cat > in0 <<EOF
			&parameters
			prb='dddd'
			mat_type='cds'
			nx=xxx, ny=yyy
			nz=1,
			kx=4, ky=4,
			Lx=100.d00, Ly=800.d00,
			nldebug=f
			nldirect=f,
			icrosst=16,
			beta=-1E-2
			levels=12
			nnu=1
			nu1= 3,4,5
			nu2= 3,4,5
			mu=1, nu0=1,
			relax= 'gs', omega=0.9
			nits=15, atol=0., rtol=1.e-8
			/
			EOF
			for x in 32 64 128 256 512; do
			y=$(echo "4 * $x" \| bc)
			echo -n "Nx = $x, Ny = $y: "
			sed "s/xxx/$x/g" in0 \| sed "s/yyy/$y/g" \| ${EXE} \| grep ' 3 3 '
			mv poisson_mg.h5 $TMP/dddd_${x}x${y}.h5
			done