Method MakeDefaultImageSystemFromMesh

src/gui/frame.cpp:2641–2696

// ---------------------------------------------------------------------

void MyFrame::MakeDefaultImageSystemFromMesh(vtkUnstructuredGrid* ug)
{
    size_t num_chemicals = 2;
    size_t target_chemical = 1;
    size_t largest_dimension = 64;
    float value_inside = 0.0f;
    float value_outside = 1.0f;

    unique_ptr<ImageRD> image_sys;
    if (this->is_opencl_available)
    {
        IntegerDialog nc_dlg(this, _("Number of chemicals in new volume image"), _("Number of chemicals:"),
            3, 1, 20, wxDefaultPosition, wxDefaultSize);
        if (nc_dlg.ShowModal() != wxID_OK) return;
        num_chemicals = nc_dlg.GetValue();

        wxArrayString choices;
        for (size_t i = 0; i<num_chemicals; i++)
            choices.Add(GetChemicalName(i));
        wxSingleChoiceDialog tc_dlg(this, _("Select the chemical to write the volume into:"), _("Select target chemical"),
            choices);
        tc_dlg.SetSelection(2);
        if (tc_dlg.ShowModal() != wxID_OK) return;
        target_chemical = tc_dlg.GetSelection();

        IntegerDialog ld_dlg(this, _("Largest dimension of the new volume image"), _("Largest dimension:"),
            32, 1, 1024, wxDefaultPosition, wxDefaultSize);
        if (ld_dlg.ShowModal() != wxID_OK) return;
        largest_dimension = ld_dlg.GetValue();

        FloatDialog inval_dlg(this, _("Value to set inside the mesh"), _("Value inside:"),
            0.1f, wxDefaultPosition, wxDefaultSize);
        if (inval_dlg.ShowModal() != wxID_OK) return;
        value_inside = inval_dlg.GetValue();

        FloatDialog outval_dlg(this, _("Value to set outside the mesh"), _("Value outside:"),
            0.5f, wxDefaultPosition, wxDefaultSize);
        if (outval_dlg.ShowModal() != wxID_OK) return;
        value_outside = outval_dlg.GetValue();

        // at some point we would want the user to decide what data type to use in the image
        const int data_type = VTK_FLOAT;

        image_sys = make_unique<FormulaOpenCLImageRD>(opencl_platform, opencl_device, data_type);
        image_sys->SetFormula("delta_a = D_a * laplacian_a - a*b*b + F*(1.0f-a);\ndelta_b = D_b * laplacian_b + a*b*b - (F+K+c*0.035f)*b;");
    }
    else {
        image_sys = make_unique<GrayScottImageRD>();
    }

    image_sys->CopyFromMesh(ug, num_chemicals, target_chemical, largest_dimension, value_inside, value_outside);
    image_sys->CreateDefaultInitialPatternGenerator(2);
    image_sys->GenerateInitialPattern();
    this->render_settings.GetProperty("timesteps_per_render").SetInt(16);
    this->SetCurrentRDSystem(move(image_sys));
}

// ---------------------------------------------------------------------