From d79db77e067cd369ed49fc6bed41270c89bac869 Mon Sep 17 00:00:00 2001
From: Andrew E Slaughter <andrew.slaughter@inl.gov>
Date: Mon, 15 Jan 2018 16:46:25 -0700
Subject: [PATCH] Create transfer for controlling layered material

(refs #107)
---
 include/controls/LayeredMaterialControl.h    |  31 +++++++
 src/base/MastodonApp.C                       |   6 ++
 src/controls/LayeredMaterialControl.C        |  27 ++++++
 src/materials/LinearSoilMaterial.C           |   2 +
 test/tests/controls/gold/layer_control_out.e | Bin 0 -> 27888 bytes
 test/tests/controls/layer_control.i          |  84 +++++++++++++++++++
 test/tests/controls/tests                    |  13 +++
 7 files changed, 163 insertions(+)
 create mode 100644 include/controls/LayeredMaterialControl.h
 create mode 100644 src/controls/LayeredMaterialControl.C
 create mode 100644 test/tests/controls/gold/layer_control_out.e
 create mode 100644 test/tests/controls/layer_control.i
 create mode 100644 test/tests/controls/tests

diff --git a/include/controls/LayeredMaterialControl.h b/include/controls/LayeredMaterialControl.h
new file mode 100644
index 0000000000..96f2bf234d
--- /dev/null
+++ b/include/controls/LayeredMaterialControl.h
@@ -0,0 +1,31 @@
+#ifndef LAYEREDMATERIALCONTROL_H
+#define LAYEREDMATERIALCONTROL_H
+
+// MOOSE includes
+#include "Control.h"
+
+// Forward Declarations
+class LayeredMaterialControl;
+
+template <>
+InputParameters validParams<LayeredMaterialControl>();
+
+/**
+ * A control for changing parameters associated with layer ids.
+ */
+class LayeredMaterialControl : public Control
+{
+public:
+  LayeredMaterialControl(const InputParameters & parameters);
+
+protected:
+  virtual void execute() override;
+
+  /// The values of the parameter to control
+  const std::vector<Real> & _current_values;
+
+  /// The values to use from this control
+  const std::vector<Real> & _control_values;
+};
+
+#endif
diff --git a/src/base/MastodonApp.C b/src/base/MastodonApp.C
index a14acce1fb..d95e57efed 100644
--- a/src/base/MastodonApp.C
+++ b/src/base/MastodonApp.C
@@ -45,6 +45,9 @@
 #include "ResponseSpectraCalculator.h"
 #include "HousnerSpectrumIntensity.h"
 
+// Controls
+#include "LayeredMaterialControl.h"
+
 // Testing
 #include "TestLayeredMaterialInterface.h"
 
@@ -123,6 +126,9 @@ MastodonApp::registerObjects(Factory & factory)
   registerVectorPostprocessor(ResponseSpectraCalculator);
   registerVectorPostprocessor(HousnerSpectrumIntensity);
 
+  // Controls
+  registerControl(LayeredMaterialControl);
+
   // Testing
   registerMaterial(TestLayeredMaterialInterfaceDocString);
   registerKernel(TestLayeredMaterialInterfaceTypeError);
diff --git a/src/controls/LayeredMaterialControl.C b/src/controls/LayeredMaterialControl.C
new file mode 100644
index 0000000000..8e631893bc
--- /dev/null
+++ b/src/controls/LayeredMaterialControl.C
@@ -0,0 +1,27 @@
+// Mastodon includes
+#include "LayeredMaterialControl.h"
+
+template <>
+InputParameters
+validParams<LayeredMaterialControl>()
+{
+  InputParameters params = validParams<Control>();
+  params.addRequiredParam<std::string>("parameter", "The name of the parameter to control.");
+  params.addRequiredParam<std::vector<Real>>("values", "The values to replace on the specified material object.");
+  return params;
+}
+
+LayeredMaterialControl::LayeredMaterialControl(const InputParameters & parameters)
+  : Control(parameters),
+  _current_values(getControllableValue<std::vector<Real>>("parameter")),
+  _control_values(getParam<std::vector<Real>>("values"))
+{
+  if (_current_values.size() != _control_values.size())
+    mooseError("The layer values specified in the ", name(), " Control object and the values being controlled must be the same size.");
+}
+
+void
+LayeredMaterialControl::execute()
+{
+  setControllableValue<std::vector<Real>>("parameter", _control_values);
+}
diff --git a/src/materials/LinearSoilMaterial.C b/src/materials/LinearSoilMaterial.C
index ca7ed2b5bf..90dec543da 100644
--- a/src/materials/LinearSoilMaterial.C
+++ b/src/materials/LinearSoilMaterial.C
@@ -21,6 +21,8 @@ validParams<LinearSoilMaterial>()
   params.addClassDescription("Material for computing the shear wave speed "
                              "and minimum element size as a function "
                              "of shear modulus and density.");
+
+  params.declareControllable("shear_modulus density");
   return params;
 }
 
diff --git a/test/tests/controls/gold/layer_control_out.e b/test/tests/controls/gold/layer_control_out.e
new file mode 100644
index 0000000000000000000000000000000000000000..79b8d6af2e16673c91deb52385f2300f09216342
GIT binary patch
literal 27888
zcmeHQPjlS172h=fB*$r*riUDQKqr|=I;rH{wIxfQZevB3UA3~5#)>bV4$CF65^FBG
zLvn4&mtOk;dhF!XnO=J6p@$xN=%I%m`w{XD+HcVQejvF^?vh+$#WdE$m{}17fgd06
z9v&Wn2Wa=s-EXy8t(Wk-g5Nv5y3#YkC~&+J+_hRiz;A_VuH#9>e}UgC1AiRwir;GF
z3?&N45pI&8@_Xix`C0yTZ#*<?XUNAd<G0Ft&$ngBlsIjLY0{O<M4GO|_eS4+%q)t=
zH`8)PD5In-^Q!dFu|Y>J{o4tK;}NE#aG3wC1jiW30Pj=}n0}S%F;2)two5jxo20#w
zmiJUPl`oUVY}2x=lkvo4d+{Bro{l&0k!1__SMhrtzc<X0V?2>T==dIPTd(5x((ON9
z>EZf%N&343*Eb_@${6{M7a33ezy|%#*zV*F>0iNbb>O&YoA^QInW`TOZr%U`&1lRf
zP0JzLdNpz)mtY8_{d=x?jux|gFA99u@W;`*diu`LJadNRDR^!`SWco-mPxS$;&0&h
zef+K{d3k#3C%~kEGLeSY5FyP7?mzKoRraUE;XbVsGcfQS{2{Ec+zRrK0)Le9k7;aU
zmU|_UN!I1%{X6)ZFVB3dUO4x@Rh2XT!*?b0t$c@Z<j=m3c^ADeE&nd^u7>mPV_!@|
z_I(10E6AUHAM@6I--pQin&m@E&Jc^+gzcRM#z_;xN*jOTyyq%`?N^rfr<1(T7UumG
za*#)41?5fk30xB%A{M8;uVSqgQpz^Y=4Ttv=D(BV|2^}gjk5XKMrHZG#Iq}GlSzsA
z{p^3UWo)D&tiF#wSE;`YDuSU&)i>KdEK9y1to{~%R?uraX+@Mv)|vf*&r{y5Dm@XM
zk23x3wY@v{_LMe>X`#F!to{yvuB!a`_fwCgf$vY`Z1oR}tNh_auf3UO#FcsZLs<P2
z{;bUCz_Pr5X5Iuh)s5dKN)ocTe?eM`yD0BpnKy3V3B7Nnptxpn|CY&{;wJ5+N@Vl?
z9Vu^8Fw{VkC}pIPA8+HAv;ksq+8q^Zt^Nbox1X49I4ReWeKHb3=@MMY&or=~Bs94H
zLOfr#Tf6L+C>!$`<A)KMecbbX<SXe-1L;n9eTCPn9AAt;TIlI%ANwxSDWhosPM__J
zaD}v{&v>1UXYW~{_3@YL()##6@%5|t>&ih?i5oa4r`0)HO<+{g(A-0|aVkq`eD<E{
z*=LNutj6RjI>J*4>3M2CfeA2@(k{-cu-_{P#o247z9)7AX-3i({c~~u=;*<oc&}#B
z;%o&3UKmj9`op2=+2Q~uuGk4q#zW~vp}@O}iKt8!YmUg-N$ZYSTMJM9r)!Df6arX0
zj7>eUkFHw%rn><&7EHUxDsM1^#pipsDru;G81o?^;Fg$%6DJZD+5(r&&2GQHxzXO}
z3*Qsnc4up?-C65w6e1ru{rfUJP3ts;-0Jt+o3?G<+}gfwZFjHB?M`2Anr7GRY_zY-
z?xx+huWwj$sGsf~Kd|QEgt)$bbA7XrhU&!KaS(u~;us$*Q}R&r`>5T|Oiy$+MCazM
zZueHVqiS0X*J|yZNox#^>`H|R^LaQ@z%zjDTiwlCfQw%E^*^W7?t7zgB<|vKqFQFZ
z9%>)`vf=j%;drJMIZx!YljLGOq4!21Xf8n#4GhgP$YH6Qp4j>P_^3Xp5G6n<v_7CR
z{>oAH0cq5x<P=<c;@-~Z4<6L6bK_jD{pPwEAMK3K9-4te^$<2KlXZ!x%qPU7^=oRB
zYKT1YweR}Y;{pb&q@LJ6e7JLve05h6ryjD<z#Ka6c?n_<_Qc`Q{r#PT#bC@qn9AO)
z-<JvKM2_ilO2kGCgGnYVQ|0o02jHIAJvw}FymNS5zsi`TfyxduWM!huP?%PE5coqw
zO)POnhZBH`!h@aTy$3bvOPLb4;S>;GJ1LR+6h2yCzt)gYvz0G_w)tfD`BX9!A&s>S
zzUP|a3#Fs)k3HKAiq%6;Ed2DFP{#-X(`n7fJelt<`f`Y&-c$A=vFAfiBohcZI?8W|
z;0#BuW5GHv3QGZ2SvO1R!|<s8Ys_-(z-Y=VLy)`-g_X8LPuy_=$2xUo^vP~{(@mbR
z*~X#7Qbag%O-tGa6xV5{F?*st$y_JGqEJ?9pwGZeqdhUMgEPs)lD1=>7$HB`o!oGu
zP?wRzGOXnRF^ixI7yVhT4c&36iKs6&4G#em?lcuq`cO?@HKY*9hXwku#%x)pdIWr-
z57`RS@IoJE1f)31#(JR-ry|4HSD+7d2dRej;UzGMcjG1M=hrS)LqMuwZRv$hbbg`1
zfTGs?=k5R5oI9X@fwb>q#Q;W?n{csM*VHoFQ#mMu$fvUnZ8?&j9U61iD<3XneKr0S
zAl1~fET?e)xb~z-i?wR}vtnd&$pSni`fR|wO~T%5bSW}}Yl-;hc`><Td91yT<@`rk
z4BZnSbZ)lWn8I(jKbRCNkFe`CIz>?aswJP~l*7=s$L{!Inv4zgLOI0ZlK#bHS~5t~
z;}XlgP!4J9#aI1LGa)pUL%e#Obb>ifk&H#k;l_G*Lw?lW5}ozyjBHjZiN=+~B`}ua
zwM@*{$lyXINXnfBx`xAOnzh~)fK`UbGBZn+^GXCf4onzjSgET;q7Y`^53s5S-&x<q
z)Xl(j+-zUb{|dlDJSOfr2^3F{2NDZ&W2`mw<-iXloB$=fU5&mK*80R)@IAi_Q|#i2
z$h6w}%E6A{NSYR(zXDKxbiMw^5A=03ESN)^AooO9u61tI;xIEUn9^NF4q<B=TtAu6
zl#l3F*)_~QJxb-si2R%gD&(G6)0%e{!i9t5nc<W#A$8ga8g1wn5;mzftGw2pa|tPZ
z0VsqS3D0jcu(ETx89-YO%&}XqCzH#p@Yw)v+vAZ9TXlkzM0(;zcl*Zn=FRSQEiwT(
z7=X>C%iulnbW~cuf(Js-Mc)jf@knuc&N9~5HeR*dD5h-aa5(`il0i1(Qh?;kVC+HT
z%?AK#${{GMlI)()s4tiVrnvCAbW8i*QdK5kWgE7)3IwD$2^9LqK!Qk)Df3NOYaGc!
z1CSdgplgaaf)}95+58LRunw_-!!rC4myZfS^*_LKfEY~iVVbS^MFExZc|kzOGO=%i
zYP*Qj`eHMR0n3DY38ut=WwNL!)RIFWnP<w65r&c#`7;v*`V=ZatmmbZ&oiQ<Z*a!O
zRjGxeuV4X4jU9VT^vy6Qh09o9p>+yRvaIOk;=ReV&YB#%GFa62goq8-(osskr(luQ
ziMTx}N7&?M3Bb@mrLD(>{UjMvwB!)`M5-)AeJ;m5E;Z=%0SU!zYL0IkoR^%mUL4VH
z!z`epHV=6rz&zANC`|2FiLWp#FXZ}0d&EyU1H4wkGt+_;oBDPQdr|B|FoZ$rn4^&~
zuOsvUvA?hH&#Htep}K(I$T<gEHrN;km&@U_<)$&T;n~fBl-=%Dr*osT&}tGQ18B@E
zY>JH*;3Z(n{^f8aQ=3p7tBJWP%EzUj7{7~8BHGLQb$p;&ta9U?x(l37^j&|`rO`5H
zyj|{Vs~cG5#+TvvAzI+%#&>Z{aRwir{kXlp@uSZ6RvUkAZ2z!=T$XW8JkINleIs6G
zke)SlZqywD2$&kfjL!vlNU|SQ(3F4lO`I72xePq%UTo<t---?2sf9w@84Q}xg`|yQ
za_RNq9W=-|Jx5_-i?hm&>%#@=D-o=Pw8a@nVdGAr(mBDQ9*>oT^98xEjJ9c9!8Y!e
z@cjG1KJU)S4LeR;0+TrPRe5d)0W}1~jx@~-0#pQ-ehA=BZZmk`J8sIrDqYJeH$DV>
z1tzgt5GJwnL4V3rEeK-<h<j1z_Zh%~Jv=(xt7^6UGm^Am)Js`W&vYkX;u*&1e4yVK
zT)wbxyRP5YSMU|k_#hP@9pKE(^r`x@%J`fR{G6!rVo6!ULNQ>;w}Adv05%+}9fj%;
zA{j(Z-Y-vo{xAhBnmh#gJt5EF0as6EP;XW_Y<ho6&IyTc$ObYc@hlRluTd3qShnGT
z{u?+hW8+pc<6YcW7L0k6!?aV?yV7KjcbUb9iDV^4#r1oN8lnuRZ0pssfU`@Y6H1o=
zO!j1!5P1Q>Vj}ecF+|~T!?6{`=lr=&h$T7hWLC#b@fMt}zaAjO$*GUU3foAA`x$iz
zFx|PYc|fdM4_vic3fAZYHtstY;$SJbEtPFm*J(fYCFtwN&d7jYD2^M%H0Iobp*~>g
zB4#YxIB`P4*^%cG@!Xn)9UA)DbKrD0gzLnBX9Ud#eF{z2^`8O`4itT??#-@|>YH<Y
z#|({f=eZ>TrxTF#Rn~RwP!v9!#hKjeQ)t0EGw>aIDmplv(spp}nVOf??<rH|Fwq9F
za!^O2B+8+FwHFw~%Z(|8i(6?U>bIdfy|28x=(l0zG!Uy?zx`t6aAJ<JbnE(kGm3oU
zB$yfh^t+x?u*r$39;U%g;klbJK9>Ykr?8rlb*j!;O&P+gejCzr$Cc(_4l0I48`dA8
zl!Eo(DH1xFN@nRJg<5#6*zYSHsMH5W9$TNwUb*TbtyuJ5sHL+RKw-J+iUD=y&?{G6
z%K>sDgzM73uh20Rha2qO-TC~W?pZ2oXGJ`MmTkP!^@AN1$}?-L0{mdL<c1WogT!97
z@Eha^e3Ic)Pk0ip5Lrf6#)UXTS4qh!xa!N{S*64H5<TU-X$sGqrtn$Q6h3Q~!pksA
z-0;`vC8+XIe6Cbu-{9CmGvdx$qsZf)U0Z*U$^aU)-p#zvv;MfO{KwQ$>jk=I844}m
zo%MaqBI~O!TCg}-JqwbJF9pc6-t)uz6jGllR`_IEd_JZRX!`@4z&($Yr04@y?yn}7
zkN5B6h_k&<jlIL<gy^bK3X{%i0<(fbT~kDzx>n$lHmS1qoWrnTrkr@=k?ZICTc!c@
znR1H%?vMePjq0DXRJp`OE?4n%{u$>FvEd!t0cJ@{?)AG~<3@WjaqkS2JEQtH4yH5R
zWudPRN0_*X)ISEn0n_RJ#Hp&#uW}5?L%HCqYkDV`SL0FUX_M-UO$-SC=>fhV=l#(2
z0eOydN>0I5zYU?T$DFP(PBHiB&mCvC1cm3GrM47K8IdJaU$7%-HQSe;M5U*Ml2ZLD
zpV?-%inH-$0B3fj<RNtdFM&Dt@VP#&7W!8CIgA@KH{*rE#uxpSS5;nErgyqx!;%Y*
zz&+Tx_u0<jy_)Ov1gsqQ^pyhv<NfJe6HxQ_d`lYi_2Gg8t$X6`{%7^imu-=d#*e@a
zEfXF$^TUa%wf?9#ZNupixZI}0O)!A$i5F+;&Kq30l=T0Q-2UR<|7E!SMHpXyb6e%R
J{gXKR{{Z7R5CQ-I

literal 0
HcmV?d00001

diff --git a/test/tests/controls/layer_control.i b/test/tests/controls/layer_control.i
new file mode 100644
index 0000000000..9f617c26fe
--- /dev/null
+++ b/test/tests/controls/layer_control.i
@@ -0,0 +1,84 @@
+[Mesh]
+  type = GeneratedMesh
+  dim = 1
+  nx = 2
+[]
+
+[Variables]
+  [./u]
+  [../]
+[]
+
+[Kernels]
+  [./diff]
+    type = Diffusion
+    variable = u
+  [../]
+[]
+
+[BCs]
+  [./left]
+    type = DirichletBC
+    boundary = left
+    variable = u
+    value = 1
+  [../]
+  [./right]
+    type = DirichletBC
+    boundary = right
+    variable = u
+    value = 0
+  [../]
+[]
+
+[AuxVariables]
+  [./layer]
+    order = CONSTANT
+    family = MONOMIAL
+  [../]
+[]
+
+[ICs]
+  [./layer_ic]
+    type = FunctionIC
+    function = 'if(x<0.5,1980,1949)'
+    variable = layer
+  [../]
+[]
+
+[Materials]
+  [./linear]
+    type = LinearSoilMaterial
+    layer_variable = layer
+    layer_ids     = '1980   1949'
+    shear_modulus = '1.35e6 4.25e8'
+    density       = '1500   1700'
+    outputs = exodus
+  [../]
+[]
+
+[Controls]
+  [./control_density]
+    type = LayeredMaterialControl
+    parameter = Material/linear/density
+    values = '1700 1500'
+    execute_on = 'initial'
+  [../]
+  [./control_modulus]
+    type = LayeredMaterialControl
+    parameter = Material/linear/shear_modulus
+    values = '4.25e8 1.35e6'
+    execute_on = 'initial'
+  [../]
+[]
+
+[Executioner]
+  type = Transient
+  num_steps = 1
+[]
+
+[Outputs]
+  execute_on = 'TIMESTEP_END'
+  hide = 'u'
+  exodus = true
+[]
diff --git a/test/tests/controls/tests b/test/tests/controls/tests
new file mode 100644
index 0000000000..12b7ae5756
--- /dev/null
+++ b/test/tests/controls/tests
@@ -0,0 +1,13 @@
+[Tests]
+  [./control_layer_values]
+    type = Exodiff
+    input = layer_control.i
+    exodiff = layer_control_out.e
+  [../]
+  [./invalid_value_size]
+    type = RunException
+    input = layer_control.i
+    cli_args = Controls/control_density/values='1500'
+    expect_err = "The layer values specified in the control_density Control object and the values being controlled must be the same size."
+  [../]
+[]