NetworkManager.C 223 KB
Newer Older
hrchilds's avatar
hrchilds committed
1 2
/*****************************************************************************
*
brugger's avatar
 
brugger committed
3
* Copyright (c) 2000 - 2011, Lawrence Livermore National Security, LLC
hrchilds's avatar
hrchilds committed
4
* Produced at the Lawrence Livermore National Laboratory
brugger's avatar
 
brugger committed
5
* LLNL-CODE-442911
hrchilds's avatar
hrchilds committed
6 7
* All rights reserved.
*
brugger's avatar
 
brugger committed
8
* This file is  part of VisIt. For  details, see https://visit.llnl.gov/.  The
hrchilds's avatar
hrchilds committed
9 10 11 12 13 14 15 16 17 18
* full copyright notice is contained in the file COPYRIGHT located at the root
* of the VisIt distribution or at http://www.llnl.gov/visit/copyright.html.
*
* Redistribution  and  use  in  source  and  binary  forms,  with  or  without
* modification, are permitted provided that the following conditions are met:
*
*  - Redistributions of  source code must  retain the above  copyright notice,
*    this list of conditions and the disclaimer below.
*  - Redistributions in binary form must reproduce the above copyright notice,
*    this  list of  conditions  and  the  disclaimer (as noted below)  in  the
brugger's avatar
 
brugger committed
19 20 21
*    documentation and/or other materials provided with the distribution.
*  - Neither the name of  the LLNS/LLNL nor the names of  its contributors may
*    be used to endorse or promote products derived from this software without
hrchilds's avatar
hrchilds committed
22 23 24 25 26
*    specific prior written permission.
*
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT  HOLDERS AND CONTRIBUTORS "AS IS"
* AND ANY EXPRESS OR  IMPLIED WARRANTIES, INCLUDING,  BUT NOT  LIMITED TO, THE
* IMPLIED WARRANTIES OF MERCHANTABILITY AND  FITNESS FOR A PARTICULAR  PURPOSE
brugger's avatar
 
brugger committed
27 28 29
* ARE  DISCLAIMED. IN  NO EVENT  SHALL LAWRENCE  LIVERMORE NATIONAL  SECURITY,
* LLC, THE  U.S.  DEPARTMENT OF  ENERGY  OR  CONTRIBUTORS BE  LIABLE  FOR  ANY
* DIRECT,  INDIRECT,   INCIDENTAL,   SPECIAL,   EXEMPLARY,  OR   CONSEQUENTIAL
hrchilds's avatar
hrchilds committed
30 31 32 33 34 35 36 37 38
* DAMAGES (INCLUDING, BUT NOT  LIMITED TO, PROCUREMENT OF  SUBSTITUTE GOODS OR
* SERVICES; LOSS OF  USE, DATA, OR PROFITS; OR  BUSINESS INTERRUPTION) HOWEVER
* CAUSED  AND  ON  ANY  THEORY  OF  LIABILITY,  WHETHER  IN  CONTRACT,  STRICT
* LIABILITY, OR TORT  (INCLUDING NEGLIGENCE OR OTHERWISE)  ARISING IN ANY  WAY
* OUT OF THE  USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH
* DAMAGE.
*
*****************************************************************************/

hrchilds's avatar
hrchilds committed
39 40 41 42
// ************************************************************************* //
//                              NetworkManager.C                             //
// ************************************************************************* //

hrchilds's avatar
hrchilds committed
43
#include <snprintf.h>
hrchilds's avatar
hrchilds committed
44
#include <AttributeSubject.h>
hrchilds's avatar
hrchilds committed
45 46
#include <NetworkManager.h>
#include <DataNetwork.h>
hrchilds's avatar
hrchilds committed
47
#include <ClonedDataNetwork.h>
48
#include <CumulativeQueryNamedSelectionExtension.h>
hrchilds's avatar
hrchilds committed
49 50 51
#include <DebugStream.h>
#include <avtDatabaseFactory.h>
#include <LoadBalancer.h>
52
#include <ConstructDataBinningAttributes.h>
hrchilds's avatar
hrchilds committed
53 54
#include <DBOptionsAttributes.h>
#include <ExportDBAttributes.h>
hrchilds's avatar
hrchilds committed
55
#include <MaterialAttributes.h>
hrchilds's avatar
hrchilds committed
56
#include <MeshManagementAttributes.h>
hrchilds's avatar
hrchilds committed
57
#include <avtExpressionEvaluatorFilter.h>
58 59
#include <avtNamedSelectionFilter.h>
#include <avtNamedSelectionManager.h>
hrchilds's avatar
hrchilds committed
60
#include <ImproperUseException.h>
hrchilds's avatar
hrchilds committed
61
#include <InvalidVariableException.h>
hrchilds's avatar
hrchilds committed
62 63 64 65 66
#include <DatabaseException.h>
#include <NoInputException.h>
#include <avtPluginFilter.h>
#include <OperatorPluginManager.h>
#include <OperatorPluginInfo.h>
67
#include <AnnotationObject.h>
hrchilds's avatar
hrchilds committed
68
#include <AnnotationObjectList.h>
fogal1's avatar
fogal1 committed
69
#include <AnnotationAttributes.h>
hrchilds's avatar
hrchilds committed
70
#include <PickAttributes.h>
71
#include <SelectionProperties.h>
hrchilds's avatar
hrchilds committed
72 73
#include <VisualCueInfo.h>
#include <VisualCueList.h>
74
#include <avtApplyDataBinningExpression.h>
hrchilds's avatar
hrchilds committed
75
#include <avtCallback.h>
hrchilds's avatar
hrchilds committed
76
#include <avtColorTables.h>
77
#include <avtDebugDumpOptions.h>
hrchilds's avatar
hrchilds committed
78
#include <avtExtents.h>
hrchilds's avatar
hrchilds committed
79
#include <avtNullData.h>
80
#include <avtDatabaseMetaData.h>
hrchilds's avatar
hrchilds committed
81
#include <avtDataObjectQuery.h>
hrchilds's avatar
hrchilds committed
82
#include <avtMultipleInputQuery.h>
hrchilds's avatar
hrchilds committed
83
#include <avtAreaBetweenCurvesQuery.h>
hrchilds's avatar
hrchilds committed
84
#include <avtFileWriter.h>
hrchilds's avatar
hrchilds committed
85
#include <avtDatabaseWriter.h>
hrchilds's avatar
hrchilds committed
86
#include <avtL2NormBetweenCurvesQuery.h>
hrchilds's avatar
hrchilds committed
87
#include <avtLocateQuery.h>
hrchilds's avatar
hrchilds committed
88
#include <avtLocateCellQuery.h>
hrchilds's avatar
hrchilds committed
89 90 91 92
#include <avtLocateNodeQuery.h>
#include <avtActualCoordsQuery.h>
#include <avtActualNodeCoordsQuery.h>
#include <avtActualZoneCoordsQuery.h>
hrchilds's avatar
hrchilds committed
93
#include <avtPickQuery.h>
94
#include <avtNamedSelection.h>
hrchilds's avatar
hrchilds committed
95
#include <avtNodePickQuery.h>
hrchilds's avatar
hrchilds committed
96
#include <avtParallel.h>
hrchilds's avatar
hrchilds committed
97 98
#include <avtPickByNodeQuery.h>
#include <avtPickByZoneQuery.h>
99
#include <avtTransparencyActor.h>
hrchilds's avatar
hrchilds committed
100
#include <avtZonePickQuery.h>
hrchilds's avatar
hrchilds committed
101
#include <avtCurvePickQuery.h>
hrchilds's avatar
hrchilds committed
102
#include <avtSoftwareShader.h>
103 104
#include <avtDataBinning.h>
#include <avtDataBinningConstructor.h>
hrchilds's avatar
hrchilds committed
105 106 107
#include <avtSourceFromAVTImage.h>
#include <avtSourceFromImage.h>
#include <avtSourceFromNullData.h>
hrchilds's avatar
hrchilds committed
108
#include <avtTiledImageCompositor.h>
hrchilds's avatar
hrchilds committed
109 110
#include <avtWholeImageCompositerWithZ.h>
#include <avtWholeImageCompositerNoZ.h>
hrchilds's avatar
hrchilds committed
111
#include <avtPlot.h>
hrchilds's avatar
hrchilds committed
112 113
#include <avtQueryOverTimeFilter.h>
#include <avtQueryFactory.h>
hrchilds's avatar
hrchilds committed
114
#include <CompactSILRestrictionAttributes.h>
hrchilds's avatar
hrchilds committed
115
#include <VisWindow.h>
hrchilds's avatar
hrchilds committed
116
#include <VisWinRendering.h> // for SetStereoEnabled
hrchilds's avatar
hrchilds committed
117
#include <ParsingExprList.h>
hrchilds's avatar
hrchilds committed
118
#include <avtExprNode.h>
hrchilds's avatar
hrchilds committed
119 120
#include <DatabasePluginManager.h>
#include <DatabasePluginInfo.h>
hrchilds's avatar
hrchilds committed
121 122
#include <PlotPluginManager.h>
#include <PlotPluginInfo.h>
fogal1's avatar
fogal1 committed
123
#include <StackTimer.h>
124 125 126

#include <vtkImageData.h>

hrchilds's avatar
hrchilds committed
127 128 129
#ifdef PARALLEL
#include <mpi.h>
#endif
hrchilds's avatar
hrchilds committed
130
#include <visit-config.h>
hrchilds's avatar
hrchilds committed
131
#include <TimingsManager.h>
132
#include <SaveWindowAttributes.h>
hrchilds's avatar
hrchilds committed
133

bonnell's avatar
bonnell committed
134
#include <algorithm>
135
#include <climits>
136 137 138
#include <map>
#include <string>
#include <vector>
hrchilds's avatar
hrchilds committed
139

hrchilds's avatar
hrchilds committed
140 141 142
//
// Static functions.
//
hrchilds's avatar
hrchilds committed
143 144
static ref_ptr<avtDatabase> GetDatabase(void *, const std::string &,
                                        int, const char *);
145
static avtDataBinning *GetDataBinningCallbackBridge(void *arg, const char *name);
bonnell's avatar
bonnell committed
146
static bool OnlyRootNodeHasData(avtImage_p &);
147
static void BroadcastImage(avtImage_p &, bool, int);
148
#ifdef PARALLEL
149
static bool IsBlankImage(avtImage_p img);
150 151
static std::vector<int> BuildBlankImageVector(avtImage_p img);
#endif
hrchilds's avatar
hrchilds committed
152

hrchilds's avatar
hrchilds committed
153 154 155 156 157 158 159 160 161 162
//
// Static data members of the NetworkManager class.
//
InitializeProgressCallback NetworkManager::initializeProgressCallback = NULL;
void                      *NetworkManager::initializeProgressCallbackArgs=NULL;
ProgressCallback           NetworkManager::progressCallback = NULL;
void                      *NetworkManager::progressCallbackArgs = NULL;



hrchilds's avatar
hrchilds committed
163 164 165 166 167 168 169 170 171
// ****************************************************************************
//  Method: NetworkManager default constructor
//
//  Programmer: Jeremy Meredith
//  Creation:   September 29, 2000
//
//  Modifications:
//
//    Kathleen Bonnell, Tue Nov 20 12:35:54 PST 2001
172
//    Initialize requireOriginalCells.
hrchilds's avatar
hrchilds committed
173
//
hrchilds's avatar
hrchilds committed
174 175 176
//    Hank Childs, Mon Jan  5 16:04:57 PST 2004
//    Initialize uniqueNetworkId.
//
hrchilds's avatar
hrchilds committed
177 178 179
//    Hank Childs, Thu Mar  4 09:11:46 PST 2004
//    Change instantiation of vis window, since it is now assumed to be nowin.
//
hrchilds's avatar
hrchilds committed
180 181 182
//    Mark C. Miller, Thu Apr 29 16:15:25 PDT 2004
//    Added missing initialization of dumpRenders
//
hrchilds's avatar
hrchilds committed
183 184 185
//    Mark C. Miller, Tue May 25 20:44:10 PDT 2004
//    Removed unnecessary scope. Set annotationAttributes data member
//
186 187
//    Kathleen Bonnell, Wed Jun  2 09:48:29 PDT 2004
//    Initialize requireOriginalNodes.
hrchilds's avatar
hrchilds committed
188
//
hrchilds's avatar
hrchilds committed
189 190 191 192
//    Mark C. Miller, Tue Jan  4 10:23:19 PST 2005
//    Made the viswin part of a map of EngineVisWinInfo objects
//    Removed some unnecessary initialization
//
hrchilds's avatar
hrchilds committed
193 194 195
//    Mark C. Miller, Wed Jan  5 10:14:21 PST 2005
//    Added call to NewVisWindow
//
hrchilds's avatar
hrchilds committed
196 197 198
//    Hank Childs, Mon Feb 28 16:56:20 PST 2005
//    Initialize inQueryMode.
//
hrchilds's avatar
hrchilds committed
199 200 201
//    Hank Childs, Fri Aug 26 15:44:48 PDT 2005
//    Register a callback to get databases.
//
hrchilds's avatar
hrchilds committed
202 203 204
//    Hank Childs, Sat Feb 18 11:31:23 PST 2006
//    Added a call back for the apply ddf expression.
//
205 206 207
//    Cyrus Harrison, Tue Feb 19 08:42:51 PST 2008
//    Removed dumpRenders (now controled by avtDebugDumpOptions)
//
208 209 210
//    Brad Whitlock, Tue Jun 24 15:40:29 PDT 2008
//    Added plugin managers.
//
fogal1's avatar
fogal1 committed
211 212 213
//    Tom Fogal, Wed Dec  9 14:10:01 MST 2009
//    Remove creation of window 0; we'll do it dynamically when needed.
//
214 215 216
//    Hank Childs, Sat Aug 21 14:35:47 PDT 2010
//    Rename DDF to DataBinning.
//
hrchilds's avatar
hrchilds committed
217 218 219 220 221 222
// ****************************************************************************
NetworkManager::NetworkManager(void) : virtualDatabases()
{
    workingNet = NULL;
    loadBalancer = NULL;
    requireOriginalCells = false;
hrchilds's avatar
hrchilds committed
223
    requireOriginalNodes = false;
hrchilds's avatar
hrchilds committed
224
    inQueryMode = false;
hrchilds's avatar
hrchilds committed
225
    uniqueNetworkId = 0;
hrchilds's avatar
hrchilds committed
226

hrchilds's avatar
hrchilds committed
227
    avtCallback::RegisterGetDatabaseCallback(GetDatabase, this);
228 229 230
    avtApplyDataBinningExpression::RegisterGetDataBinningCallback(GetDataBinningCallbackBridge, this);
    avtExpressionEvaluatorFilter::RegisterGetDataBinningCallback(
                                                  GetDataBinningCallbackBridge, this);
231 232 233 234

    databasePlugins = new DatabasePluginManager;
    operatorPlugins = new OperatorPluginManager;
    plotPlugins = new PlotPluginManager;
hrchilds's avatar
hrchilds committed
235 236 237 238 239 240 241 242 243
}

// ****************************************************************************
//  Method: NetworkManager destructor
//
//  Programmer: Jeremy Meredith
//  Creation:   September 29, 2000
//
//  Modifications:
hrchilds's avatar
hrchilds committed
244 245 246 247 248 249
//    Sean Ahern, Fri May 24 11:12:25 PDT 2002
//    Got rid of deletion of the database cache, since all DBs in there
//    will be deleted by deleting the network cache.
//
//    Hank Childs, Mon Jan  5 16:30:04 PST 2004
//    Don't assume the caches are valid.
hrchilds's avatar
hrchilds committed
250
//
hrchilds's avatar
hrchilds committed
251 252 253
//    Mark C. Miller, Tue Jan  4 10:23:19 PST 2005
//    Added code to delete VisWindow objects
//
hrchilds's avatar
hrchilds committed
254 255 256
//    Hank Childs, Mon Feb 13 23:14:23 PST 2006
//    Delete the DDFs.
//
257 258 259
//    Brad Whitlock, Tue Jun 24 15:41:08 PDT 2008
//    Added plugin managers.
//
260 261 262
//    Hank Childs, Sat Aug 21 14:35:47 PDT 2010
//    Rename DDF to DataBinning.
//
263 264 265 266
//    David Camp, Thu Jan 13 11:15:00 PST 2011
//    Added the call to DeleteInstance of the avtColorTables. 
//    Help debug memory leaks.
//
hrchilds's avatar
hrchilds committed
267
// ****************************************************************************
268

hrchilds's avatar
hrchilds committed
269 270
NetworkManager::~NetworkManager(void)
{
271
    for (size_t i = 0; i < networkCache.size(); i++)
hrchilds's avatar
hrchilds committed
272 273
        if (networkCache[i] != NULL)
            delete networkCache[i];
hrchilds's avatar
hrchilds committed
274 275 276 277

    std::map<int, EngineVisWinInfo>::iterator it;
    for (it = viswinMap.begin(); it != viswinMap.end(); it++)
        delete it->second.viswin;
hrchilds's avatar
hrchilds committed
278

279 280
    for (size_t d = 0 ; d < dataBinnings.size() ; d++)
        delete dataBinnings[d];
281 282 283 284

    delete databasePlugins;
    delete operatorPlugins;
    delete plotPlugins;
285 286

    avtColorTables::Instance()->DeleteInstance();
287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349
}

// ****************************************************************************
// Method: NetworkManager::GetDatabasePluginManager
//
// Purpose: 
//   Return the database plugin manager.
//
// Returns:    the database plugin manager.
//
// Programmer: Brad Whitlock
// Creation:   Tue Jun 24 15:25:41 PDT 2008
//
// Modifications:
//   
// ****************************************************************************

DatabasePluginManager *
NetworkManager::GetDatabasePluginManager() const
{
    return databasePlugins;
}

// ****************************************************************************
// Method: NetworkManager::GetOperatorPluginManager
//
// Purpose: 
//   Return the operator plugin manager.
//
// Returns:    the operator plugin manager.
//
// Programmer: Brad Whitlock
// Creation:   Tue Jun 24 15:26:07 PDT 2008
//
// Modifications:
//   
// ****************************************************************************

OperatorPluginManager *
NetworkManager::GetOperatorPluginManager() const
{
    return operatorPlugins;
}

// ****************************************************************************
// Method: NetworkManager::GetPlotPluginManager
//
// Purpose: 
//   Return the plot plugin manager.
//
// Returns:    the plot plugin manager.
//
// Programmer: Brad Whitlock
// Creation:   Tue Jun 24 15:26:28 PDT 2008
//
// Modifications:
//   
// ****************************************************************************

PlotPluginManager *
NetworkManager::GetPlotPluginManager() const
{
    return plotPlugins;
hrchilds's avatar
hrchilds committed
350 351 352 353 354 355 356 357 358 359 360 361
}

// ****************************************************************************
//  Method: NetworkManager::ClearAllNetworks
//
//  Purpose:
//      Clears out the caches.  Really only used for checking for memory leaks.
//
//  Programmer: Hank Childs
//  Creation:   January 7, 2001
//
//  Modifications:
hrchilds's avatar
hrchilds committed
362 363 364 365 366
//    Sean Ahern, Tue May 14 14:35:42 PDT 2002
//    Retrofitted for my new network management techniques.
//
//    Sean Ahern, Tue May 21 16:09:56 PDT 2002
//    Told the underlying objects to free their data.
hrchilds's avatar
hrchilds committed
367
//
hrchilds's avatar
hrchilds committed
368 369 370
//    Hank Childs, Mon Jan  5 16:30:04 PST 2004
//    Don't clear out the vectors, so this can be consistent with how the
//    ClearNetworksWithDatabase works.
hrchilds's avatar
hrchilds committed
371
//
hrchilds's avatar
hrchilds committed
372 373 374
//    Brad Whitlock, Thu Feb 26 11:58:58 PDT 2004
//    I replaced the commented out cerr with debug3.
//
hrchilds's avatar
hrchilds committed
375 376 377 378 379
//    Hank Childs, Sun Mar  7 16:05:03 PST 2004
//    Delete the networks before the databases, since the networks have
//    references to the databases (and the previous ordering yielded dangling
//    pointers).
//
hrchilds's avatar
hrchilds committed
380 381 382 383
//    Mark C. Miller, Wed Sep  8 17:06:25 PDT 2004
//    Added code to clear the vis window and the list of plots currently in
//    the window.
//
hrchilds's avatar
hrchilds committed
384 385 386 387
//    Mark C. Miller, Mon Sep 13 18:30:26 PDT 2004
//    Protected call to clear viswin and plotsCurrentlyInWindow with test
//    for emptiness
//
hrchilds's avatar
hrchilds committed
388 389 390
//    Mark C. Miller, Tue Jan  4 10:23:19 PST 2005
//    Added code to delete VisWindow objects
//
hrchilds's avatar
hrchilds committed
391 392 393
//    Hank Childs, Thu Mar  2 10:06:33 PST 2006
//    Added support for image based plots.
//
394 395 396
//    Brad Whitlock, Wed Sep  7 13:31:27 PDT 2011
//    Tell the NSM to clear its cache.
//
hrchilds's avatar
hrchilds committed
397
// ****************************************************************************
398

hrchilds's avatar
hrchilds committed
399 400 401
void
NetworkManager::ClearAllNetworks(void)
{
hrchilds's avatar
hrchilds committed
402
    debug3 << "NetworkManager::ClearAllNetworks(void)" << endl;
hrchilds's avatar
hrchilds committed
403

404
    for (size_t i = 0; i < networkCache.size(); i++)
hrchilds's avatar
hrchilds committed
405 406 407 408 409 410
    {
        if (networkCache[i] != NULL)
            delete networkCache[i];
        networkCache[i] = NULL;
    }

411
    for (size_t i = 0; i < databaseCache.size(); i++)
hrchilds's avatar
hrchilds committed
412 413 414 415 416 417
    {
        if (databaseCache[i] != NULL)
            delete databaseCache[i];
        databaseCache[i] = NULL;
    }

418
    for (size_t i = 0 ; i < globalCellCounts.size() ; i++)
hrchilds's avatar
hrchilds committed
419
    {
hrchilds's avatar
hrchilds committed
420
        globalCellCounts[i] = -1;
hrchilds's avatar
hrchilds committed
421
    }
hrchilds's avatar
hrchilds committed
422 423 424 425 426 427

    std::map<int, EngineVisWinInfo>::iterator it;
    for (it = viswinMap.begin(); it != viswinMap.end(); it++)
    {
        it->second.viswin->ClearPlots();
        it->second.plotsCurrentlyInWindow.clear();
hrchilds's avatar
hrchilds committed
428
        it->second.imageBasedPlots.clear();
hrchilds's avatar
hrchilds committed
429 430 431 432
        delete it->second.viswin;
    }
    viswinMap.clear();
    NewVisWindow(0);
433 434

    avtNamedSelectionManager::GetInstance()->ClearCache();
hrchilds's avatar
hrchilds committed
435 436 437 438 439 440 441 442 443 444 445 446 447 448
}

// ****************************************************************************
//  Method: NetworkManager::ClearNetworksWithDatabase
//
//  Purpose:
//      Clears out the caches that reference a specific database.
//
//  Arguments:
//      db      The name of a database.
//
//  Programmer: Hank Childs
//  Creation:   January 5, 2004
//
hrchilds's avatar
hrchilds committed
449 450 451 452
//  Modifications:
//    Brad Whitlock, Thu Feb 26 11:58:58 PDT 2004
//    I replaced the commented out cerr with debug3.
//
hrchilds's avatar
hrchilds committed
453 454 455
//    Mark C. Miller, Wed Sep  8 17:06:25 PDT 2004
//    Added code to clear the vis window and the list of plots currently in
//    the window for the first plot id it encounters in this db
hrchilds's avatar
hrchilds committed
456 457 458 459
//
//    Mark C. Miller, Mon Sep 13 18:30:26 PDT 2004
//    Moved code to delete the networkCache entry to after the clearning
//    of the viswin.
hrchilds's avatar
hrchilds committed
460 461 462
//
//    Mark C. Miller, Tue Jan  4 10:23:19 PST 2005
//    Modified for viswinMap object
463
//
hrchilds's avatar
hrchilds committed
464 465 466 467
// ****************************************************************************
void
NetworkManager::ClearNetworksWithDatabase(const std::string &db)
{
hrchilds's avatar
hrchilds committed
468
    debug3 << "NetworkManager::ClearNetworksWithDatabase()" << endl;
hrchilds's avatar
hrchilds committed
469

470
    //
hrchilds's avatar
hrchilds committed
471 472 473
    // Clear out the networks before the databases.  This is because if we
    // delete the databases first, the networks will have dangling pointers.
    //
474
    for (size_t i = 0; i < networkCache.size(); i++)
hrchilds's avatar
hrchilds committed
475 476 477 478 479 480 481 482
    {
        if (networkCache[i] != NULL)
        {
            NetnodeDB *ndb = networkCache[i]->GetNetDB();
            if (ndb != NULL)
            {
                if (ndb->GetFilename() == db)
                {
hrchilds's avatar
hrchilds committed
483
                    DoneWithNetwork(i);
hrchilds's avatar
hrchilds committed
484 485 486 487
                }
            }
        }
    }
hrchilds's avatar
hrchilds committed
488

489
    // Remove the database from the cache and delete it.
490
    for (size_t i = 0; i < databaseCache.size(); i++)
hrchilds's avatar
hrchilds committed
491 492 493 494 495 496 497 498 499 500
    {
        if (databaseCache[i] != NULL)
        {
            if (databaseCache[i]->GetFilename() == db)
            {
                delete databaseCache[i];
                databaseCache[i] = NULL;
            }
        }
    }
hrchilds's avatar
hrchilds committed
501 502 503
}

// ****************************************************************************
hrchilds's avatar
hrchilds committed
504
//  Method: NetworkManager::GetDBFromCache
hrchilds's avatar
hrchilds committed
505 506
//
//  Purpose:
hrchilds's avatar
hrchilds committed
507 508 509 510 511 512 513 514 515 516 517 518 519
//      Get a DB from the cache.  If it's not there already, create it and
//      return it.
//
//  Programmer: Sean Ahern
//  Creation:   Mon Dec 23 12:28:30 PST 2002
//
//  Modifications:
//
//    Brad Whitlock, Wed May 14 12:56:16 PDT 2003
//    I made it get the metadata and the SIL if the time is greater than 0.
//    This ensures that we pick up the metadata and SIL for the time state
//    that we're intersted in.
//
hrchilds's avatar
hrchilds committed
520 521 522 523
//    Brad Whitlock, Mon Aug 25 16:51:35 PST 2003
//    Changed the return statement to CATCH_RETURN2 so bad things don't happen
//    when we use fake exceptions.
//
hrchilds's avatar
hrchilds committed
524 525 526
//    Mark C. Miller Thu Oct  9 10:59:27 PDT 2003
//    I modified to force GetMetaData and GetSIL if they're invariant 
//
hrchilds's avatar
hrchilds committed
527 528 529
//    Hank Childs, Tue Nov  4 14:19:05 PST 2003
//    Checked in work-around (HACK) to avoid apparent compiler bug on AIX.
//
hrchilds's avatar
hrchilds committed
530 531 532
//    Hank Childs, Mon Jan  5 16:49:28 PST 2004
//    Do not assume that all of the cached databases are non-NULL.
//
hrchilds's avatar
hrchilds committed
533 534 535
//    Hank Childs, Mon Mar  1 08:48:26 PST 2004
//    Send the time to the database factory as well.
//
hrchilds's avatar
hrchilds committed
536 537 538
//    Hank Childs, Mon Mar 22 11:10:43 PST 2004
//    Allow for the DB type to be explicitly specified.
//
hrchilds's avatar
hrchilds committed
539 540 541 542 543
//    Mark C. Miller, Tue Sep 28 19:57:42 PDT 2004
//    Added code to pass avtDatabaseMetaData to LoadBalancer->AddDatabase
//
//    Hank Childs, Sun Feb 27 12:30:17 PST 2005
//    Added avtDatabase argument to AddDatabase.
hrchilds's avatar
hrchilds committed
544
//
hrchilds's avatar
hrchilds committed
545 546 547
//    Hank Childs, Tue Sep 20 13:03:00 PDT 2005
//    Don't assume that filename lengths are bigger than 6 characters ['6616].
//
hrchilds's avatar
hrchilds committed
548 549 550 551 552 553 554
//    Hank Childs, Fri Oct  7 17:22:01 PDT 2005
//    Tell database what it's full name is.  Also, remove some antiquated code
//    that the variable name in the DB info to the empty string.  I believe
//    the code was no longer doing anything useful, besides torpedoing SR
//    and CMFE when used in combination.  The test suite passes without it
//    and there is really no explanation for why we were doing it.
//
hrchilds's avatar
hrchilds committed
555 556
//    Mark C. Miller, Wed Nov 16 10:46:36 PST 2005
//    Changed interface to lb->AddDatabase 
hrchilds's avatar
hrchilds committed
557 558 559 560
//
//    Hank Childs, Thu Jan 11 16:10:12 PST 2007
//    Added argument to DatabaseFactory calls.
//
561 562
//    Mark C. Miller, Thu Jun 14 10:26:37 PDT 2007
//    Added support to treat all databases as time varying
563 564 565
//
//    Mark C. Miller, Wed Aug 22 20:16:59 PDT 2007
//    Added treatAllDBsAsTimeVarying to GetSIL call
566 567 568 569
//
//    Hank Childs, Fri Feb  1 15:45:28 PST 2008
//    Added argument for loading database plugins.
//
570 571
//    Mark C. Miller, Tue Jun 10 22:36:25 PDT 2008
//    Added support for ignoring bad extents from dbs.
572 573 574 575
//
//    Brad Whitlock, Tue Jun 24 15:45:43 PDT 2008
//    Pass the database plugin manager to the factory.
//
hrchilds's avatar
hrchilds committed
576
// ****************************************************************************
hrchilds's avatar
hrchilds committed
577

hrchilds's avatar
hrchilds committed
578
NetnodeDB *
579
NetworkManager::GetDBFromCache(const std::string &filename, int time,
580
    const char *format, bool treatAllDBsAsTimeVarying, 
581
    bool fileMayHaveUnloadedPlugin, bool ignoreExtents)
hrchilds's avatar
hrchilds committed
582 583 584 585
{
    // If we don't have a load balancer, we're dead.
    if (loadBalancer == NULL)
    {
hrchilds's avatar
hrchilds committed
586 587
        debug1 << "Internal error: A load balancer was never registered." 
               << endl;
hrchilds's avatar
hrchilds committed
588 589 590
        EXCEPTION0(ImproperUseException);
    }

hrchilds's avatar
hrchilds committed
591 592
    // Find a matching database 
    NetnodeDB* cachedDB = NULL;
593
    for (size_t i = 0; i < databaseCache.size(); i++)
hrchilds's avatar
hrchilds committed
594
    {
hrchilds's avatar
hrchilds committed
595 596
        if (databaseCache[i] == NULL)
            continue;
hrchilds's avatar
hrchilds committed
597
        if (databaseCache[i]->GetFilename() == filename)
hrchilds's avatar
hrchilds committed
598
        {
hrchilds's avatar
hrchilds committed
599 600
            cachedDB = databaseCache[i];
            break;
hrchilds's avatar
hrchilds committed
601
        }
hrchilds's avatar
hrchilds committed
602 603
    }

604 605 606 607
    // use of these makes calls to GetMetaData a little clearer
    const bool forceReadAllCyclesAndTimes = false;
    const bool forceReadThisCycleAndTime  = false;

hrchilds's avatar
hrchilds committed
608 609 610 611 612
    // got a match
    if (cachedDB != NULL)
    {
        // even if we found the DB in the cache,
        // we need to update the metadata if its time-varying
613
        if (treatAllDBsAsTimeVarying ||
614
            !cachedDB->GetDB()->MetaDataIsInvariant() ||
hrchilds's avatar
hrchilds committed
615
            !cachedDB->GetDB()->SILIsInvariant())
hrchilds's avatar
hrchilds committed
616
        {
617
            cachedDB->GetDB()->GetMetaData(time,
618 619 620
                                           forceReadAllCyclesAndTimes,
                                           forceReadThisCycleAndTime,
                                           treatAllDBsAsTimeVarying);
621
            cachedDB->GetDB()->GetSIL(time, treatAllDBsAsTimeVarying);
hrchilds's avatar
hrchilds committed
622
        }
hrchilds's avatar
hrchilds committed
623 624

        return cachedDB;
hrchilds's avatar
hrchilds committed
625 626 627 628 629 630
    }

    // No match.  Load a new DB.
    debug3 << "Loading new database" << endl;
    TRY
    {
631
        if (fileMayHaveUnloadedPlugin)
632 633
             GetDatabasePluginManager()->LoadPluginsNow();

hrchilds's avatar
hrchilds committed
634 635 636
        avtDatabase *db = NULL;
        NetnodeDB *netDB = NULL;
        const char *filename_c = filename.c_str();
hrchilds's avatar
hrchilds committed
637
        std::vector<std::string> plugins;  // unused
hrchilds's avatar
hrchilds committed
638 639
        if (filename.length() >= 6 &&
            filename.substr(filename.length() - 6) == ".visit")
640 641 642 643 644 645 646 647
        {
            db = avtDatabaseFactory::VisitFile(
                GetDatabasePluginManager(),
                filename_c,
                time,
                plugins,
                format);
        }
hrchilds's avatar
hrchilds committed
648
        else
649 650 651 652 653 654 655 656 657 658
        {
            db = avtDatabaseFactory::FileList(
                GetDatabasePluginManager(),
                &filename_c,
                1,
                time,
                plugins, 
                format);
        }

hrchilds's avatar
hrchilds committed
659
        db->SetFullDBName(filename);
660
        db->SetIgnoreExtents(ignoreExtents);
hrchilds's avatar
hrchilds committed
661 662

        // If we want to open the file at a later timestep, get the
hrchilds's avatar
hrchilds committed
663 664
        // SIL so that it contains the right data.
        if ((time > 0) ||
665
            treatAllDBsAsTimeVarying ||
hrchilds's avatar
hrchilds committed
666 667
            (!db->MetaDataIsInvariant()) ||
            (!db->SILIsInvariant()))
hrchilds's avatar
hrchilds committed
668 669 670
        {
            debug2 << "NetworkManager::AddDB: We were instructed to open "
                   << filename.c_str() << " at timestate=" << time
hrchilds's avatar
hrchilds committed
671
                   << " so we're reading the SIL early."
hrchilds's avatar
hrchilds committed
672
                   << endl;
673
            db->GetMetaData(time,
674 675 676
                            forceReadAllCyclesAndTimes,
                            forceReadThisCycleAndTime,
                            treatAllDBsAsTimeVarying);
677
            db->GetSIL(time, treatAllDBsAsTimeVarying);
hrchilds's avatar
hrchilds committed
678 679 680 681 682 683
        }

        netDB = new NetnodeDB(db);
        databaseCache.push_back(netDB);

        netDB->SetDBInfo(filename, "", time);
hrchilds's avatar
hrchilds committed
684
        loadBalancer->AddDatabase(filename, db, time);
hrchilds's avatar
hrchilds committed
685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720
      
        // The code here should be:
        // CATCH_RETURN2(1, netDB);
        //
        // However: this is crashing on AIX when you do a re-open (see '3984).
        // After investigating this bug, I have determined the following:
        // (1) A valid pointer is being returned from this function.  However
        // the calling function receives NULL (leading to a seg fault).
        // (2) When compiling with fake exceptions and running under purify
        // on the Suns, there are no warnings.
        // (3) The catch return unwinds into 3 calls.  If you put a print
        // statement right before the last one (the return), everything is
        // returned correctly and there is no seg fault.
        //
        // I'm always very leery to call something a compiler bug.  However,
        // on different platforms there is no problem.  I think it is even
        // more fishy that things work if there is a print statement and don't
        // if there is not.  It certainly sounds like a compiler bug.
        //
        // Below is code that "makes this work".  It is essentially a
        // CATCH_RETURN2 with a print statement in the middle.  If CATCH_RETURN
        // were to ever change, this coding would have to change as well
        // (thus, I have introduced a maintenance problem).  As such, I put
        // a comment in VisItException.h referencing this comment.  If the
        // coding below can ever be replaced with an everyday CATCH_RETURN,
        // the comment in VisItException.h should be removed.
        //
#ifdef FAKE_EXCEPTIONS
           exception_delete(exception_caught); 
           jump_stack_top -= (1); 
           debug5 << "This debug statement (from the NetworkManager) is "
                  << "being issued to get around a compiler bug." << endl;
           return (netDB); 
#else
           return netDB;
#endif
hrchilds's avatar
hrchilds committed
721 722 723 724 725 726 727 728 729 730 731 732 733 734 735
    }
    CATCH(DatabaseException)
    {
        debug1 << "ERROR - could not create database " << filename.c_str()
            << endl;
        RETHROW;
    }
    ENDTRY
}

// ****************************************************************************
//  Method: NetworkManager::StartNetwork
//
//  Purpose:
//      Start building the pipeline.
hrchilds's avatar
hrchilds committed
736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781
//
//  Programmer: Jeremy Meredith
//  Creation:   September 29, 2000
//
//  Modifications:
//
//    Jeremy Meredith, Thu Nov  2 14:07:34 PST 2000
//    Changed the file initialization.  Apparently substr() may not portably 
//    be a const method of strings....
//
//    Hank Childs, Tue Jun 12 20:01:55 PDT 2001
//    Incorporate data specifications and the load balancer.
//
//    Jeremy Meredith, Thu Jul 26 12:08:42 PDT 2001
//    Fixed a bug passing the data file name to the load balancer.
//
//    Hank Childs, Thu Jul 26 20:17:56 PDT 2001
//    Re-enable caching by checking data specifications.
//
//    Brad Whitlock, Mon Oct 22 18:33:37 PST 2001
//    Changed the exception keywords to macros.
//
//    Jeremy Meredith, Fri Nov  9 10:24:23 PST 2001
//    Replaced netPool with netMRU, and simply made it release its
//    data, not delete the whole network. Made is only set the data
//    specification, not the pipeline specification yet.
//
//    Hank Childs, Fri Dec 14 17:39:36 PST 2001
//    Add support for using a more compact version of sil attributes.
//
//    Sean Ahern, Mon May 13 17:24:25 PDT 2002
//    Refactored to start support for a tree-based network scheme.
//
//    Sean Ahern, Tue May 28 09:47:21 PDT 2002
//    Added tree-based networks.
//
//    Brad Whitlock, Fri Jun 28 14:03:49 PST 2002
//    I renamed Network to DataNetwork to make the Windows compiler happy.
//
//    Hank Childs, Thu Sep 12 11:11:16 PDT 2002
//    There was an issue with multiple networks sharing NetnodeDBs, so I
//    made them all have their own reference.
//
//    Jeremy Meredith, Thu Oct 24 16:03:39 PDT 2002
//    Added material options.
//
hrchilds's avatar
hrchilds committed
782 783 784 785 786 787
//    Sean Ahern, Mon Dec 23 12:38:17 PST 2002
//    Changed this dramatically.  Broke out a good portion of the logic
//    into GetDBFromCache.
//
//    Sean Ahern, Sat Mar  8 00:59:16 America/Los_Angeles 2003
//    Cleaned up unnecessarily complicated code.
hrchilds's avatar
hrchilds committed
788
//
hrchilds's avatar
hrchilds committed
789 790 791
//    Jeremy Meredith, Wed Jul 30 10:45:45 PDT 2003
//    Added the check for requiring full connectivity.
//
hrchilds's avatar
hrchilds committed
792 793 794
//    Jeremy Meredith, Fri Sep  5 15:23:13 PDT 2003
//    Added a flag for the new MIR algorithm.
//
hrchilds's avatar
hrchilds committed
795 796 797
//    Jeremy Meredith, Thu Oct 30 16:09:32 PST 2003
//    Added code to make sure varLeaves was non-empty before accessing it.
//
hrchilds's avatar
hrchilds committed
798 799 800
//    Jeremy Meredith, Fri Oct 31 13:05:26 PST 2003
//    Made the error message for no-real-variables more informative.
//
hrchilds's avatar
hrchilds committed
801 802 803 804
//    Hank Childs, Sun Mar  7 16:05:03 PST 2004
//    Do not make the database be a node for the network.  Also add the EEF
//    as a node.
//
hrchilds's avatar
hrchilds committed
805 806 807
//    Jeremy Meredith, Tue Mar 23 14:40:20 PST 2004
//    Added the file format as an argument.  Made use of it in GetDBFromCache.
//
hrchilds's avatar
hrchilds committed
808 809 810
//    Mark C. Miller, Thu Aug 19 10:28:45 PDT 2004
//    Removed #ifdef 0 code block
//
hrchilds's avatar
hrchilds committed
811 812 813 814 815 816
//    Jeremy Meredith, Tue Dec 14 14:02:35 PST 2004
//    The code to get the real variable name was a duplicate of two
//    other spots in VisIt, but this one was out of date.  I
//    refactored the best one into ParsingExprList::GetRealVariable
//    and made this one point to it.
//
hrchilds's avatar
hrchilds committed
817 818 819 820
//    Brad Whitlock, Tue Feb 22 12:13:44 PDT 2005
//    Made format argument come first. I also made it possible for an empty
//    format string to make the code behave as through no format was specified.
//
hrchilds's avatar
hrchilds committed
821 822 823
//    Hank Childs, Wed Jul 13 10:38:02 PDT 2005
//    Added call to SetTime to prevent UMR when exporting databases.
//
hrchilds's avatar
hrchilds committed
824 825 826
//    Hank Childs, Wed Aug 17 21:55:55 PDT 2005
//    Pass along the material options to simplify the interface reconstruction.
//
hrchilds's avatar
hrchilds committed
827 828 829
//    Jeremy Meredith, Thu Aug 18 17:54:51 PDT 2005
//    Added a new isovolume algorithm, with adjustable VF cutoff.
//
hrchilds's avatar
hrchilds committed
830 831
//    Mark C. Miller, Wed Nov 16 10:46:36 PST 2005
//    Added mesh management attributes 
hrchilds's avatar
hrchilds committed
832 833 834
//
//    Mark C. Miller, Sun Dec  3 12:20:11 PST 2006
//    Added setting of flatness tolerance to data spec for CSG stuff
hrchilds's avatar
hrchilds committed
835 836 837 838 839
//
//    Hank Childs, Wed Jan 17 10:42:21 PST 2007
//    Tell the output of the expression evaluator filter that it is *not*
//    transient, since picks sometimes like to look at these.
//
840 841
//    Mark C. Miller, Thu Jun 14 10:26:37 PDT 2007
//    Added support to treat all databases as time varying
842 843 844
//
//    Mark C. Miller, Wed Aug 22 20:16:59 PDT 2007
//    Added treatAllDBsAsTimeVarying to GetSIL call.
845 846 847 848
//
//    Hank Childs, Fri Feb  1 15:48:01 PST 2008
//    Add new Boolean argument to GetDBFromCache. 
//
849 850
//    Mark C. Miller, Tue Jun 10 22:36:25 PDT 2008
//    Added support for ignoring bad extents from dbs.
851 852 853 854
//
//    Hank Childs, Fri Jan 16 15:30:15 PST 2009
//    Store the variable with the working network.
//
855 856 857
//    John C. Anderson, Thu Jan 15 10:20:20 2009
//    Added annealing time for Discrete MIR.
//
858 859 860
//    Jeremy Meredith, Fri Feb 13 11:22:39 EST 2009
//    Added MIR iteration capability.
//
861 862 863
//    Brad Whitlock, Mon Aug 22 13:34:36 PDT 2011
//    I added named selection creation here instead of at the end of the pipeline.
//
864 865 866 867 868 869
//    Brad Whitlock, Thu Sep  1 11:07:55 PDT 2011
//    Stash the selection name into the data request in case we need it 
//    when creating other pipelines based on this one. For instance, the 
//    avtExecuteThenTimeLoopFilter creates its own pipeline from the original
//    source so we need the selection name there to apply it.
//
hrchilds's avatar
hrchilds committed
870
// ****************************************************************************
hrchilds's avatar
hrchilds committed
871

hrchilds's avatar
hrchilds committed
872
void
873 874 875
NetworkManager::StartNetwork(const std::string &format,
                             const std::string &filename,
                             const std::string &var,
hrchilds's avatar
hrchilds committed
876 877
                             int time,
                             const CompactSILRestrictionAttributes &atts,
hrchilds's avatar
hrchilds committed
878
                             const MaterialAttributes &matopts,
879
                             const MeshManagementAttributes &meshopts,
880
                             bool treatAllDBsAsTimeVarying,
881 882
                             bool ignoreExtents,
                             const std::string &selName)
hrchilds's avatar
hrchilds committed
883
{
hrchilds's avatar
hrchilds committed
884 885
    // If the variable is an expression, we need to find a "real" variable
    // name to work with.
886
    std::string leaf = ParsingExprList::GetRealVariable(var);
hrchilds's avatar
hrchilds committed
887

hrchilds's avatar
hrchilds committed
888 889 890 891 892
    // Make empty strings behave as though no format was specified.
    const char *defaultFormat = 0;
    if(format.size() > 0)
        defaultFormat = format.c_str();

hrchilds's avatar
hrchilds committed
893 894
    // Start up the DataNetwork and add the database to it.
    workingNet = new DataNetwork;
895
    bool fileMayHaveUnloadedPlugin = false;
896
    NetnodeDB *netDB = GetDBFromCache(filename, time, defaultFormat,
897
                                      treatAllDBsAsTimeVarying, 
898 899
                                      fileMayHaveUnloadedPlugin,
                                      ignoreExtents);
hrchilds's avatar
hrchilds committed
900
    workingNet->SetNetDB(netDB);
901
    workingNet->SetVariable(leaf);
hrchilds's avatar
hrchilds committed
902
    netDB->SetDBInfo(filename, leaf, time);
903
    Netnode *input = netDB;
hrchilds's avatar
hrchilds committed
904

905
    // Add ExpressionEvaluatorFilter to handle expressions that come up the pipe.
hrchilds's avatar
hrchilds committed
906 907
    avtExpressionEvaluatorFilter *f = new avtExpressionEvaluatorFilter();
    NetnodeFilter *filt = new NetnodeFilter(f, "ExpressionEvaluator");
908
    filt->GetInputNodes().push_back(input);
hrchilds's avatar
hrchilds committed
909
    f->GetOutput()->SetTransientStatus(false);
hrchilds's avatar
hrchilds committed
910

911 912 913 914 915 916 917 918 919 920 921 922 923 924
    input = filt;
    if(!selName.empty())
    {
        // Add the EEF to the network.
        workingNet->AddNode(filt);

        avtNamedSelectionFilter *f = new avtNamedSelectionFilter();
        f->SetSelectionName(selName);
        filt = new NetnodeFilter(f, "NamedSelection");
        filt->GetInputNodes().push_back(input);
        workingNet->SetSelectionName(selName);
    }

    // Push the last filter onto the working list.
hrchilds's avatar
hrchilds committed
925
    workingNetnodeList.push_back(filt);
hrchilds's avatar
hrchilds committed
926

hrchilds's avatar
hrchilds committed
927
    workingNet->AddNode(filt);
hrchilds's avatar
hrchilds committed
928
    // Push the variable name onto the name stack.
hrchilds's avatar
hrchilds committed
929
    nameStack.push_back(var);
930
    debug4 << "NetworkManager::AddDB: Adding " << var.c_str()
hrchilds's avatar
hrchilds committed
931 932 933
           << " to the name stack" << endl;

    // Set up the data spec.
hrchilds's avatar
hrchilds committed
934
    avtSILRestriction_p silr =
935
        new avtSILRestriction(workingNet->GetNetDB()->GetDB()->
936
            GetSIL(time, treatAllDBsAsTimeVarying), atts);
hrchilds's avatar
hrchilds committed
937
    avtDataRequest *dataRequest = new avtDataRequest(var.c_str(), time, silr);
hrchilds's avatar
hrchilds committed
938 939

    // Set up some options from the data specification
hrchilds's avatar
hrchilds committed
940 941 942 943 944
    dataRequest->SetNeedMixedVariableReconstruction(matopts.GetForceMIR());
    dataRequest->SetNeedSmoothMaterialInterfaces(matopts.GetSmoothing());
    dataRequest->SetNeedCleanZonesOnly(matopts.GetCleanZonesOnly());
    dataRequest->SetNeedValidFaceConnectivity(matopts.GetNeedValidConnectivity());
    dataRequest->SetMIRAlgorithm(matopts.GetAlgorithm());
945 946 947
    dataRequest->SetMIRNumIterations(matopts.GetIterationEnabled() ?
                                     matopts.GetNumIterations() : 0);
    dataRequest->SetMIRIterationDamping(matopts.GetIterationDamping());
hrchilds's avatar
hrchilds committed
948 949 950
    dataRequest->SetSimplifyHeavilyMixedZones(matopts.GetSimplifyHeavilyMixedZones());
    dataRequest->SetMaxMaterialsPerZone(matopts.GetMaxMaterialsPerZone());
    dataRequest->SetIsovolumeMIRVF(matopts.GetIsoVolumeFraction());
951
    dataRequest->SetAnnealingTime(matopts.GetAnnealingTime());
hrchilds's avatar
hrchilds committed
952 953 954 955 956
    dataRequest->SetDiscTol(meshopts.GetDiscretizationTolerance()[0]);
    dataRequest->SetFlatTol(meshopts.GetDiscretizationTolerance()[1]);
    dataRequest->SetDiscMode(meshopts.GetDiscretizationMode());
    dataRequest->SetDiscBoundaryOnly(meshopts.GetDiscretizeBoundaryOnly());
    dataRequest->SetPassNativeCSG(meshopts.GetPassNativeCSG());
957
    dataRequest->SetSelectionName(selName);
hrchilds's avatar
hrchilds committed
958 959
    workingNet->SetDataSpec(dataRequest);
    workingNet->SetTime(dataRequest->GetTimestep());
hrchilds's avatar
hrchilds committed
960 961 962 963 964 965 966 967 968 969 970 971 972 973 974

    // The plot starts out as NULL.
    workingNet->SetPlot(NULL);
}

// ****************************************************************************
// Method: NetworkManager::DefineDB
//
// Purpose: 
//   Defines a new database with the specified name and list of files.
//
// Arguments:
//    dbName   : The database name.
//    files    : The list of files that make up the database.
//    time     : The timestep that we want to examine.
hrchilds's avatar
hrchilds committed
975
//    format   : The file format type of the DB.
hrchilds's avatar
hrchilds committed
976 977 978 979 980 981 982 983 984 985
//   
// Programmer: Brad Whitlock
// Creation:   Tue Mar 25 13:41:08 PST 2003
//
// Modifications:
//   Brad Whitlock, Wed May 14 12:56:16 PDT 2003
//   I made it get the metadata and the SIL if the time is greater than 0.
//   This ensures that we pick up the metadata and SIL for the time state
//   that we're intersted in.
//
hrchilds's avatar
hrchilds committed
986 987 988
//   Hank Childs, Mon Jan  5 16:49:28 PST 2004
//   Account for databases in the cached lists that have been cleared.
//
hrchilds's avatar
hrchilds committed
989 990 991
//   Hank Childs, Mon Mar  1 08:48:26 PST 2004
//   Send the time to the database factory as well.
//
hrchilds's avatar
hrchilds committed
992 993
//   Hank Childs, Mon Mar 22 11:10:43 PST 2004
//   Send the file format type to the database factory.
hrchilds's avatar
hrchilds committed
994 995 996 997
//
//   Mark C. Miller, Tue Sep 28 19:57:42 PDT 2004
//   Added code to pass avtDatabaseMetaData to LoadBalancer->AddDatabase
//
hrchilds's avatar
hrchilds committed
998 999 1000 1001
//   Brad Whitlock, Tue Feb 22 12:17:27 PDT 2005
//   I changed the code so it an empty format string behaves as though no
//   format was specified.
//
hrchilds's avatar
hrchilds committed
1002 1003 1004
//   Hank Childs, Sun Feb 27 12:30:17 PST 2005
//   Added avtDatabase argument to AddDatabase.
//
hrchilds's avatar
hrchilds committed
1005 1006 1007
//   Hank Childs, Fri Oct  7 09:30:57 PDT 2005
//   Tell the database what its full name is.
//
hrchilds's avatar
hrchilds committed
1008 1009 1010 1011 1012 1013
//   Mark C. Miller, Wed Nov 16 10:46:36 PST 2005
//   Changed interface to lb->AddDatabase 
//
//   Hank Childs, Thu Jan 11 16:10:12 PST 2007
//   Added argument to DatabaseFactory calls.
//
1014 1015 1016
//   Brad Whitlock, Tue Jun 24 15:48:27 PDT 2008
//   Pass the database plugin manager to the database factory.
//
hrchilds's avatar
hrchilds committed
1017
// ****************************************************************************
hrchilds's avatar
hrchilds committed
1018

hrchilds's avatar
hrchilds committed
1019
void
1020 1021
NetworkManager::DefineDB(const std::string &dbName, const std::string &dbPath,
    const stringVector &files, int time, const std::string &format)
hrchilds's avatar
hrchilds committed
1022 1023 1024 1025 1026 1027 1028 1029 1030 1031 1032 1033 1034 1035
{
    //
    // We gotta have a load balancer
    //
    if (loadBalancer == NULL)
    {
        debug1 << "Internal error: A load balancer was never registered."
               << endl;
        EXCEPTION0(ImproperUseException);
    }

    // Prepend the path to the files.
    stringVector filesWithPath;
    filesWithPath.reserve(files.size());
1036
    for(size_t i = 0; i < files.size(); ++i)
hrchilds's avatar
hrchilds committed
1037 1038 1039 1040 1041 1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060
        filesWithPath.push_back(dbPath + files[i]);

    //
    // If the database is in the virtualDatabases map, then check the
    // new list of files against the old list of files.
    //
    StringVectorMap::const_iterator it = virtualDatabases.find(dbName);
    if(it != virtualDatabases.end())
    {
        // If the list of files is the same, return because we don't
        // need to change the list of files. If the list of files differ
        // then replace the list of files.
        if(filesWithPath == it->second)
        {
            debug3 << "NetworkManager::DefineDB: No need to change "
                   << "virtual database definition for "
                   << dbName.c_str() << endl;

            // The definition is the same so we don't need to change it
            // but we have to make sure that it exists in the database
            // cache or else we'll have to define it. It's possible that
            // a clear cache rpc could have purged it from the database
            // cache.
            bool found = false;
1061
            for(size_t i = 0; i < databaseCache.size() && !found; ++i)
hrchilds's avatar
hrchilds committed
1062 1063 1064
            {
                if (databaseCache[i] == NULL)
                    continue;
hrchilds's avatar
hrchilds committed
1065
                found = (databaseCache[i]->GetFilename() == dbName);
hrchilds's avatar
hrchilds committed
1066
            }
hrchilds's avatar
hrchilds committed
1067 1068 1069 1070 1071 1072 1073 1074 1075 1076 1077 1078 1079 1080 1081 1082 1083 1084

            // If the database was in the database cache, we can return.
            if(found)
                return;
        }
        else
        {
            debug3 << "NetworkManager::DefineDB: Changed "
                   << "virtual database definition for "
                   << dbName.c_str() << endl;
            virtualDatabases[dbName] = filesWithPath;
        }
    }

    //
    // Remove from the database cache databases that have the same name.
    //
    std::vector<NetnodeDB*> databaseBadMatches;
1085
    for (size_t i = 0; i < databaseCache.size(); i++)
hrchilds's avatar
hrchilds committed
1086
    {
hrchilds's avatar
hrchilds committed
1087 1088
        if (databaseCache[i] == NULL)
            continue;
hrchilds's avatar
hrchilds committed
1089 1090 1091 1092 1093 1094 1095 1096 1097 1098 1099 1100
        if (databaseCache[i]->GetFilename() != dbName)
            databaseBadMatches.push_back(databaseCache[i]);
    }
    databaseCache = databaseBadMatches;

    const char **names = 0;
    TRY
    {
        // Open the database
        avtDatabase *db = NULL;
        const char *dbName_c = dbName.c_str();

hrchilds's avatar
hrchilds committed
1101 1102 1103 1104 1105
        // Make empty strings behave as though no format was specified.
        const char *defaultFormat = 0;
        if(format.size() > 0)
            defaultFormat = format.c_str();

hrchilds's avatar
hrchilds committed
1106
        std::vector<std::string> plugins;  // unused
hrchilds's avatar
hrchilds committed
1107 1108 1109 1110 1111 1112 1113 1114
        if (filesWithPath.size() > 0)
        {
            // Make an array of pointers that we can pass to the database
            // factory so we can create a database based on a virtual
            // file.
            names = new const char *[filesWithPath.size()];
            for(int i = 0; i < filesWithPath.size(); ++i)
                names[i] = filesWithPath[i].c_str();
1115 1116
            db = avtDatabaseFactory::FileList(GetDatabasePluginManager(),
                  names, filesWithPath.size(), time, plugins, defaultFormat);
hrchilds's avatar
hrchilds committed
1117 1118 1119 1120 1121 1122 1123 1124 1125
            delete [] names;
            names = 0;

            // Add the virtual database to the virtual database map.
            virtualDatabases[dbName] = filesWithPath;
            debug3 << "NetworkManager::DefineDB: Added new virtual database "
                   << "definition for " << dbName.c_str() << endl;
        }
        else if (dbName.substr(dbName.length() - 6) == ".visit")
1126 1127 1128 1129
        {
            db = avtDatabaseFactory::VisitFile(GetDatabasePluginManager(),
                dbName_c, time, plugins, defaultFormat);
        }
hrchilds's avatar
hrchilds committed
1130
        else
1131 1132 1133 1134
        {
            db = avtDatabaseFactory::FileList(GetDatabasePluginManager(),
                &dbName_c, 1, time, plugins, defaultFormat);
        }
hrchilds's avatar
hrchilds committed
1135
        db->SetFullDBName(dbName);
hrchilds's avatar
hrchilds committed
1136 1137 1138

        // If we want to open the file at a later timestep, get the
        // metadata and the SIL so that it contains the right data.
hrchilds's avatar
hrchilds committed
1139 1140 1141
        if ((time > 0) ||
            (!db->MetaDataIsInvariant()) || 
            (!db->SILIsInvariant()))
hrchilds's avatar
hrchilds committed
1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153
        {
            debug2 << "NetworkManager::DefineDB: We were instructed to define "
                       << dbName.c_str() << " at timestate=" << time
                       << " so we're reading the MetaData and SIL early."
                       << endl;
            db->GetMetaData(time);
            db->GetSIL(time);
        }

        // Create a netnode using the opened database and stick it in
        // the database cache.
        NetnodeDB *netDB = new NetnodeDB(db);
hrchilds's avatar
hrchilds committed
1154
        netDB->SetDBInfo(dbName, "", time);
hrchilds's avatar
hrchilds committed
1155 1156 1157
        databaseCache.push_back(netDB);

        // Add the database to the load balancer.
hrchilds's avatar
hrchilds committed
1158
        loadBalancer->AddDatabase(dbName, db, time);
hrchilds's avatar
hrchilds committed
1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 1185 1186 1187 1188 1189 1190 1191 1192 1193 1194 1195 1196 1197
    }
    CATCH(DatabaseException)
    {
        delete [] names;
        debug1 << "ERROR - could not define database " << dbName.c_str()
               << endl;
        RETHROW;
    }
    ENDTRY 
}

// ****************************************************************************
//  Method: NetworkManager::AddFilter
//
//  Purpose:
//      Create a new filter and add it to the network
//
//  Programmer: Jeremy Meredith
//  Creation:   February 28, 2001
//
//  Modifications:
//    Jeremy Meredith, Sun Mar  4 16:58:13 PST 2001
//    Made it use (a new) PluginManager to create the operator.
//
//    Hank Childs, Thu Apr 12 15:34:22 PDT 2001
//    Added check to make sure that there were enough filters to compare to.
//
//    Jeremy Meredith, Thu Jul 26 03:39:21 PDT 2001
//    Added support for the new (real) operator plugin manager.
//
//    Jeremy Meredith, Fri Sep 28 13:32:56 PDT 2001
//    Made use of plugin manager keying off of ids instead of indexes.
//
//    Jeremy Meredith, Fri Nov  9 10:25:46 PST 2001
//    Added an error check, and renamed type to filtertype.
//
//    Sean Ahern, Tue May 28 08:52:48 PDT 2002
//    Added binary operators.
//
1198 1199 1200
//    Brad Whitlock, Tue Jun 24 16:08:45 PDT 2008
//    Changed how the operator plugin manager is accessed.
//
hrchilds's avatar
hrchilds committed
1201 1202
// ****************************************************************************
void
1203
NetworkManager::AddFilter(const std::string &filtertype,
hrchilds's avatar
hrchilds committed
1204 1205 1206 1207 1208 1209 1210 1211 1212 1213
                          const AttributeGroup *atts,
                          const unsigned int nInputs)
{
    // Check that we have a network to work on.
    if (workingNet == NULL)
    {
        debug1 << "Adding a filter to a non-existent network." << endl;
        EXCEPTION0(ImproperUseException);
    }

1214
    avtPluginFilter *f = GetOperatorPluginManager()->
hrchilds's avatar
hrchilds committed
1215 1216 1217 1218 1219 1220 1221 1222 1223 1224 1225
        GetEnginePluginInfo(filtertype)->AllocAvtPluginFilter();
    f->SetAtts(atts);
    if (nInputs != 1)
    {
        debug1 << "Unable to handle filters with nInputs (" << nInputs << ") > 1"
            << endl;
        EXCEPTION0(ImproperUseException);
    }

    NetnodeFilter *filt = new NetnodeFilter(f, filtertype);
    std::vector<Netnode*> &filtInputs = filt->GetInputNodes();
fogal1's avatar
fogal1 committed
1226
    for (unsigned int i = 0; i < nInputs; i++)
hrchilds's avatar
hrchilds committed
1227
    {
hrchilds's avatar
hrchilds committed
1228
        // Pull a node off the working list and push it onto the filter
hrchilds's avatar
hrchilds committed
1229 1230 1231 1232 1233 1234 1235 1236 1237 1238 1239 1240 1241 1242 1243 1244 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255 1256 1257 1258 1259 1260 1261 1262 1263 1264 1265
        // inputs.
        Netnode *n = workingNetnodeList.back();
        workingNetnodeList.pop_back();
        filtInputs.push_back(n);
    }
    // Push the filter onto the working list.
    workingNetnodeList.push_back(filt);
    workingNet->AddNode(filt);
}

// ****************************************************************************
//  Method:  NetworkManager::MakePlot
//
//  Purpose:
//    Create the plot for the network.
//
//  Programmer:  Jeremy Meredith
//  Creation:    March  1, 2001
//
//  Modifications:
//    Jeremy Meredith, Sun Mar  4 16:58:13 PST 2001
//    Made it use (a new) PluginManager to create the plot.
//
//    Jeremy Meredith, Tue Mar 20 12:05:50 PST 2001
//    Made it use the new plugin manager.  For now we do a reverse lookup
//    of the name to an index.
//
//    Jeremy Meredith, Fri Sep 28 13:32:56 PDT 2001
//    Made use of plugin manager keying off of ids instead of indexes.
//
//    Jeremy Meredith, Fri Nov  9 10:26:18 PST 2001
//    Allowed it to reuse plots so we can cache whole networks again.
//
//    Jeremy Meredith, Fri Feb  8 10:37:59 PST 2002
//    Made sure the total number of filters was the same when matching
//    networks (since when a plot is applied, no more filter can be).
//
hrchilds's avatar
hrchilds committed
1266 1267 1268
//    Eric Brugger, Fri Mar 19 15:25:20 PST 2004
//    Modified the rpc to set the data extents in the engine.
//
hrchilds's avatar
hrchilds committed
1269 1270 1271
//    Brad Whitlock, Wed Mar 21 23:00:00 PST 2007
//    Added plotName.
//
1272 1273 1274
//    Gunther H. Weber, Thu Apr 12 10:52:36 PDT 2007
//    Add filter to beginning of pipeline if necessary
//
1275 1276 1277
//    Brad Whitlock, Tue Jun 24 16:09:19 PDT 2008
//    Changed how the plot plugin manager is accessed.
//
hrchilds's avatar
hrchilds committed
1278
// ****************************************************************************
1279

hrchilds's avatar
hrchilds committed
1280
void
1281 1282
NetworkManager::MakePlot(const std::string &plotName, const std::string &pluginID, 
    const AttributeGroup *atts, const std::vector<double> &dataExtents)
hrchilds's avatar
hrchilds committed
1283 1284 1285 1286 1287 1288 1289 1290 1291 1292 1293 1294 1295 1296 1297 1298
{
    if (workingNet == NULL)
    {
        debug1 << "Adding a plot to a non-existent network." << endl;
        EXCEPTION0(ImproperUseException);
    }

    if (*workingNet->GetPlot() != NULL)
    {
        debug1 << "Trying to stipulate plot type more than once." << endl;
        EXCEPTION0(ImproperUseException);
    }

    // Checking to see if the network has been built successfully.
    if (workingNetnodeList.size() != 1)
    {
hrchilds's avatar
hrchilds committed
1299
        debug1 << "Network building still in progress.  Filter required to "
hrchilds's avatar
hrchilds committed
1300
               << "absorb " << workingNetnodeList.size() << " nodes."  << endl;
hrchilds's avatar
hrchilds committed
1301 1302 1303
        EXCEPTION0(ImproperUseException);
    }

1304 1305
    avtPlot *p = GetPlotPluginManager()->
        GetEnginePluginInfo(pluginID)->AllocAvtPlot();
hrchilds's avatar
hrchilds committed
1306

1307 1308
    // Check, whether plot wants to place a filter at the beginning of
    // the pipeline
1309 1310 1311 1312
    if (avtFilter *f = p->GetFilterForTopOfPipeline()) 
    {
        debug4 << "NetworkManager::MakePlot(): Inserting filter on top of "
               << "pipeline." << endl;
1313 1314 1315 1316 1317 1318 1319 1320 1321 1322 1323 1324 1325 1326 1327
        NetnodeFilter *filt = new NetnodeFilter(f, "InsertedPlotFilter");
        //f->GetOutput()->SetTransientStatus(true);

        if (workingNet->GetNodeList().size() > 1)
            workingNet->AddFilterNodeAfterExpressionEvaluator(filt);
        else
        {
            std::vector<Netnode*> &filtInputs = filt->GetInputNodes();
            Netnode *n = workingNetnodeList.back();
            workingNetnodeList.pop_back();
            filtInputs.push_back(n);
            // Push the filter onto the working list.
            workingNetnodeList.push_back(filt);
            workingNet->AddNode(filt);
        }
1328 1329
        debug4 << "NetworkManager::MakePlot(): Added filter after expression "
               << "evaluator." << endl;
1330 1331
    }

hrchilds's avatar
hrchilds committed
1332
    p->SetDataExtents(dataExtents);
hrchilds's avatar
hrchilds committed
1333 1334
    workingNet->SetPlot(p);
    workingNet->GetPlot()->SetAtts(atts);
hrchilds's avatar
hrchilds committed
1335 1336
    workingNet->SetPlottype(pluginID);
    workingNet->SetPlotName(plotName);
1337

1338 1339
    debug4 << "NetworkManager::MakePlot(): Leaving NetworkManager::MakePlot()." 
           << endl;
hrchilds's avatar
hrchilds committed
1340 1341 1342 1343 1344 1345 1346 1347 1348 1349 1350 1351
}

// ****************************************************************************
//  Method:  NetworkManager::EndNetwork
//
//  Purpose:
//    Closes off the new network, assigns a pipeline specification, and a
//    network id.
//
//  Programmer:  Jeremy Meredith
//  Creation:    November  9, 2001
//
hrchilds's avatar
hrchilds committed
1352 1353 1354 1355 1356
//  Modifications:
//
//    Hank Childs, Mon Jan  5 16:04:57 PST 2004
//    Make sure that the network id is always unique for each network.
//
hrchilds's avatar
hrchilds committed
1357 1358 1359
//    Mark C. Miller, Tue Jan  4 10:23:19 PST 2005
//    Added code to set network's window ID and manage multiple VisWindow
//    objects
1360 1361 1362
//
//    Gunther H. Weber, Thu Apr 12 11:00:57 PDT 2007
//    Added missing space to debug message
1363 1364 1365 1366
//
//    Hank Childs, Mon Apr  6 13:06:22 PDT 2009
//    Add support for sending named selection info to all filters.
//
1367 1368 1369
//    Brad Whitlock, Tue Aug 10 16:11:25 PDT 2010
//    Use find() method.
//
1370 1371 1372
//    Brad Whitlock, Mon Aug 22 10:39:12 PDT 2011
//    Changed how named selections are finally set up.
//
hrchilds's avatar
hrchilds committed
1373
// ****************************************************************************
1374

hrchilds's avatar
hrchilds committed
1375
int
hrchilds's avatar
hrchilds committed
1376
NetworkManager::EndNetwork(int windowID)
hrchilds's avatar
hrchilds committed
1377
{
1378 1379
    // Checking to see if the network has been built successfully.
    if (workingNetnodeList.size() != 1)
1380
    {
1381 1382
        debug1 << "Network building still in progress.  Filter required to "
               << "absorb " << workingNetnodeList.size() << " nodes."  << endl;
1383

1384 1385
        EXCEPTION0(ImproperUseException);
    }
1386

1387 1388 1389 1390 1391
    // If a named selection has been applied, make sure all of the filters
    // know about it.
    if (!workingNet->GetSelectionName().empty())
    {
        const std::string &selName = workingNet->GetSelectionName();
1392 1393