def get_model_syntest1():
model_ref = SeismicModel.ak135()
types = [ParameterType.VSH]
radii = 6371. - np.array([493.33, 410.])
model_params = ModelParameters(types, radii, mesh_type='boxcar')
model, mesh = model_ref.boxcar_mesh(model_params)
values = np.array(
[.2 * (-1)**i for i in range(model_params._n_grd_params)])
values_dict = {param_type: values for param_type in types}
values_mat = model_params.get_values_matrix(values_dict)
mesh_mul = mesh.multiply(model_params.get_nodes(), values_mat)
model_mul = model + mesh_mul
return model_mul
[0. if i % 2 == 1 else 0. for i in range(model_params._n_grd_params)])
values = np.array(
[0. if i % 2 == 1 else 0. for i in range(model_params._n_grd_params)])
values[2] = -0.1
values[3] = 0
values_r = np.zeros(model_params._n_grd_params, dtype=np.float64)
values_r[2] = -20
values_r[3] = -0
values_dict = {ParameterType.VSH: values, ParameterType.RADIUS: values_r}
values_dict_m = {
ParameterType.VSH: values_m,
ParameterType.RADIUS: values_r
}
values_mat = model_params.get_values_matrix(values_dict)
values_mat_m = model_params.get_values_matrix(values_dict_m)
model_ = model_.multiply(model_params.get_nodes(), values_mat,
values_mat_m)
# # mesh
# model, mesh = ak135.boxcar_mesh(model_params)
# # multiply mesh with values
# values = np.array(
# [0.1 * (-1)**i for i in range(model_params._n_grd_params)])
# values_dict = {
# ParameterType.VSH: values}
# values_mat = model_params.get_values_matrix(values_dict)
# mesh_ = mesh.multiply(model_params.get_nodes(), values_mat)
# model_ = model + mesh_
# # retrieve model perturbations