347 lines
No EOL
39 KiB
HTML
347 lines
No EOL
39 KiB
HTML
<!-- -*- mode: jinja2 -*- -->
|
|
|
|
<!DOCTYPE html>
|
|
<html lang="en">
|
|
<head>
|
|
<meta charset="utf-8">
|
|
<meta http-equiv="X-UA-Compatible" content="IE=edge">
|
|
<meta name="viewport" content="width=device-width, initial-scale=1">
|
|
<meta name="description" content="Tridiagonal compact finite difference solver, mathematical operators, and MPI-IO drivers.">
|
|
<meta name="author" content="Google DeepMind Team & Ignis" >
|
|
<link rel="icon" href="../favicon.png">
|
|
|
|
<title>test_compact.f90 – HPC DNS Post-Processing Fortran Core</title>
|
|
|
|
<!-- Bootstrap -->
|
|
<link href="https://cdn.jsdelivr.net/npm/bootstrap@5.3.3/dist/css/bootstrap.min.css" rel="stylesheet"
|
|
integrity="sha384-QWTKZyjpPEjISv5WaRU9OFeRpok6YctnYmDr5pNlyT2bRjXh0JMhjY6hW+ALEwIH" crossorigin="anonymous">
|
|
<script src="https://cdn.jsdelivr.net/npm/bootstrap@5.3.3/dist/js/bootstrap.bundle.min.js"
|
|
integrity="sha384-YvpcrYf0tY3lHB60NNkmXc5s9fDVZLESaAA55NDzOxhy9GkcIdslK1eN7N6jIeHz" crossorigin="anonymous"></script>
|
|
<!-- Font Awesome -->
|
|
<link href="../css/fontawesome.min.css" rel="stylesheet">
|
|
<link href="../css/brands.min.css" rel="stylesheet">
|
|
<link href="../css/regular.min.css" rel="stylesheet">
|
|
<link href="../css/solid.min.css" rel="stylesheet">
|
|
<link href="../css/v4-font-face.min.css" rel="stylesheet">
|
|
<link href="../css/v4-shims.min.css" rel="stylesheet">
|
|
<!-- MathJax -->
|
|
<script type="text/x-mathjax-config">
|
|
MathJax.Hub.Config({
|
|
TeX: { equationNumbers: { autoNumber: "AMS" } }
|
|
});
|
|
</script>
|
|
<script src="https://cdn.jsdelivr.net/npm/mathjax@2.7.9/MathJax.js?config=TeX-AMS-MML_HTMLorMML" async
|
|
integrity="sha256-DViIOMYdwlM/axqoGDPeUyf0urLoHMN4QACBKyB58Uw=" crossorigin="anonymous"></script>
|
|
<!-- Other scripts and stylesheets -->
|
|
<link href="../css/local.css" rel="stylesheet">
|
|
<link href="../css/pygments.css" rel="stylesheet">
|
|
<script src="../js/svg-pan-zoom.min.js"></script>
|
|
</head>
|
|
|
|
<body>
|
|
|
|
<!-- Fixed navbar -->
|
|
<div class="container-fluid mb-sm-4 mb-xl-2">
|
|
<nav class="navbar navbar-expand-lg navbar-dark bg-dark fixed-top">
|
|
<div class="container">
|
|
<a class="navbar-brand" href="../index.html">HPC DNS Post-Processing Fortran Core <small>1.0.0</small></a>
|
|
<button type="button" class="navbar-toggler" data-bs-toggle="collapse" data-bs-target="#navbar"
|
|
aria-expanded="false" aria-controls="navbar" aria-label="Toggle navigation">
|
|
<span class="navbar-toggler-icon">
|
|
</button>
|
|
|
|
<div id="navbar" class="navbar-collapse collapse">
|
|
<ul class="navbar-nav">
|
|
<li class="nav-item">
|
|
<a class="nav-link" href="../lists/files.html">Source Files</a>
|
|
</li>
|
|
<li class="nav-item">
|
|
<a class="nav-link" href="../lists/modules.html">Modules</a>
|
|
</li>
|
|
<li class="nav-item">
|
|
<a class="nav-link" href="../lists/procedures.html">Procedures</a>
|
|
</li>
|
|
<li class="nav-item">
|
|
<a class="nav-link" href="../lists/programs.html">Programs</a>
|
|
</li>
|
|
</ul>
|
|
<div class="d-flex align-items-end flex-grow-1">
|
|
<form action="../search.html" role="search" class="ms-auto">
|
|
<input type="text" class="form-control" aria-label="Search" placeholder="Search" name="q" id="tipue_search_input" autocomplete="off" required>
|
|
</form>
|
|
</div>
|
|
</div><!--/.nav-collapse -->
|
|
</div>
|
|
</nav>
|
|
</div>
|
|
|
|
<div class="container">
|
|
<div class="row">
|
|
<h1>test_compact.f90
|
|
<small>Source File</small>
|
|
|
|
</h1>
|
|
<div class="container p-2 mb-4 bg-light border rounded-3">
|
|
<div class="row align-items-center justify-content-between" id="info-bar">
|
|
<div class="col">
|
|
<ul class="list-inline" style="margin-bottom:0px;display:inline">
|
|
|
|
<li class="list-inline-item" id="statements"><i class="fa fa-list-ol"></i>
|
|
<a data-bs-toggle="tooltip" data-bs-placement="bottom" data-bs-html="true"
|
|
title=" 6.3% of total for source files.">99 statements</a>
|
|
</li>
|
|
|
|
<li class="list-inline-item" id="source-file">
|
|
<i class="fa fa-code"></i>
|
|
<a href="../src/test_compact.f90"> Source File</a>
|
|
</li>
|
|
</ul>
|
|
</div>
|
|
<div class="col">
|
|
<nav aria-label="breadcrumb">
|
|
<ol class="breadcrumb justify-content-end mb-0">
|
|
<li class="breadcrumb-item active" aria-current="page">test_compact.f90</li>
|
|
</ol>
|
|
</nav>
|
|
</div>
|
|
</div>
|
|
</div>
|
|
<script>
|
|
// Enable Bootstrap tooltips
|
|
(function () {
|
|
const tooltipTriggerList = document.querySelectorAll('[data-bs-toggle="tooltip"]')
|
|
const tooltipList = [...tooltipTriggerList].map(tooltipTriggerEl => new bootstrap.Tooltip(tooltipTriggerEl))
|
|
})();
|
|
</script>
|
|
|
|
</div>
|
|
<div class="row">
|
|
<div class="col-md-3 hidden-xs hidden-sm visible-md visible-lg">
|
|
<div id="sidebar">
|
|
<h3>Contents</h3>
|
|
|
|
|
|
|
|
|
|
<div class="card mb-4">
|
|
<a data-bs-toggle="collapse" href="#progs-0"
|
|
aria-expanded="false" aria-controls="progs-0">
|
|
<h4 class="card-header bg-primary text-white">Programs</h4>
|
|
</a>
|
|
<div id="progs-0" class="collapse">
|
|
<div class="list-group list-group-flush">
|
|
<a class="list-group-item" href="../program/test_calculate.html">test_calculate</a>
|
|
</div>
|
|
</div>
|
|
</div>
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
<div class="card card-primary">
|
|
<div class="card-header text-left"><h3 class="card-title">Source Code</h3></div>
|
|
<div class="list-group">
|
|
<a class="list-group-item" href="../sourcefile/test_compact.f90.html#src">test_compact.f90</a>
|
|
</div>
|
|
</div>
|
|
|
|
|
|
</div>
|
|
|
|
</div>
|
|
<div class="col-md-9" id='text'>
|
|
|
|
<br>
|
|
|
|
<section>
|
|
<h2><span class="anchor" id="src"></span>Source Code</h2>
|
|
<div class="hl codehilite"><pre><span></span><a id="ln-1" name="ln-1" href="#ln-1"></a><span class="k">program </span><span class="n">test_calculate</span>
|
|
<a id="ln-2" name="ln-2" href="#ln-2"></a>
|
|
<a id="ln-3" name="ln-3" href="#ln-3"></a><span class="k">use </span><span class="n">Compact</span>
|
|
<a id="ln-4" name="ln-4" href="#ln-4"></a>
|
|
<a id="ln-5" name="ln-5" href="#ln-5"></a><span class="k">implicit none</span>
|
|
<a id="ln-6" name="ln-6" href="#ln-6"></a>
|
|
<a id="ln-7" name="ln-7" href="#ln-7"></a><span class="kt">real</span><span class="o">*</span><span class="mi">8</span><span class="p">,</span><span class="w"> </span><span class="k">parameter</span><span class="w"> </span><span class="kd">::</span><span class="w"> </span><span class="n">pi</span><span class="o">=</span><span class="mf">3.14159265358979323846d0</span>
|
|
<a id="ln-8" name="ln-8" href="#ln-8"></a><span class="kt">real</span><span class="o">*</span><span class="mi">8</span><span class="p">,</span><span class="w"> </span><span class="k">parameter</span><span class="w"> </span><span class="kd">::</span><span class="w"> </span><span class="n">me</span><span class="o">=</span><span class="mf">1.00d-20</span>
|
|
<a id="ln-9" name="ln-9" href="#ln-9"></a>
|
|
<a id="ln-10" name="ln-10" href="#ln-10"></a><span class="kt">integer</span><span class="w"> </span><span class="kd">::</span><span class="w"> </span><span class="n">nxp</span><span class="p">,</span><span class="n">nyp</span><span class="p">,</span><span class="n">nzp</span>
|
|
<a id="ln-11" name="ln-11" href="#ln-11"></a>
|
|
<a id="ln-12" name="ln-12" href="#ln-12"></a><span class="kt">real</span><span class="o">*</span><span class="mi">8</span><span class="w"> </span><span class="kd">::</span><span class="w"> </span><span class="n">hxp</span><span class="p">,</span><span class="n">hyp</span><span class="p">,</span><span class="n">hzp</span>
|
|
<a id="ln-13" name="ln-13" href="#ln-13"></a>
|
|
<a id="ln-14" name="ln-14" href="#ln-14"></a><span class="kt">real</span><span class="o">*</span><span class="mi">8</span><span class="w"> </span><span class="kd">::</span><span class="w"> </span><span class="n">l_0</span>
|
|
<a id="ln-15" name="ln-15" href="#ln-15"></a>
|
|
<a id="ln-16" name="ln-16" href="#ln-16"></a><span class="kt">real</span><span class="o">*</span><span class="mi">8</span><span class="p">,</span><span class="w"> </span><span class="k">allocatable</span><span class="p">,</span><span class="w"> </span><span class="k">dimension</span><span class="p">(:,:,:)</span><span class="w"> </span><span class="kd">::</span><span class="w"> </span><span class="n">aaa</span><span class="p">,</span><span class="n">bbb</span><span class="p">,</span><span class="n">ccc</span><span class="p">,</span><span class="n">ddd</span>
|
|
<a id="ln-17" name="ln-17" href="#ln-17"></a>
|
|
<a id="ln-18" name="ln-18" href="#ln-18"></a><span class="kt">real</span><span class="o">*</span><span class="mi">8</span><span class="p">,</span><span class="w"> </span><span class="k">allocatable</span><span class="p">,</span><span class="w"> </span><span class="k">dimension</span><span class="p">(:,:)</span><span class="w"> </span><span class="kd">::</span><span class="w"> </span><span class="n">aa</span><span class="p">,</span><span class="n">bb</span><span class="p">,</span><span class="n">cc</span><span class="p">,</span><span class="n">dd</span><span class="p">,</span><span class="n">ee</span>
|
|
<a id="ln-19" name="ln-19" href="#ln-19"></a>
|
|
<a id="ln-20" name="ln-20" href="#ln-20"></a><span class="kt">integer</span><span class="w"> </span><span class="kd">::</span><span class="w"> </span><span class="n">ierr</span>
|
|
<a id="ln-21" name="ln-21" href="#ln-21"></a>
|
|
<a id="ln-22" name="ln-22" href="#ln-22"></a><span class="kt">real</span><span class="o">*</span><span class="mi">8</span><span class="w"> </span><span class="kd">::</span><span class="w"> </span><span class="n">xx</span><span class="p">,</span><span class="n">yy</span><span class="p">,</span><span class="n">zz</span><span class="p">,</span><span class="n">fxyz</span>
|
|
<a id="ln-23" name="ln-23" href="#ln-23"></a>
|
|
<a id="ln-24" name="ln-24" href="#ln-24"></a><span class="kt">integer</span><span class="w"> </span><span class="kd">::</span><span class="w"> </span><span class="n">i</span><span class="p">,</span><span class="n">j</span><span class="p">,</span><span class="n">k</span>
|
|
<a id="ln-25" name="ln-25" href="#ln-25"></a>
|
|
<a id="ln-26" name="ln-26" href="#ln-26"></a><span class="kt">real</span><span class="o">*</span><span class="mi">4</span><span class="p">,</span><span class="w"> </span><span class="k">dimension</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span><span class="w"> </span><span class="kd">::</span><span class="w"> </span><span class="n">startt</span>
|
|
<a id="ln-27" name="ln-27" href="#ln-27"></a><span class="kt">real</span><span class="o">*</span><span class="mi">4</span><span class="p">,</span><span class="w"> </span><span class="k">dimension</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span><span class="w"> </span><span class="kd">::</span><span class="w"> </span><span class="n">endt</span>
|
|
<a id="ln-28" name="ln-28" href="#ln-28"></a><span class="kt">real</span><span class="o">*</span><span class="mi">4</span><span class="w"> </span><span class="kd">::</span><span class="w"> </span><span class="k">result</span>
|
|
<a id="ln-29" name="ln-29" href="#ln-29"></a>
|
|
<a id="ln-30" name="ln-30" href="#ln-30"></a><span class="n">nxp</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="mi">512</span>
|
|
<a id="ln-31" name="ln-31" href="#ln-31"></a><span class="n">nyp</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="mi">256</span>
|
|
<a id="ln-32" name="ln-32" href="#ln-32"></a><span class="n">nzp</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="mi">256</span>
|
|
<a id="ln-33" name="ln-33" href="#ln-33"></a>
|
|
<a id="ln-34" name="ln-34" href="#ln-34"></a><span class="n">l_0</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="mf">2.0d0</span>
|
|
<a id="ln-35" name="ln-35" href="#ln-35"></a><span class="n">hyp</span><span class="o">=</span><span class="n">l_0</span><span class="o">*</span><span class="n">pi</span><span class="o">/</span><span class="kt">REAL</span><span class="p">(</span><span class="n">nyp</span><span class="p">)</span>
|
|
<a id="ln-36" name="ln-36" href="#ln-36"></a><span class="n">hxp</span><span class="o">=</span><span class="n">hyp</span>
|
|
<a id="ln-37" name="ln-37" href="#ln-37"></a><span class="n">hzp</span><span class="o">=</span><span class="n">hyp</span>
|
|
<a id="ln-38" name="ln-38" href="#ln-38"></a>
|
|
<a id="ln-39" name="ln-39" href="#ln-39"></a><span class="k">allocate</span><span class="p">(</span><span class="n">aa</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span><span class="w"> </span><span class="n">nxp</span><span class="p">),</span><span class="w"> </span><span class="nb">stat</span><span class="o">=</span><span class="n">ierr</span><span class="p">);</span><span class="w"> </span><span class="n">aa</span><span class="o">=</span><span class="mf">0.0</span>
|
|
<a id="ln-40" name="ln-40" href="#ln-40"></a><span class="k">allocate</span><span class="p">(</span><span class="n">bb</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span><span class="w"> </span><span class="n">nxp</span><span class="p">),</span><span class="w"> </span><span class="nb">stat</span><span class="o">=</span><span class="n">ierr</span><span class="p">);</span><span class="w"> </span><span class="n">bb</span><span class="o">=</span><span class="mf">0.0</span>
|
|
<a id="ln-41" name="ln-41" href="#ln-41"></a><span class="k">allocate</span><span class="p">(</span><span class="n">cc</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span><span class="w"> </span><span class="n">nxp</span><span class="p">),</span><span class="w"> </span><span class="nb">stat</span><span class="o">=</span><span class="n">ierr</span><span class="p">);</span><span class="w"> </span><span class="n">cc</span><span class="o">=</span><span class="mf">0.0</span>
|
|
<a id="ln-42" name="ln-42" href="#ln-42"></a><span class="k">allocate</span><span class="p">(</span><span class="n">dd</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span><span class="w"> </span><span class="n">nxp</span><span class="p">),</span><span class="w"> </span><span class="nb">stat</span><span class="o">=</span><span class="n">ierr</span><span class="p">);</span><span class="w"> </span><span class="n">dd</span><span class="o">=</span><span class="mf">0.0</span>
|
|
<a id="ln-43" name="ln-43" href="#ln-43"></a><span class="k">allocate</span><span class="p">(</span><span class="n">ee</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span><span class="w"> </span><span class="n">nxp</span><span class="p">),</span><span class="w"> </span><span class="nb">stat</span><span class="o">=</span><span class="n">ierr</span><span class="p">);</span><span class="w"> </span><span class="n">ee</span><span class="o">=</span><span class="mf">0.0</span>
|
|
<a id="ln-44" name="ln-44" href="#ln-44"></a>
|
|
<a id="ln-45" name="ln-45" href="#ln-45"></a>
|
|
<a id="ln-46" name="ln-46" href="#ln-46"></a><span class="k">allocate</span><span class="p">(</span><span class="n">aaa</span><span class="p">(</span><span class="n">nxp</span><span class="p">,</span><span class="w"> </span><span class="n">nyp</span><span class="p">,</span><span class="w"> </span><span class="n">nzp</span><span class="p">),</span><span class="w"> </span><span class="nb">stat</span><span class="o">=</span><span class="n">ierr</span><span class="p">);</span><span class="w"> </span><span class="n">aaa</span><span class="o">=</span><span class="mf">0.0</span>
|
|
<a id="ln-47" name="ln-47" href="#ln-47"></a><span class="k">allocate</span><span class="p">(</span><span class="n">bbb</span><span class="p">(</span><span class="n">nxp</span><span class="p">,</span><span class="w"> </span><span class="n">nyp</span><span class="p">,</span><span class="w"> </span><span class="n">nzp</span><span class="p">),</span><span class="w"> </span><span class="nb">stat</span><span class="o">=</span><span class="n">ierr</span><span class="p">);</span><span class="w"> </span><span class="n">bbb</span><span class="o">=</span><span class="mf">0.0</span>
|
|
<a id="ln-48" name="ln-48" href="#ln-48"></a><span class="k">allocate</span><span class="p">(</span><span class="n">ccc</span><span class="p">(</span><span class="n">nxp</span><span class="p">,</span><span class="w"> </span><span class="n">nyp</span><span class="p">,</span><span class="w"> </span><span class="n">nzp</span><span class="p">),</span><span class="w"> </span><span class="nb">stat</span><span class="o">=</span><span class="n">ierr</span><span class="p">);</span><span class="w"> </span><span class="n">ccc</span><span class="o">=</span><span class="mf">0.0</span>
|
|
<a id="ln-49" name="ln-49" href="#ln-49"></a><span class="k">allocate</span><span class="p">(</span><span class="n">ddd</span><span class="p">(</span><span class="n">nxp</span><span class="p">,</span><span class="w"> </span><span class="n">nyp</span><span class="p">,</span><span class="w"> </span><span class="n">nzp</span><span class="p">),</span><span class="w"> </span><span class="nb">stat</span><span class="o">=</span><span class="n">ierr</span><span class="p">);</span><span class="w"> </span><span class="n">ddd</span><span class="o">=</span><span class="mf">0.0</span>
|
|
<a id="ln-50" name="ln-50" href="#ln-50"></a>
|
|
<a id="ln-51" name="ln-51" href="#ln-51"></a>
|
|
<a id="ln-52" name="ln-52" href="#ln-52"></a><span class="k">call </span><span class="n">ludcmp</span><span class="p">(</span><span class="n">nxp</span><span class="p">,</span><span class="w"> </span><span class="n">nyp</span><span class="p">,</span><span class="w"> </span><span class="n">nzp</span><span class="p">,</span><span class="w"> </span><span class="mi">1</span><span class="p">,</span><span class="w"> </span><span class="mi">0</span><span class="p">,</span><span class="w"> </span><span class="mi">0</span><span class="p">)</span>
|
|
<a id="ln-53" name="ln-53" href="#ln-53"></a>
|
|
<a id="ln-54" name="ln-54" href="#ln-54"></a><span class="k">do </span><span class="n">i</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="mi">1</span><span class="p">,</span><span class="n">nxp</span>
|
|
<a id="ln-55" name="ln-55" href="#ln-55"></a><span class="n">xx</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="n">dfloat</span><span class="p">(</span><span class="n">i</span><span class="o">-</span><span class="mi">1</span><span class="p">)</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="n">hxp</span>
|
|
<a id="ln-56" name="ln-56" href="#ln-56"></a><span class="n">aa</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span><span class="n">i</span><span class="p">)</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="nb">sin</span><span class="p">(</span><span class="mf">1.1d0</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="n">xx</span><span class="p">)</span>
|
|
<a id="ln-57" name="ln-57" href="#ln-57"></a><span class="n">dd</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span><span class="n">i</span><span class="p">)</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="mf">1.1d0</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="nb">cos</span><span class="p">(</span><span class="mf">1.1d0</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="n">xx</span><span class="p">)</span>
|
|
<a id="ln-58" name="ln-58" href="#ln-58"></a><span class="k">end do</span>
|
|
<a id="ln-59" name="ln-59" href="#ln-59"></a>
|
|
<a id="ln-60" name="ln-60" href="#ln-60"></a><span class="k">call </span><span class="n">rhs1np</span><span class="p">(</span><span class="n">nxp</span><span class="p">,</span><span class="w"> </span><span class="n">hxp</span><span class="p">,</span><span class="w"> </span><span class="n">aa</span><span class="p">,</span><span class="w"> </span><span class="n">bb</span><span class="p">,</span><span class="w"> </span><span class="mi">1</span><span class="p">)</span>
|
|
<a id="ln-61" name="ln-61" href="#ln-61"></a>
|
|
<a id="ln-62" name="ln-62" href="#ln-62"></a><span class="k">call </span><span class="n">dfnonp</span><span class="p">(</span><span class="n">nxp</span><span class="p">,</span><span class="w"> </span><span class="n">hxp</span><span class="p">,</span><span class="w"> </span><span class="n">aa</span><span class="p">,</span><span class="w"> </span><span class="n">cc</span><span class="p">,</span><span class="w"> </span><span class="mi">1</span><span class="p">,</span><span class="w"> </span><span class="mi">1</span><span class="p">)</span>
|
|
<a id="ln-63" name="ln-63" href="#ln-63"></a>
|
|
<a id="ln-64" name="ln-64" href="#ln-64"></a>
|
|
<a id="ln-65" name="ln-65" href="#ln-65"></a><span class="k">print</span><span class="w"> </span><span class="o">*</span><span class="p">,</span><span class="w"> </span><span class="n">lxf</span>
|
|
<a id="ln-66" name="ln-66" href="#ln-66"></a>
|
|
<a id="ln-67" name="ln-67" href="#ln-67"></a><span class="k">print</span><span class="w"> </span><span class="o">*</span><span class="p">,</span><span class="w"> </span><span class="n">lxs</span>
|
|
<a id="ln-68" name="ln-68" href="#ln-68"></a>
|
|
<a id="ln-69" name="ln-69" href="#ln-69"></a><span class="k">print</span><span class="w"> </span><span class="o">*</span><span class="p">,</span><span class="w"> </span><span class="n">aa</span>
|
|
<a id="ln-70" name="ln-70" href="#ln-70"></a>
|
|
<a id="ln-71" name="ln-71" href="#ln-71"></a><span class="k">print</span><span class="w"> </span><span class="o">*</span><span class="p">,</span><span class="w"> </span><span class="n">bb</span>
|
|
<a id="ln-72" name="ln-72" href="#ln-72"></a>
|
|
<a id="ln-73" name="ln-73" href="#ln-73"></a><span class="k">print</span><span class="w"> </span><span class="o">*</span><span class="p">,</span><span class="w"> </span><span class="n">cc</span>
|
|
<a id="ln-74" name="ln-74" href="#ln-74"></a>
|
|
<a id="ln-75" name="ln-75" href="#ln-75"></a><span class="k">print</span><span class="w"> </span><span class="o">*</span><span class="p">,</span><span class="w"> </span><span class="n">dd</span>
|
|
<a id="ln-76" name="ln-76" href="#ln-76"></a>
|
|
<a id="ln-77" name="ln-77" href="#ln-77"></a><span class="n">ee</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="p">((</span><span class="n">cc</span><span class="w"> </span><span class="o">-</span><span class="w"> </span><span class="n">dd</span><span class="p">)</span><span class="o">/</span><span class="n">dd</span><span class="p">)</span>
|
|
<a id="ln-78" name="ln-78" href="#ln-78"></a>
|
|
<a id="ln-79" name="ln-79" href="#ln-79"></a>
|
|
<a id="ln-80" name="ln-80" href="#ln-80"></a><span class="k">write</span><span class="p">(</span><span class="o">*</span><span class="p">,</span><span class="o">*</span><span class="p">)</span><span class="w"> </span><span class="s2">"relerr_min_max"</span><span class="p">,</span><span class="w"> </span><span class="nb">minval</span><span class="p">(</span><span class="n">ee</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span><span class="n">nxp</span><span class="o">/</span><span class="mi">4</span><span class="p">:</span><span class="mi">3</span><span class="o">*</span><span class="n">nxp</span><span class="o">/</span><span class="mi">4</span><span class="p">)),</span><span class="w"> </span><span class="nb">maxval</span><span class="p">(</span><span class="n">ee</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span><span class="n">nxp</span><span class="o">/</span><span class="mi">4</span><span class="p">:</span><span class="mi">3</span><span class="o">*</span><span class="n">nxp</span><span class="o">/</span><span class="mi">4</span><span class="p">))</span>
|
|
<a id="ln-81" name="ln-81" href="#ln-81"></a><span class="k">write</span><span class="p">(</span><span class="o">*</span><span class="p">,</span><span class="o">*</span><span class="p">)</span><span class="w"> </span><span class="s2">"relerr_min_max"</span><span class="p">,</span><span class="w"> </span><span class="nb">minval</span><span class="p">(</span><span class="n">ee</span><span class="p">),</span><span class="w"> </span><span class="nb">maxval</span><span class="p">(</span><span class="n">ee</span><span class="p">)</span>
|
|
<a id="ln-82" name="ln-82" href="#ln-82"></a><span class="k">write</span><span class="p">(</span><span class="o">*</span><span class="p">,</span><span class="o">*</span><span class="p">)</span><span class="w"> </span><span class="s2">" exact_min_max"</span><span class="p">,</span><span class="w"> </span><span class="nb">minval</span><span class="p">(</span><span class="n">dd</span><span class="p">),</span><span class="w"> </span><span class="nb">maxval</span><span class="p">(</span><span class="n">dd</span><span class="p">)</span>
|
|
<a id="ln-83" name="ln-83" href="#ln-83"></a><span class="k">write</span><span class="p">(</span><span class="o">*</span><span class="p">,</span><span class="o">*</span><span class="p">)</span><span class="w"> </span><span class="s2">"compact_min_max"</span><span class="p">,</span><span class="w"> </span><span class="nb">minval</span><span class="p">(</span><span class="n">cc</span><span class="p">),</span><span class="w"> </span><span class="nb">maxval</span><span class="p">(</span><span class="n">cc</span><span class="p">)</span>
|
|
<a id="ln-84" name="ln-84" href="#ln-84"></a>
|
|
<a id="ln-85" name="ln-85" href="#ln-85"></a>
|
|
<a id="ln-86" name="ln-86" href="#ln-86"></a><span class="k">write</span><span class="p">(</span><span class="o">*</span><span class="p">,</span><span class="o">*</span><span class="p">)</span><span class="w"> </span><span class="s2">"test ddx"</span>
|
|
<a id="ln-87" name="ln-87" href="#ln-87"></a>
|
|
<a id="ln-88" name="ln-88" href="#ln-88"></a><span class="k">do </span><span class="n">k</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="mi">1</span><span class="p">,</span><span class="n">nzp</span>
|
|
<a id="ln-89" name="ln-89" href="#ln-89"></a><span class="k">do </span><span class="n">j</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="mi">1</span><span class="p">,</span><span class="n">nyp</span>
|
|
<a id="ln-90" name="ln-90" href="#ln-90"></a><span class="k">do </span><span class="n">i</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="mi">1</span><span class="p">,</span><span class="n">nxp</span>
|
|
<a id="ln-91" name="ln-91" href="#ln-91"></a>
|
|
<a id="ln-92" name="ln-92" href="#ln-92"></a><span class="n">xx</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="n">i</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="n">hxp</span>
|
|
<a id="ln-93" name="ln-93" href="#ln-93"></a><span class="n">yy</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="n">j</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="n">hyp</span>
|
|
<a id="ln-94" name="ln-94" href="#ln-94"></a><span class="n">zz</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="n">k</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="n">hzp</span>
|
|
<a id="ln-95" name="ln-95" href="#ln-95"></a>
|
|
<a id="ln-96" name="ln-96" href="#ln-96"></a><span class="n">aaa</span><span class="p">(</span><span class="n">i</span><span class="p">,</span><span class="n">j</span><span class="p">,</span><span class="n">k</span><span class="p">)</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="nb">sin</span><span class="p">(</span><span class="mf">1.1</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="n">xx</span><span class="p">)</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="nb">sin</span><span class="p">(</span><span class="mf">3.0</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="n">yy</span><span class="p">)</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="nb">sin</span><span class="p">(</span><span class="mf">2.0</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="n">zz</span><span class="p">)</span>
|
|
<a id="ln-97" name="ln-97" href="#ln-97"></a>
|
|
<a id="ln-98" name="ln-98" href="#ln-98"></a><span class="k">end do</span>
|
|
<a id="ln-99" name="ln-99" href="#ln-99"></a><span class="k">end do</span>
|
|
<a id="ln-100" name="ln-100" href="#ln-100"></a><span class="k">end do</span>
|
|
<a id="ln-101" name="ln-101" href="#ln-101"></a>
|
|
<a id="ln-102" name="ln-102" href="#ln-102"></a>
|
|
<a id="ln-103" name="ln-103" href="#ln-103"></a>
|
|
<a id="ln-104" name="ln-104" href="#ln-104"></a>
|
|
<a id="ln-105" name="ln-105" href="#ln-105"></a>
|
|
<a id="ln-106" name="ln-106" href="#ln-106"></a><span class="k">do </span><span class="n">k</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="mi">1</span><span class="p">,</span><span class="n">nzp</span>
|
|
<a id="ln-107" name="ln-107" href="#ln-107"></a><span class="k">do </span><span class="n">j</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="mi">1</span><span class="p">,</span><span class="n">nyp</span>
|
|
<a id="ln-108" name="ln-108" href="#ln-108"></a>
|
|
<a id="ln-109" name="ln-109" href="#ln-109"></a><span class="k">call </span><span class="n">dfnonp</span><span class="p">(</span><span class="n">nxp</span><span class="p">,</span><span class="w"> </span><span class="n">hxp</span><span class="p">,</span><span class="w"> </span><span class="n">aaa</span><span class="p">(:,</span><span class="n">j</span><span class="p">,</span><span class="n">k</span><span class="p">),</span><span class="w"> </span><span class="n">bbb</span><span class="p">(:,</span><span class="n">j</span><span class="p">,</span><span class="n">k</span><span class="p">),</span><span class="w"> </span><span class="mi">1</span><span class="p">,</span><span class="w"> </span><span class="mi">1</span><span class="p">)</span>
|
|
<a id="ln-110" name="ln-110" href="#ln-110"></a>
|
|
<a id="ln-111" name="ln-111" href="#ln-111"></a><span class="k">end do</span>
|
|
<a id="ln-112" name="ln-112" href="#ln-112"></a><span class="k">end do</span>
|
|
<a id="ln-113" name="ln-113" href="#ln-113"></a>
|
|
<a id="ln-114" name="ln-114" href="#ln-114"></a><span class="k">do </span><span class="n">k</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="mi">1</span><span class="p">,</span><span class="n">nzp</span>
|
|
<a id="ln-115" name="ln-115" href="#ln-115"></a><span class="k">do </span><span class="n">j</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="mi">1</span><span class="p">,</span><span class="n">nyp</span>
|
|
<a id="ln-116" name="ln-116" href="#ln-116"></a><span class="k">do </span><span class="n">i</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="mi">1</span><span class="p">,</span><span class="n">nxp</span><span class="w"> </span><span class="c">! nxp/4,3*nxp/4</span>
|
|
<a id="ln-117" name="ln-117" href="#ln-117"></a>
|
|
<a id="ln-118" name="ln-118" href="#ln-118"></a><span class="n">xx</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="n">i</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="n">hxp</span>
|
|
<a id="ln-119" name="ln-119" href="#ln-119"></a><span class="n">yy</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="n">j</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="n">hyp</span>
|
|
<a id="ln-120" name="ln-120" href="#ln-120"></a><span class="n">zz</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="n">k</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="n">hzp</span>
|
|
<a id="ln-121" name="ln-121" href="#ln-121"></a><span class="n">fxyz</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="p">(</span><span class="mf">1.1</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="nb">cos</span><span class="p">(</span><span class="mf">1.1</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="n">xx</span><span class="p">)</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="nb">sin</span><span class="p">(</span><span class="mf">3.0</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="n">yy</span><span class="p">)</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="nb">sin</span><span class="p">(</span><span class="mf">2.0</span><span class="w"> </span><span class="o">*</span><span class="w"> </span><span class="n">zz</span><span class="p">))</span>
|
|
<a id="ln-122" name="ln-122" href="#ln-122"></a><span class="n">ccc</span><span class="p">(</span><span class="n">i</span><span class="p">,</span><span class="n">j</span><span class="p">,</span><span class="n">k</span><span class="p">)</span><span class="w"> </span><span class="o">=</span><span class="w"> </span><span class="p">(</span><span class="n">bbb</span><span class="p">(</span><span class="n">i</span><span class="p">,</span><span class="n">j</span><span class="p">,</span><span class="n">k</span><span class="p">)</span><span class="w"> </span><span class="o">-</span><span class="w"> </span><span class="n">fxyz</span><span class="p">)</span><span class="w"> </span><span class="o">/</span><span class="w"> </span><span class="p">(</span><span class="n">fxyz</span><span class="p">)</span>
|
|
<a id="ln-123" name="ln-123" href="#ln-123"></a>
|
|
<a id="ln-124" name="ln-124" href="#ln-124"></a><span class="k">if</span><span class="w"> </span><span class="p">((</span><span class="nb">abs</span><span class="p">(</span><span class="n">ccc</span><span class="p">(</span><span class="n">i</span><span class="p">,</span><span class="n">j</span><span class="p">,</span><span class="n">k</span><span class="p">))</span><span class="w"> </span><span class="o">></span><span class="w"> </span><span class="mf">0.01</span><span class="w"> </span><span class="p">)</span><span class="w"> </span><span class="p">.</span><span class="nb">and</span><span class="p">.</span><span class="w"> </span><span class="p">(</span><span class="nb">abs</span><span class="p">(</span><span class="n">bbb</span><span class="p">(</span><span class="n">i</span><span class="p">,</span><span class="n">j</span><span class="p">,</span><span class="n">k</span><span class="p">))</span><span class="w"> </span><span class="o">></span><span class="w"> </span><span class="mf">1.0e-14</span><span class="w"> </span><span class="p">))</span><span class="w"> </span><span class="k">write</span><span class="p">(</span><span class="o">*</span><span class="p">,</span><span class="o">*</span><span class="p">)</span><span class="w"> </span><span class="n">bbb</span><span class="p">(</span><span class="n">i</span><span class="p">,</span><span class="n">j</span><span class="p">,</span><span class="n">k</span><span class="p">),</span><span class="w"> </span><span class="n">fxyz</span>
|
|
<a id="ln-125" name="ln-125" href="#ln-125"></a>
|
|
<a id="ln-126" name="ln-126" href="#ln-126"></a><span class="k">end do</span>
|
|
<a id="ln-127" name="ln-127" href="#ln-127"></a><span class="k">end do</span>
|
|
<a id="ln-128" name="ln-128" href="#ln-128"></a><span class="k">end do</span>
|
|
<a id="ln-129" name="ln-129" href="#ln-129"></a>
|
|
<a id="ln-130" name="ln-130" href="#ln-130"></a><span class="k">write</span><span class="p">(</span><span class="o">*</span><span class="p">,</span><span class="o">*</span><span class="p">)</span><span class="w"> </span><span class="s2">"relerr_min_max"</span><span class="p">,</span><span class="w"> </span><span class="nb">minval</span><span class="p">(</span><span class="n">ccc</span><span class="p">),</span><span class="w"> </span><span class="nb">maxval</span><span class="p">(</span><span class="n">ccc</span><span class="p">)</span>
|
|
<a id="ln-131" name="ln-131" href="#ln-131"></a>
|
|
<a id="ln-132" name="ln-132" href="#ln-132"></a>
|
|
<a id="ln-133" name="ln-133" href="#ln-133"></a>
|
|
<a id="ln-134" name="ln-134" href="#ln-134"></a>
|
|
<a id="ln-135" name="ln-135" href="#ln-135"></a>
|
|
<a id="ln-136" name="ln-136" href="#ln-136"></a>
|
|
<a id="ln-137" name="ln-137" href="#ln-137"></a>
|
|
<a id="ln-138" name="ln-138" href="#ln-138"></a><span class="k">deallocate</span><span class="p">(</span><span class="n">aa</span><span class="p">)</span>
|
|
<a id="ln-139" name="ln-139" href="#ln-139"></a><span class="k">deallocate</span><span class="p">(</span><span class="n">bb</span><span class="p">)</span>
|
|
<a id="ln-140" name="ln-140" href="#ln-140"></a><span class="k">deallocate</span><span class="p">(</span><span class="n">cc</span><span class="p">)</span>
|
|
<a id="ln-141" name="ln-141" href="#ln-141"></a><span class="k">deallocate</span><span class="p">(</span><span class="n">dd</span><span class="p">)</span>
|
|
<a id="ln-142" name="ln-142" href="#ln-142"></a><span class="k">deallocate</span><span class="p">(</span><span class="n">ee</span><span class="p">)</span>
|
|
<a id="ln-143" name="ln-143" href="#ln-143"></a>
|
|
<a id="ln-144" name="ln-144" href="#ln-144"></a><span class="k">deallocate</span><span class="p">(</span><span class="n">aaa</span><span class="p">)</span>
|
|
<a id="ln-145" name="ln-145" href="#ln-145"></a><span class="k">deallocate</span><span class="p">(</span><span class="n">bbb</span><span class="p">)</span>
|
|
<a id="ln-146" name="ln-146" href="#ln-146"></a><span class="k">deallocate</span><span class="p">(</span><span class="n">ccc</span><span class="p">)</span>
|
|
<a id="ln-147" name="ln-147" href="#ln-147"></a><span class="k">deallocate</span><span class="p">(</span><span class="n">ddd</span><span class="p">)</span>
|
|
<a id="ln-148" name="ln-148" href="#ln-148"></a>
|
|
<a id="ln-149" name="ln-149" href="#ln-149"></a><span class="k">end program</span>
|
|
</pre></div>
|
|
|
|
</section>
|
|
</div>
|
|
</div>
|
|
|
|
<hr>
|
|
</div> <!-- /container -->
|
|
<footer>
|
|
<div class="container">
|
|
<div class="row justify-content-between">
|
|
<div class="col">
|
|
<p>
|
|
HPC DNS Post-Processing Fortran Core
|
|
was developed by Google DeepMind Team & Ignis<br> © 2026
|
|
</p>
|
|
</div>
|
|
<div class="col">
|
|
<p class="text-end">
|
|
Documentation generated by
|
|
<a href="https://github.com/Fortran-FOSS-Programmers/ford">FORD</a>
|
|
</p>
|
|
</div>
|
|
</div>
|
|
<br>
|
|
</div> <!-- /container -->
|
|
</footer>
|
|
</body>
|
|
</html> |