Check for ac

lcapy/analysis.py

@@ -13,6 +13,8 @@ def __init__(self, cct, components):
         self.has_ic = False
         self.zeroic = True
         self.has_s = False
+        self.has_ac = False
+        self.has_dc = False
         self.has_transient = False
         self.ac_count = 0
         self.dc_count = 0
@@ -46,6 +48,10 @@ def __init__(self, cct, components):
                     self.dc_count += 1
                 if not elt.is_causal:
                     self.causal = False
+                if elt.has_ac:
+                    self.has_ac = True
+                if elt.has_dc:
+                    self.has_dc = True
             if elt.is_dependent_source:
                 self.dependent_sources.append(eltname)
             if elt.is_reactive:
@@ -56,9 +62,6 @@ def __init__(self, cct, components):
 
         num_sources = len(self.independent_sources)
 
-        self.has_dc = self.dc_count > 0
-        self.has_ac = self.ac_count > 0
-
         self.ac = self.ac_count > 0 and \
             (num_sources == self.ac_count) and not self.has_ic
         self.dc = self.dc_count > 0 and \

lcapy/netlist.py

@@ -655,16 +655,17 @@ def laplace(self):
 
         """
 
+        lap =  self.select('laplace')
         if self.is_causal:
-            return self.select('laplace')
+            return lap
 
         # FIXME
         if self.has_ac:
             warn('Ignoring initial conditions due to AC source')
+            return
 
         # Calculate initial conditions due to DC source
         dc = self.dc()
-        lap = self.select('laplace')
 
         new = self._new()