<?xml version="1.0" encoding="ISO-8859-1"?><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">
<front>
<journal-meta>
<journal-id>0122-5383</journal-id>
<journal-title><![CDATA[CT&F - Ciencia, Tecnología y Futuro]]></journal-title>
<abbrev-journal-title><![CDATA[C.T.F Cienc. Tecnol. Futuro]]></abbrev-journal-title>
<issn>0122-5383</issn>
<publisher>
<publisher-name><![CDATA[Instituto Colombiano del Petróleo (ICP) - ECOPETROL S.A.]]></publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id>S0122-53832019000100119</article-id>
<article-id pub-id-type="doi">10.29047/01225383.159</article-id>
<title-group>
<article-title xml:lang="en"><![CDATA[SOLUTION OF A P AND S WAVE PROPAGATION MODEL USING HIGH PERFORMANCE COMPUTATION]]></article-title>
<article-title xml:lang="es"><![CDATA[SOLUCIÓN DE UN MODELO DE PROPAGACIÓN DE ONDAS P Y S, USANDO COMPUTACIÓN DE ALTO RENDIMIENTO]]></article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname><![CDATA[Amado]]></surname>
<given-names><![CDATA[Jonathan]]></given-names>
</name>
<xref ref-type="aff" rid="Aff"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname><![CDATA[Piedrahita]]></surname>
<given-names><![CDATA[Carlos]]></given-names>
</name>
<xref ref-type="aff" rid="Aff"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname><![CDATA[Vera-Ciro]]></surname>
<given-names><![CDATA[Carlos]]></given-names>
</name>
<xref ref-type="aff" rid="Aff"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname><![CDATA[Parra]]></surname>
<given-names><![CDATA[Carlos]]></given-names>
</name>
<xref ref-type="aff" rid="Aff"/>
</contrib>
</contrib-group>
<aff id="Af1">
<institution><![CDATA[,Universidad de Medellín  ]]></institution>
<addr-line><![CDATA[Medellín ]]></addr-line>
<country>Colombia</country>
</aff>
<aff id="Af2">
<institution><![CDATA[,Universidad de Pamplona  ]]></institution>
<addr-line><![CDATA[Pamplona ]]></addr-line>
<country>Colombia</country>
</aff>
<pub-date pub-type="pub">
<day>00</day>
<month>06</month>
<year>2019</year>
</pub-date>
<pub-date pub-type="epub">
<day>00</day>
<month>06</month>
<year>2019</year>
</pub-date>
<volume>9</volume>
<numero>1</numero>
<fpage>119</fpage>
<lpage>130</lpage>
<copyright-statement/>
<copyright-year/>
<self-uri xlink:href="http://www.scielo.org.co/scielo.php?script=sci_arttext&amp;pid=S0122-53832019000100119&amp;lng=en&amp;nrm=iso"></self-uri><self-uri xlink:href="http://www.scielo.org.co/scielo.php?script=sci_abstract&amp;pid=S0122-53832019000100119&amp;lng=en&amp;nrm=iso"></self-uri><self-uri xlink:href="http://www.scielo.org.co/scielo.php?script=sci_pdf&amp;pid=S0122-53832019000100119&amp;lng=en&amp;nrm=iso"></self-uri><abstract abstract-type="short" xml:lang="en"><p><![CDATA[ABSTRACT The propagation of seismic waves is affected by the type of transmission media. Therefore, it is necessary to solve a differential equation system in partial derivatives allowing for identification of waves propagating into an elastic media. This paper summarizes a research using a partial differential equation system representing the wave equation using the finite differences method to obtain the elastic media response, using an staggered grid. To prevent reflections in the computational regions, absorbent boundaries were used with the PML method. The implementation of the numerical scheme was made on two computational architectures (CPU and GPU) that share the same type of memory distribution. Finally, different code versions were created to take advantage of the architecture in the GPU memory, performing a detailed analysis of variables such as usage of bandwidth of the GPU internal memory, added to a version that is not limited by the internal memory in the graphic processing unit, but rather by the memory of the whole computational system.]]></p></abstract>
<abstract abstract-type="short" xml:lang="es"><p><![CDATA[RESUMEN La propagación de las ondas sísmicas se ven afectadas por el tipo de medio donde se transmiten. Por lo tanto se requiere solucionar un sistema de ecuaciones diferenciales en derivada parcial que permitan identificar las ondas que se propagan en el medio elástico, Este documento resume una investigación donde se usó un sistema de ecuaciones diferenciales parciales que representan la ecuación de onda utilizando el método de diferencias finitas para obtener la respuesta del medio elástico, usando una malla intercalada. Para evitar las reflexiones en las fronteras computaciones se utilizaron fronteras absorbentes con el método PML. La implementación del esquema numérico se realizó en dos arquitecturas de cómputo (CPU y GPU) que comparten el mismo tipo de distribución de memoria. Finalmente, se crearon diferentes versiones del código para aprovechar la arquitectura en la memoria de la GPU, se realiza una análisis detallado de variables como utilización y ancho de banda de la memoria interna de la GPU, además de una versión la cual no se ve limitada por la cantidad de memoria interna en la unidad de procesamiento grafico sino por la cantidad de memoria que tiene todo el sistema de cómputo.]]></p></abstract>
<kwd-group>
<kwd lng="en"><![CDATA[Modelling]]></kwd>
<kwd lng="en"><![CDATA[Elastic media]]></kwd>
<kwd lng="en"><![CDATA[PML]]></kwd>
<kwd lng="en"><![CDATA[GPU constant memory]]></kwd>
<kwd lng="en"><![CDATA[GPU shared memory]]></kwd>
<kwd lng="en"><![CDATA[Asynchronous copies and executions]]></kwd>
<kwd lng="es"><![CDATA[Modelado]]></kwd>
<kwd lng="es"><![CDATA[Medio elástico]]></kwd>
<kwd lng="es"><![CDATA[PML]]></kwd>
<kwd lng="es"><![CDATA[Memoria constante GPU]]></kwd>
<kwd lng="es"><![CDATA[Memoria compartida GPU]]></kwd>
<kwd lng="es"><![CDATA[Copias y ejecuciones no síncronas]]></kwd>
</kwd-group>
</article-meta>
</front><back>
<ref-list>
<ref id="B1">
<label>[1]</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Komatitsch]]></surname>
<given-names><![CDATA[D.]]></given-names>
</name>
<name>
<surname><![CDATA[Erlebacher]]></surname>
<given-names><![CDATA[G.]]></given-names>
</name>
<name>
<surname><![CDATA[Goddeke]]></surname>
<given-names><![CDATA[D.]]></given-names>
</name>
<name>
<surname><![CDATA[Michea]]></surname>
<given-names><![CDATA[D]]></given-names>
</name>
</person-group>
<article-title xml:lang=""><![CDATA[High-order finite-element seismic wave propagation modeling with MPI on a large CPL cluster]]></article-title>
<source><![CDATA[Journal of Computational Physics]]></source>
<year>2010</year>
<numero>229</numero>
<issue>229</issue>
<page-range>7692-714</page-range></nlm-citation>
</ref>
<ref id="B2">
<label>[2]</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Weiss]]></surname>
<given-names><![CDATA[R.]]></given-names>
</name>
<name>
<surname><![CDATA[Shragge]]></surname>
<given-names><![CDATA[J]]></given-names>
</name>
</person-group>
<article-title xml:lang=""><![CDATA[Solving 3C anisotropic elastic wave equations on paralell GP0 devices]]></article-title>
<source><![CDATA[Geophysics]]></source>
<year>2013</year>
<volume>78</volume>
<page-range>1-9</page-range></nlm-citation>
</ref>
<ref id="B3">
<label>[3]</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Das]]></surname>
<given-names><![CDATA[S]]></given-names>
</name>
<name>
<surname><![CDATA[Chen]]></surname>
<given-names><![CDATA[X]]></given-names>
</name>
<name>
<surname><![CDATA[Hobson]]></surname>
<given-names><![CDATA[M]]></given-names>
</name>
</person-group>
<article-title xml:lang=""><![CDATA[Fast GP0-cased seismogram simulation from microseismic events n marine environments using heterogeneous velocity models.]]></article-title>
<source><![CDATA[IEEE transactions on computational imaging]]></source>
<year>2017</year>
<volume>3</volume>
<numero>1 2</numero>
<issue>1 2</issue>
<page-range>316-29</page-range></nlm-citation>
</ref>
<ref id="B4">
<label>[4]</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Virieux]]></surname>
<given-names><![CDATA[J]]></given-names>
</name>
</person-group>
<article-title xml:lang=""><![CDATA[A p-sv wave propagation in heterogeneous media: Velocity-stress finite difference method]]></article-title>
<source><![CDATA[Geophysics]]></source>
<year>1986</year>
<volume>51</volume>
<page-range>899-904</page-range></nlm-citation>
</ref>
<ref id="B5">
<label>[5]</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Fornberg]]></surname>
<given-names><![CDATA[B]]></given-names>
</name>
</person-group>
<article-title xml:lang=""><![CDATA[Generation of finite difference formulas on arbitrarily spaced grids]]></article-title>
<source><![CDATA[Mathematics of Computation]]></source>
<year>1988</year>
<volume>51</volume>
<page-range>699-706</page-range></nlm-citation>
</ref>
<ref id="B6">
<label>[6]</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Bamberger]]></surname>
<given-names><![CDATA[A]]></given-names>
</name>
<name>
<surname><![CDATA[Chavent]]></surname>
<given-names><![CDATA[G]]></given-names>
</name>
<name>
<surname><![CDATA[Lailly]]></surname>
<given-names><![CDATA[P]]></given-names>
</name>
</person-group>
<article-title xml:lang=""><![CDATA[Étude de schémas numérigues de l'élastodynamigue linéaire.]]></article-title>
<source><![CDATA[Technical report, INRIA]]></source>
<year>2006</year>
</nlm-citation>
</ref>
<ref id="B7">
<label>[7]</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Berenger]]></surname>
</name>
</person-group>
<article-title xml:lang=""><![CDATA[A perfectly matched layer for the absorption of electromagnetic waves]]></article-title>
<source><![CDATA[Journal of Computational Geophysics]]></source>
<year>1994</year>
<volume>114</volume>
<page-range>185-200</page-range></nlm-citation>
</ref>
<ref id="B8">
<label>[8]</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Zeng]]></surname>
</name>
<name>
<surname><![CDATA[Liu]]></surname>
</name>
</person-group>
<article-title xml:lang=""><![CDATA[The application of the perfectly matched layer in numerical modeling of wave propagation in poroelastic media]]></article-title>
<source><![CDATA[Geophysics]]></source>
<year>2001</year>
<volume>66</volume>
<page-range>1258-66</page-range></nlm-citation>
</ref>
<ref id="B9">
<label>[9]</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Mahrer]]></surname>
</name>
</person-group>
<article-title xml:lang=""><![CDATA[An empirical study of instability and improvement of absorbing boundary conditions for elastic wave equation]]></article-title>
<source><![CDATA[Geophysics]]></source>
<year>1986</year>
<volume>51</volume>
<page-range>1499-507</page-range></nlm-citation>
</ref>
<ref id="B10">
<label>[10]</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Stacey]]></surname>
</name>
</person-group>
<article-title xml:lang=""><![CDATA[Improved transparent boundary formulations for the elastic wave equation]]></article-title>
<source><![CDATA[Bulletin of Seismological Society of America]]></source>
<year>1988</year>
<volume>78</volume>
<page-range>2080-97</page-range></nlm-citation>
</ref>
<ref id="B11">
<label>[11]</label><nlm-citation citation-type="">
<collab>NVIDIA Corporation</collab>
<source><![CDATA[CUDA NVIDIA visual profiler]]></source>
<year>2017</year>
</nlm-citation>
</ref>
</ref-list>
</back>
</article>
