Documentation update

Author: mabuni1998

docs/src/example_lodahl.md

@@ -9,7 +9,7 @@ In many of the examples considered so far, we only consider a single waveguide t
 ![`alt text`](./illustrations/two_waveguide_lodahl.png)
 
 
-A way to model this scenario is to have two waveguides: a waveguide to the left and the right, describing the first half of the waveguide and the latter half. For this, we use [`WaveguideBasis`](@ref) but with an extra argument specifying that we need two waveguides (see [`twowaveguide`](@ref) for an introduction). We initialize [`WaveguideBasis`](@ref) with two waveguides and a basis for the atom (note that a fockbasis with only one excitation allowed is the same as a two-level-system):
+A way to model this scenario is to have two waveguides: a waveguide to the left and the right, describing the first half of the waveguide and the latter half. For this, we use [`WaveguideBasis`](@ref) but with an extra argument specifying that we need two waveguides (see [`Multiple Waveguides`](@ref twowaveguide) for an introduction). We initialize [`WaveguideBasis`](@ref) with two waveguides and a basis for the atom (note that a fockbasis with only one excitation allowed is the same as a two-level-system):
 
 ```@example lodahl
 using WaveguideQED #hide
@@ -101,7 +101,7 @@ rcParams = PyPlot.PyDict(PyPlot.matplotlib."rcParams") #hide
 rcParams["font.size"] = 20 #hide
 rcParams["font.family"] = "serif" #hide
 rcParams["mathtext.fontset"] ="cm" #hide
-fig,axs = subplots(3,2,figsize=(6,17))
+fig,axs = subplots(3,2,figsize=(6,9))
 plot_list = [ψ2RightScat,ψ2LeftScat,ψ2LeftRightScat,ψ1RightScat,ψ1LeftScat,ψ1LeftRightScat]
 for (i,ax) in enumerate(axs)
     plot_twophoton!(ax,plot_list[i],times)

docs/src/theoreticalbackground.md

@@ -108,7 +108,7 @@ nothing #hide
 ```
 ![alt text](created_onephoton_continuous_fockstate.svg)
 
-We see a spike around `t = times[10] = 0.9`, where we now created an excitation. In itself, the waveguide basis, states, and operators are not particularly interesting, but when combined with other quantum mechanical systems such as cavities and emitters, the framework can produce powerful results. See [`combining`](@ref) for an introduction on how to combine with quantum systems defined in ['QuantumOptics.jl'](https://qojulia.org/).
+We see a spike around `t = times[10] = 0.9`, where we now created an excitation. In itself, the waveguide basis, states, and operators are not particularly interesting, but when combined with other quantum mechanical systems such as cavities and emitters, the framework can produce powerful results. See [`Combining with QuantumOptics`](@ref combining) for an introduction on how to combine with quantum systems defined in ['QuantumOptics.jl'](https://qojulia.org/).
 
 
 ## Continuous two-photon fock states
@@ -172,7 +172,7 @@ If we want to create a two-photon Gaussian state, we instead do:
 nothing #hide
 ```
 
-Here, we defined the two-photon equivalent of our single-photon Gaussian state. When we visualize it, we now need two times, and we make a contour plot. This is easily done by viewing the two-photon state and using [`plot_twophoton`](@ref): 
+Here, we defined the two-photon equivalent of our single-photon Gaussian state. When we visualize it, we now need two times, and we make a contour plot. This is easily done by viewing the two-photon state and using [`plot_twophoton!`](@ref): 
 
 ```@example theory
 viewed_state = TwoPhotonView(ψ)

src/WaveguideInteraction.jl

@@ -257,8 +257,7 @@ const TensorWaveguideInteraction = LazyTensor{B,B,F,V,T} where {B,F,V,T<:Tuple{W
 
 function expect(a::T,psi::Ket) where T<:Union{WaveguideInteraction,TensorWaveguideInteraction}
     out = 0
-    tmp_data = QuantumOpticsBase._tp_matmul_get_tmp(eltype(psi.data), (length(psi.data),), :tp_ex, psi.data)
-    tmp_ket = Ket(psi.basis,tmp_data)
+    tmp_ket = Ket(psi.basis)
     for i in 1:get_nsteps(basis(a))
         set_waveguidetimeindex!(a,i)
         mul!(tmp_ket,a,psi,1,0)

src/WaveguideOperator.jl

@@ -70,7 +70,7 @@ Annihilation operator ``w`` for [`WaveguideBasis`](@ref) ``w_i(t_k) | 1_j \\empt
 - `i` determines which waveguide the operator acts on and should be `i ≤ Nw`. If `Nw=1` then `i=1` is assumed (there is only one waveguide).
 
 # Returns
-[`WaveguideDestroy{B,B,Np,i} <: WaveguideOperator`](@ref)
+[`WaveguideDestroy`](@ref)
 """
 function destroy(basis::WaveguideBasis{Np,1}) where {Np}
     B = typeof(basis)
@@ -95,7 +95,7 @@ Creation operator ``w^\\dagger`` for [`WaveguideBasis`](@ref) ``w_i(t_k)^\\dagge
 - `i` determines which waveguide the operator acts on and should be `i ≤ Nw`. If `Nw=1` then `i=1` is assumed (there is only one waveguide).
 
 # Returns
-[`WaveguideCreate{B,B,Np,i} <: WaveguideOperator`](@ref)
+[`WaveguideCreate`](@ref)
 """
 function create(basis::WaveguideBasis{Np,1}) where {Np}
     B = typeof(basis)