SizeType dofsFromModelName(const ModelType& model) { const PsimagLite::String& modelName = model.params().model; SizeType site = 0; // FIXME : account for Hilbert spaces changing with site SizeType dofs = SizeType(log(model.hilbertSize(site))/log(2.0)); std::cerr<<"DOFS= "<<dofs<<" <------------------------------------\n"; if (modelName.find("FeAsBasedSc")!=PsimagLite::String::npos) return dofs; if (modelName.find("FeAsBasedScExtended")!=PsimagLite::String::npos) return dofs; if (modelName.find("HubbardOneBand")!=PsimagLite::String::npos) return dofs; // max number here, site dependence taken into account elsewhere if (modelName.find("Immm")!=PsimagLite::String::npos) return 4; return 0; }
bool observeOneFullSweep( IoInputType& io, const GeometryType& geometry, const ModelType& model, const PsimagLite::String& obsOptions, bool hasTimeEvolution) { bool verbose = false; typedef typename SparseMatrixType::value_type FieldType; typedef Observer<FieldType,VectorWithOffsetType,ModelType,IoInputType> ObserverType; typedef ObservableLibrary<ObserverType,TargettingType> ObservableLibraryType; SizeType n = geometry.numberOfSites(); //PsimagLite::String sSweeps = "sweeps="; //PsimagLite::String::SizeTypeype begin = obsOptions.find(sSweeps); //if (begin != PsimagLite::String::npos) { // PsimagLite::String sTmp = obsOptions.substr(begin+sSweeps.length(),PsimagLite::String::npos); //std::cout<<"sTmp="<<sTmp<<"\n"; // n = atoi(sTmp.c_str()); //} ObservableLibraryType observerLib(io,n,hasTimeEvolution,model,verbose); bool ot = false; if (obsOptions.find("ot")!=PsimagLite::String::npos || obsOptions.find("time")!=PsimagLite::String::npos) ot = true; if (hasTimeEvolution && ot) { observerLib.measureTime("superDensity"); observerLib.measureTime("nupNdown"); observerLib.measureTime("nup+ndown"); if (obsOptions.find("sz")!=PsimagLite::String::npos) observerLib.measureTime("sz"); } if (hasTimeEvolution) observerLib.setBrackets("time","time"); const PsimagLite::String& modelName = model.params().model; SizeType rows = n; // could be n/2 if there's enough symmetry // Immm supports only onepoint: if (modelName=="Immm" && obsOptions!="onepoint") { PsimagLite::String str(__FILE__); str += " " + ttos(__LINE__) + "\n"; str += "Model Immm only supports onepoint\n"; throw PsimagLite::RuntimeError(str.c_str()); } SizeType numberOfDofs = dofsFromModelName(model); if (!hasTimeEvolution && obsOptions.find("onepoint")!=PsimagLite::String::npos) { observerLib.measureTheOnePoints(numberOfDofs); } if (modelName.find("Heisenberg")==PsimagLite::String::npos) { if (obsOptions.find("cc")!=PsimagLite::String::npos) { observerLib.measure("cc",rows,n); } if (obsOptions.find("nn")!=PsimagLite::String::npos) { observerLib.measure("nn",rows,n); } } if (obsOptions.find("szsz")!=PsimagLite::String::npos) { observerLib.measure("szsz",rows,n); } if (modelName.find("FeAsBasedSc")!=PsimagLite::String::npos || modelName.find("FeAsBasedScExtended")!=PsimagLite::String::npos || modelName.find("HubbardOneBand")!=PsimagLite::String::npos) { bool dd4 = (obsOptions.find("dd4")!=PsimagLite::String::npos); if (obsOptions.find("dd")!=PsimagLite::String::npos && !dd4) { // && //geometry.label(0).find("ladder")!=PsimagLite::String::npos) { observerLib.measure("dd",rows,n); } // FOUR-POINT DELTA-DELTA^DAGGER: if (dd4 && geometry.label(0).find("ladder")!=PsimagLite::String::npos) { observerLib.measure("dd4",rows,n); } // if dd4 } if (modelName.find("HubbardOneBand")!=PsimagLite::String::npos && obsOptions.find("multi")!=PsimagLite::String::npos) { observerLib.measure("multi",rows,n); } if (obsOptions.find("s+s-")!=PsimagLite::String::npos) { observerLib.measure("s+s-",rows,n); } if (obsOptions.find("ss")!=PsimagLite::String::npos) { observerLib.measure("ss",rows,n); } return observerLib.endOfData(); }