Point source merge

.vscode/settings.json

@@ -10,7 +10,9 @@
     },
     "files.associations": {
         "*.icc": "cpp",
-        "iosfwd": "cpp"
+        "iosfwd": "cpp",
+        "*.tcc": "cpp",
+        "fstream": "cpp"
     },
     "python.formatting.autopep8Args": [
         "-max-line-length",

configure.ac

@@ -199,6 +199,7 @@ AC_CONFIG_FILES([Makefile
 		libsrc/pylith/feassemble/Makefile
 		libsrc/pylith/fekernels/Makefile
 		libsrc/pylith/faults/Makefile
+		libsrc/pylith/sources/Makefile
 		libsrc/pylith/friction/Makefile
 		libsrc/pylith/materials/Makefile
 		libsrc/pylith/meshio/Makefile
@@ -211,6 +212,7 @@ AC_CONFIG_FILES([Makefile
 		modulesrc/include/Makefile
 		modulesrc/bc/Makefile
 		modulesrc/faults/Makefile
+		modulesrc/sources/Makefile
 		modulesrc/feassemble/Makefile
 		modulesrc/friction/Makefile
 		modulesrc/materials/Makefile
@@ -236,6 +238,7 @@ AC_CONFIG_FILES([Makefile
 		tests/libtests/problems/data/Makefile
 		tests/libtests/meshio/data/Makefile
 		tests/libtests/scales/Makefile
+		tests/libtests/sources/Makefile
 		tests/libtests/topology/Makefile
 		tests/libtests/topology/data/Makefile
 		tests/libtests/testing/Makefile
@@ -250,6 +253,8 @@ AC_CONFIG_FILES([Makefile
 		tests/mmstests/incompressibleelasticity/nofaults-2d/Makefile
 		tests/mmstests/poroelasticity/Makefile
 		tests/mmstests/poroelasticity/nofaults-2d/Makefile
+		tests/mmstests/sources/Makefile
+		tests/mmstests/sources/momenttensor-2d/Makefile
 		tests/fullscale/Makefile
 		tests/fullscale/cornercases/Makefile
 		tests/fullscale/cornercases/nofaults-2d/Makefile

docs/user/components/index.md

@@ -14,6 +14,7 @@ faults/index.md
 materials/index.md
 meshio/index.md
 problems/index.md
+sources/index.md
 testing/index.md
 topology/index.md
 utils/index.md

docs/user/components/sources/AuxSubfieldsMomentTensorForce.md

@@ -0,0 +1,29 @@
+# AuxSubfieldsMomentTensorForce
+
+% WARNING: Do not edit; this is a generated file!
+:Full name: `pylith.sources.AuxSubfieldsMomentTensorForce`
+:Journal name: `auxsubfieldsmomenttensorforce`
+
+Container for moment tensor force auxiliary subfields.
+
+This component defines the discretization for the moment tensor and time delay auxiliary subfields.
+
+## Pyre Facilities
+
+* `moment_tensor`: Moment tensor subfield.
+  - **current value**: 'subfield', from {default}
+  - **configurable as**: subfield, moment_tensor
+* `time_delay`: Time delay subfield.
+  - **current value**: 'subfield', from {default}
+  - **configurable as**: subfield, time_delay
+
+## Example
+
+Example of setting `AuxSubfieldsMomentTensorForce` Pyre properties and facilities in a parameter file.
+
+:::{code-block} cfg
+[pylithapp.problem.sources.source]
+auxiliary_subfields.moment_tensor.basis_order = 0
+auxiliary_subfields.time_delay.basis_order = 0
+:::
+

docs/user/components/sources/AuxSubfieldsSourceTime.md

@@ -0,0 +1,25 @@
+# AuxSubfieldsSourceTime
+
+% WARNING: Do not edit; this is a generated file!
+:Full name: `pylith.sources.AuxSubfieldsSourceTime`
+:Journal name: `auxsubfieldssourcetime`
+
+Container for source time function auxiliary subfields.
+
+This component defines the discretization for the center frequency auxiliary subfield used by wavelet source time functions.
+
+## Pyre Facilities
+
+* `center_frequency`: Center frequency subfield.
+  - **current value**: 'subfield', from {default}
+  - **configurable as**: subfield, center_frequency
+
+## Example
+
+Example of setting `AuxSubfieldsSourceTime` Pyre properties and facilities in a parameter file.
+
+:::{code-block} cfg
+[pylithapp.problem.sources.source]
+source_time_function.auxiliary_subfields.center_frequency.basis_order = 0
+:::
+

docs/user/components/sources/GaussianWavelet.md

@@ -0,0 +1,41 @@
+# GaussianWavelet
+
+% WARNING: Do not edit; this is a generated file!
+:Full name: `pylith.sources.GaussianWavelet`
+:Journal name: `gaussianwavelet`
+
+Gaussian wavelet source time function.
+
+The Gaussian wavelet provides a smooth pulse.
+
+The wavelet is defined as:
+
+$S(t) = \frac{1}{2\pi^2 f_0^2} \exp\left(-\pi^2 f_0^2 (t-t_d)^2\right)$
+
+where $f_0$ is the center frequency and $t_d$ is the time delay.
+
+Implements `SourceTimeFunctionMomentTensorForce`.
+
+## Pyre Facilities
+
+* `auxiliary_subfields`: Discretization information for source time function parameters.
+  - **current value**: 'auxiliary_subfields', from {file='...', line=..., function='__set__'}
+  - **configurable as**: auxiliary_subfields
+
+## Example
+
+Example of setting `GaussianWavelet` Pyre properties and facilities in a parameter file.
+
+:::{code-block} cfg
+[pylithapp.problem.sources.source]
+source_time_function = pylith.sources.GaussianWavelet
+
+# Center frequency specified in the auxiliary field database
+db_auxiliary_field = spatialdata.spatialdb.UniformDB
+db_auxiliary_field.description = Source properties
+db_auxiliary_field.values = [moment_tensor_xx, moment_tensor_yy, moment_tensor_xy, moment_tensor_zz, time_delay, center_frequency]
+db_auxiliary_field.data = [1.0e12*Pa*s, 1.0e12*Pa*s, 0.0*Pa*s, 1.0e12*Pa*s, 0.0*s, 5.0]
+
+source_time_function.auxiliary_subfields.center_frequency.basis_order = 0
+:::
+

docs/user/components/sources/MomentTensorForce.md

@@ -0,0 +1,84 @@
+# MomentTensorForce
+
+% WARNING: Do not edit; this is a generated file!
+:Full name: `pylith.sources.MomentTensorForce`
+:Journal name: `momenttensorforce`
+
+Moment tensor point source.
+
+This component implements a point source using a moment tensor representation.
+The moment tensor describes the equivalent body forces of a seismic source.
+The temporal evolution of the source is controlled by a source time function.
+
+Implements `Source`.
+
+## Pyre Facilities
+
+* `auxiliary_subfields`: Discretization information for auxiliary subfields.
+  - **current value**: 'auxiliary_subfields', from {file='...', line=..., function='__set__'}
+  - **configurable as**: auxiliary_subfields
+* `db_auxiliary_field`: Database for source parameters.
+  - **current value**: 'simpledb', from {default}
+  - **configurable as**: simpledb, db_auxiliary_field
+* `derived_subfields`: Discretization of derived subfields.
+  - **current value**: 'emptybin', from {default}
+  - **configurable as**: emptybin, derived_subfields
+* `observers`: Observers (e.g., output).
+  - **current value**: 'singlephysicsobserver', from {default}
+  - **configurable as**: singlephysicsobserver, observers
+* `reader`: Reader for source points list.
+  - **current value**: 'pointslist', from {default}
+  - **configurable as**: pointslist, reader
+* `source_time_function`: Source time function for moment tensor force.
+  - **current value**: 'timehistorywavelet', from {default}
+  - **configurable as**: timehistorywavelet, source_time_function
+
+## Pyre Properties
+
+* `description`=\<str\>: Descriptive label for source.
+  - **default value**: ''
+  - **current value**: '', from {default}
+  - **validator**: <function validateDescription at ...>
+* `field`=\<str\>: Solution subfield associated with source.
+  - **default value**: 'displacement'
+  - **current value**: 'displacement', from {default}
+* `label`=\<str\>: Name of label for source.
+  - **default value**: 'source-id'
+  - **current value**: 'source-id', from {default}
+* `label_value`=\<int\>: Value of label identifying source.
+  - **default value**: 1
+  - **current value**: 1, from {default}
+
+## Example
+
+Example of setting `MomentTensorForce` Pyre properties and facilities in a parameter file.
+
+:::{code-block} cfg
+[pylithapp.problem]
+sources = [source]
+sources.source = pylith.sources.MomentTensorForce
+
+[pylithapp.problem.sources.source]
+description = Earthquake source
+label_value = 2
+
+# Specify source locations
+reader.filename = source_sites.txt
+reader.coordsys = spatialdata.geocoords.CSCart
+reader.coordsys.space_dim = 2
+
+# Use Ricker wavelet source time function
+source_time_function = pylith.sources.RickerWavelet
+
+# Source parameters (moment tensor and time delay)
+db_auxiliary_field = spatialdata.spatialdb.UniformDB
+db_auxiliary_field.description = Source properties
+db_auxiliary_field.values = [moment_tensor_xx, moment_tensor_yy, moment_tensor_xy, moment_tensor_zz, time_delay, center_frequency]
+db_auxiliary_field.data = [1.0e12*Pa*s, 1.0e12*Pa*s, 0.0*Pa*s, 1.0e12*Pa*s, 0.0*s, 5.0]
+
+# Discretization of auxiliary subfields
+auxiliary_subfields.moment_tensor.basis_order = 0
+auxiliary_subfields.time_delay.basis_order = 0
+source_time_function.auxiliary_subfields.center_frequency.basis_order = 0
+:::
+

docs/user/components/sources/RickerWavelet.md

@@ -0,0 +1,42 @@
+# RickerWavelet
+
+% WARNING: Do not edit; this is a generated file!
+:Full name: `pylith.sources.RickerWavelet`
+:Journal name: `rickerwavelet`
+
+Ricker wavelet source time function.
+
+The Ricker wavelet (Mexican hat wavelet) is the negative normalized second derivative of a Gaussian.
+It is commonly used as a source time function in seismic wave propagation studies.
+
+The wavelet is defined as:
+
+$S(t) = \left(1 - 2\pi^2 f_0^2 (t-t_d)^2\right) \exp\left(-\pi^2 f_0^2 (t-t_d)^2\right)$
+
+where $f_0$ is the center frequency and $t_d$ is the time delay.
+
+Implements `SourceTimeFunctionMomentTensorForce`.
+
+## Pyre Facilities
+
+* `auxiliary_subfields`: Discretization information for source time function parameters.
+  - **current value**: 'auxiliary_subfields', from {file='...', line=..., function='__set__'}
+  - **configurable as**: auxiliary_subfields
+
+## Example
+
+Example of setting `RickerWavelet` Pyre properties and facilities in a parameter file.
+
+:::{code-block} cfg
+[pylithapp.problem.sources.source]
+source_time_function = pylith.sources.RickerWavelet
+
+# Center frequency specified in the auxiliary field database
+db_auxiliary_field = spatialdata.spatialdb.UniformDB
+db_auxiliary_field.description = Source properties
+db_auxiliary_field.values = [moment_tensor_xx, moment_tensor_yy, moment_tensor_xy, moment_tensor_zz, time_delay, center_frequency]
+db_auxiliary_field.data = [1.0e12*Pa*s, 1.0e12*Pa*s, 0.0*Pa*s, 1.0e12*Pa*s, 0.0*s, 5.0]
+
+source_time_function.auxiliary_subfields.center_frequency.basis_order = 0
+:::
+

docs/user/components/sources/Source.md

@@ -0,0 +1,45 @@
+# Source
+
+% WARNING: Do not edit; this is a generated file!
+:Full name: `pylith.sources.Source`
+:Journal name: `source`
+
+Abstract base class for point sources.
+
+This is the base class for all point source implementations in PyLith.
+Point sources are used to specify internal forcing terms at discrete points within the domain.
+
+## Pyre Facilities
+
+* `auxiliary_subfields`: Discretization information for auxiliary subfields.
+  - **current value**: 'emptybin', from {default}
+  - **configurable as**: emptybin, auxiliary_subfields
+* `db_auxiliary_field`: Database for source parameters.
+  - **current value**: 'simpledb', from {default}
+  - **configurable as**: simpledb, db_auxiliary_field
+* `derived_subfields`: Discretization of derived subfields.
+  - **current value**: 'emptybin', from {default}
+  - **configurable as**: emptybin, derived_subfields
+* `observers`: Observers (e.g., output).
+  - **current value**: 'singlephysicsobserver', from {default}
+  - **configurable as**: singlephysicsobserver, observers
+* `reader`: Reader for source points list.
+  - **current value**: 'pointslist', from {default}
+  - **configurable as**: pointslist, reader
+
+## Pyre Properties
+
+* `description`=\<str\>: Descriptive label for source.
+  - **default value**: ''
+  - **current value**: '', from {default}
+  - **validator**: <function validateDescription at ...>
+* `field`=\<str\>: Solution subfield associated with source.
+  - **default value**: 'displacement'
+  - **current value**: 'displacement', from {default}
+* `label`=\<str\>: Name of label for source.
+  - **default value**: 'source-id'
+  - **current value**: 'source-id', from {default}
+* `label_value`=\<int\>: Value of label identifying source.
+  - **default value**: 1
+  - **current value**: 1, from {default}
+

docs/user/components/sources/SourceTimeFunctionMomentTensorForce.md

@@ -0,0 +1,17 @@
+# SourceTimeFunctionMomentTensorForce
+
+% WARNING: Do not edit; this is a generated file!
+:Full name: `pylith.sources.SourceTimeFunctionMomentTensorForce`
+:Journal name: `sourcetimefunctionmomenttensorforce`
+
+Abstract base class for source time functions.
+
+This is the base class for source time functions used with moment tensor force sources.
+Derived classes implement specific wavelet types.
+
+## Pyre Facilities
+
+* `auxiliary_subfields`: Discretization information for physical properties and state variables.
+  - **current value**: 'emptybin', from {default}
+  - **configurable as**: emptybin, auxiliary_subfields
+

docs/user/components/sources/SquareWavelet.md

@@ -0,0 +1,41 @@
+# SquareWavelet
+
+% WARNING: Do not edit; this is a generated file!
+:Full name: `pylith.sources.SquareWavelet`
+:Journal name: `squarewavelet`
+
+Square wavelet (step function) source time function.
+
+The square wavelet is a Heaviside step function that activates the source at the time delay.
+
+The wavelet is defined as:
+
+$S(t) = \begin{cases} 0 & t < t_d \\ 1 & t \geq t_d \end{cases}$
+
+where $t_d$ is the time delay.
+
+Implements `SourceTimeFunctionMomentTensorForce`.
+
+## Pyre Facilities
+
+* `auxiliary_subfields`: Discretization information for source time function parameters.
+  - **current value**: 'auxiliary_subfields', from {file='...', line=..., function='__set__'}
+  - **configurable as**: auxiliary_subfields
+
+## Example
+
+Example of setting `SquareWavelet` Pyre properties and facilities in a parameter file.
+
+:::{code-block} cfg
+[pylithapp.problem.sources.source]
+source_time_function = pylith.sources.SquareWavelet
+
+# Source parameters specified in the auxiliary field database
+db_auxiliary_field = spatialdata.spatialdb.UniformDB
+db_auxiliary_field.description = Source properties
+db_auxiliary_field.values = [moment_tensor_xx, moment_tensor_yy, moment_tensor_xy, moment_tensor_zz, time_delay, center_frequency]
+db_auxiliary_field.data = [1.0e12*Pa*s, 1.0e12*Pa*s, 0.0*Pa*s, 1.0e12*Pa*s, 0.5*s, 1.0]
+
+source_time_function.auxiliary_subfields.center_frequency.basis_order = 0
+:::
+